<Text-field style="Heading 1" layout="Heading 1">Tanksystem</Text-field>

with(LinearAlgebra): with(plots):

Der er givet et system best\303\245ende af to omr\303\270rte tanke. x(t) og y(t) er koncentrationerne af et bestemt stof A i de to tanke. LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUklbXN1YkdGJDYlLUkjbWlHRiQ2JVEiRkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnLUkjbW9HRiQ2LVEiLEYnRj4vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRjQvJSlzdHJldGNoeUdGSS8lKnN5bW1ldHJpY0dGSS8lKGxhcmdlb3BHRkkvJS5tb3ZhYmxlbGltaXRzR0ZJLyUnYWNjZW50R0ZJLyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdRLDAuMzMzMzMzM2VtRictRiw2JUYuLUYjNiMtRjs2JFEiMkYnRj5GQA== er volumenstr\303\270mningshastigheder. LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUklbXN1YkdGJDYlLUkjbWlHRiQ2JVEiVkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnLUkjbW9HRiQ2LVEiLEYnRj4vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRjQvJSlzdHJldGNoeUdGSS8lKnN5bW1ldHJpY0dGSS8lKGxhcmdlb3BHRkkvJS5tb3ZhYmxlbGltaXRzR0ZJLyUnYWNjZW50R0ZJLyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdRLDAuMzMzMzMzM2VtRictRiw2JUYuLUYjNiMtRjs2JFEiMkYnRj5GQA== er rumfangene af de to tanke. Der er intet stof A i den til tank 1 indkommende str\303\270m.

LSUnQ0FOVkFTRzYnLSUpQk9VTkRTX1hHNiMsJC0lJkZsb2F0RzYjJSlpbmZpbml0eUchIiItJSlCT1VORFNfWUdGJy0lLUJPVU5EU19XSURUSEc2IyQiJWd6Ri0tJS5CT1VORFNfSEVJR0hURzYjJCIlZ2ZGLS0lKUNISUxEUkVORzZPLSUnU1RST0tFRzYlLSUoT1VUTElORUc2Iy0lJkNPTE9SRzYmJSRSR0JHJCIiIUYtRkdGRy0lL0xJTkVfVEhJQ0tORVNTRzYjJCIjNUYtLSUlREFUQUc2YnI3JCQiJCF5Ri0kIiUhKj5GLTckRlIkIiVJP0YtNyRGUiQiJXE/Ri03JEZSJCIlP0BGLTckRlIkIiVnQEYtNyQkIiRxKEYtJCIlP0FGLTckRl1vJCIlZ0FGLTckRl1vJCIlU0JGLTckRl1vJCIlIVEjRi03JEZdbyQiJT9DRi03JEZdbyQiJSFbI0YtNyRGXW8kIiVJREYtNyRGXW8kIiVxREYtNyRGXW8kIiU1RUYtNyRGXW8kIiVdRUYtNyRGXW8kIiUrRkYtNyRGXW8kIiVdRkYtNyRGXW8kIiUrR0YtNyRGXW8kIiVdR0YtNyRGXW8kIiU1SEYtNyRGXW8kIiVnSEYtNyQkIiRnKEYtJCIlISpIRi03JEZfciQiJUlJRi03JEZfciQiJXFJRi03JCQiJF0oRi0kIiU/SkYtNyRGanIkIiVnSkYtNyRGanIkIiUrS0YtNyRGanIkIiVTS0YtNyRGanIkIiUhSCRGLTckRmpyJCIlU0xGLTckJCIkUyhGLSQiJXFMRi03JEZedCQiJUlNRi03JEZedCQiJSFbJEYtNyRGXnQkIiVJTkYtNyRGXnQkIiUhZSRGLTckRl50JCIlP09GLTckRl50JCIlZ09GLTckRl1vJCIlXU9GLTckJCIkNSlGLSQiJUlPRi03JCQiJFMpRi1GX3U3JCQiJCsqRi0kIiU1T0YtNyQkIiRdKkYtRmJ2NyQkIiUrNUYtRmJ2NyQkIiVTNUYtRmJ2NyQkIiU1NkYtRmJ2NyQkIiUhPSJGLUZidjckJCIlPzdGLUZfdTckJCIlcTdGLUZfdTckJCIlKzhGLUZqdTckJCIlUzhGLUZqdTckJCIlIVEiRi1GanU3JCQiJT85Ri0kIiVTT0YtNyQkIiVxOUYtRmV4NyQkIiU1OkYtRmV4NyQkIiVdOkYtRmV4NyQkIiUrO0YtRmV4NyQkIiVdO0YtRmV4NyQkIiUhcCJGLUZleDckJCIlXTxGLUZleDckJCIlNT1GLUZleDckJCIlcT1GLUZleDckJCIlNT5GLUZleDckJCIlcT5GLUZleDckJCIlUz9GLUZleDckJCIlNUBGLUZleDckJCIlU0BGLUZldTckJCIlIT4jRi1GZXU3JCQiJUlBRi1GZXU3JCQiJXFBRi1GZXU3JCQiJT9CRi1GZXU3JCQiJSFSI0YtRmV1NyQkIiVJQ0YtRmV1NyRGXnBGZXU3JCQiJV1ERi1GZXU3JCQiJSFmI0YtRmV1NyQkIiVJRUYtRmV1NyQkIiVxRUYtRmV1NyRGXXFGZXg3JCQiJXFGRi1GZXg3JCQiJT9HRi1GZXg3JCQiJWdHRi1GZXg3JEZoXWwkIiUrT0YtNyQkIiVxR0YtJCIlXU5GLTckRl5ebCQiJTVORi03JCQiJSEpR0YtJCIlcU1GLTckJCIlISpHRi0kIiVTTUYtNyRGW19sJCIlK01GLTckJCIlK0hGLUZgdDckRmlxJCIlSUxGLTckJCIlP0hGLSQiJStMRi03JCQiJUlIRi0kIiVdS0YtNyQkIiVdSEYtJCIlNUtGLTckRlxyRl9zNyQkIiVxSEYtJCIlNUpGLTckRmFyJCIlZ0lGLTckJCIlK0lGLSQiJT9JRi03JCQiJTVJRi0kIiUhKUhGLTckRmNhbEZpX2w3JEZjYWxGZl5sNyRGY2FsJCIlU0dGLTckRmNhbCQiJSF5I0YtNyRGY2FsJCIlU0ZGLTckRmNhbCQiJSFwI0YtNyRGZmFsJCIlU0VGLTckRmZhbCQiJStFRi03JEZmYWxGZVxsNyRGZmFsRl5wNyRGZmFsRltwNyRGZmFsJCIlZ0JGLTckRmZhbCQiJTVCRi03JEZmYWxGaFtsNyRGZmFsRmVbbDckRmNhbEZiW2w3JEZjYWxGX1tsNyRGY2FsJCIlK0BGLTckRmNhbCQiJWc/Ri03JEZjYWwkIiU/P0YtNyRGY2FsRlctRj42JUZARkktRk82SzckRl1vRmhbbDckJCIkPylGLSQiJSFHI0YtNyQkIiRnKUYtRl1lbDckJCIkSSpGLUZdZWw3JCQiJCEpKkYtRl1lbDckJCIlSTVGLUZdZWw3JCQiJSEzIkYtRmhbbDckJCIlUzZGLUZoW2w3JCQiJSE+IkYtRmJvNyQkIiVJOEYtJCIlXUFGLTckRmh4JCIlU0FGLTckJCIlNTtGLUZqZmw3JEZqeUZlW2w3JCQiJSEqPUYtRl9vNyQkIiVdPkYtJCIlNUFGLTckJCIlKz9GLUZmZ2w3JEZpeiQiJStBRi03JCQiJSEzI0YtRlxobDckJCIlSUBGLUZcaGw3JCQiJXFARi1GXGhsNyRGZVtsRlxobDckJCIlIUgjRi1GXGhsNyQkIiVdQkYtRlxobDckJCIlK0NGLUZcaGw3JCQiJV1DRi1GXGhsNyRGXnBGYltsNyQkIiU/REYtRmJbbDckRmRwRlxobDckRmdwRlxobDckRmpwRmZnbDckRl1xRl9vNyRGY2JsRmVbbDckRmBibEZqZmw3JCQiJTVHRi1GZ2ZsNyRGZnFGYm83JEZmXmxGaFtsNyRGaV9sRmhbbDckRmNgbEZdZWw3JEZhckZdZWw3JEZjYWxGaWhsRmZqbC1GPjYlRkBGSS1GTzZlcTckJCIlNVdGLSQiJV0/Ri03JCQiJStXRi1GX2hsNyQkIiUhUiVGLUZnbjckJCIlIVElRi1GZWhsNyRGZ1ttRmVbbDckJCIlcVZGLSQiJXFCRi03JEZbXG1GYVxsNyRGW1xtJCIlIVwjRi03JCQiJWdWRi0kIiVTREYtNyRGZFxtJCIlIWUjRi03JEZkXG1GW11sNyRGZFxtRmZibDckRmRcbSQiJUlGRi03JEZkXG1GYl1sNyRGZFxtRmVdbDckRmRcbUZmXmw3JEZkXG1GaV9sNyRGZFxtRlxyNyRGZFxtRmNhbDckRmRcbUZeYWw3JEZbXG0kIiUrSkYtNyRGW1xtRl9zNyRGW1xtJCIlSUtGLTckRltcbSQiJXFLRi03JEZbXG0kIiU1TEYtNyRGW1xtJCIlXUxGLTckRltcbUZgX2w3JEZbXG0kIiVdTUYtNyRGW1xtJCIlIVwkRi03JEZbXG1GaXQ3JEZnW20kIiVnTkYtNyQkIiU/V0YtRmBfbTckJCIlZ1dGLSQiJXFORi03JCQiJTVYRi1GaF9tNyQkIiVxWEYtRlx1NyQkIiU/WUYtRlx1NyQkIiVxWUYtRlx1NyQkIiU/WkYtRlx1NyQkIiVxWkYtRlx1NyQkIiU/W0YtRlx1NyQkIiVxW0YtRlx1NyQkIiU/XEYtRlx1NyQkIiUhKVxGLUZcdTckJCIlP11GLUZcdTckJCIlcV1GLUZcdTckJCIlP15GLUZcdTckJCIlZ15GLUZcdTckJCIlK19GLUZcdTckJCIlU19GLUZcdTckJCIlK2BGLUZcdTckJCIlcWBGLUZcdTckJCIlU2FGLUZcdTckJCIlIVsmRi1GXHU3JCQiJV1iRi1GXHU3JCQiJStjRi1GXHU3JCQiJVNjRi1GXHU3JCQiJSFvJkYtRlx1NyQkIiU/ZEYtJCIlIWYkRi03JCQiJWdkRi1GZWRtNyQkIiUrZUYtRltebDckJCIlU2VGLUZbXmw3JCQiJSEpZUYtRltebDckJCIlU2ZGLUZidjckJCIlISpmRi1GYnY3JCQiJUlnRi1GX3U3JCQiJXFnRi1GX3U3JCQiJTVoRi1GanU3JCQiJV1oRi1GanU3JCQiJSE+J0YtRmV4NyQkIiVJaUYtRmV4NyQkIiVxaUYtRmV4NyQkIiU1akYtRmV4NyQkIiVdakYtRmV4NyQkIiUra0YtRmV4NyQkIiU/a0YtRltebDckRmhnbSQiJVNORi03JEZoZ21GZnQ3JCQiJUlrRi0kIiU1TUYtNyRGX2htRmVebTckJCIlU2tGLUZbYGw3JCQiJV1rRi1GZXM3JEZoaG0kIiVxSkYtNyQkIiVna0YtJCIlSUpGLTckRl5pbSQiJSE0JEYtNyQkIiVxa0YtJCIlXUlGLTckRmZpbUZhYWw3JEZmaW1GXHI3JEZmaW1GaV9sNyQkIiUhWydGLUZeXmw3JEZeam1GY3E3JEZeam0kIiVnRkYtNyRGXmptRmZibDckRl5qbSQiJT9FRi03JEZeam1GaVxtNyRGXmptRmZpbDckRmZpbSQiJWdDRi03JEZmaW0kIiU1Q0YtNyRGZmltRl1cbTckRl5pbUZlY2w3JEZeaW1GZ2ZsNyRGaGhtRmJbbDckRmVobSQiJV1ARi03JEZlaG1GXFtsNyQkIiU1a0YtJCIlITQjRi1GaFtuLUY+NiVGQEZJLUZPNjo3JEZnW21GaG83JCQiJV1XRi1GaG83JCQiJStYRi1GaG83JCQiJWdZRi1GaG83JCQiJWdbRi1GaG83JCQiJWddRi1GXVxtNyQkIiVnX0YtRl1cbTckJCIlP2FGLUZdXG03JCQiJWdiRi1GXVxtNyQkIiUrZEYtRmJjbDckRl5lbUZiY2w3JCQiJSEpZkYtRmJjbDckJCIlU2dGLUZiY2w3JCQiJSEzJ0YtRmJjbDckJCIlP2hGLUZiY2w3JCQiJXFoRi1GYmNsNyQkIiU/aUYtRmJjbDckJCIlIUcnRi1GYmNsNyQkIiVTakYtRmJjbDckJCIlIVInRi1GYmNsNyRGaGdtRl1cbUZgW243JCQiJT9sRi1GaG83JCQiJVNsRi1GaG8tRj42JUZARkktRk82ODckRmFyRmFwNyQkIiVTSUYtRmFwNyQkIiUhMyRGLUZhcDckJCIlP0tGLUZhcDckJCIlZ0xGLUZhcDckRltobUZhcDckJCIlK1BGLUZhcDckJCIlZ1FGLUZhcDckJCIlK1NGLUZhcDckJCIlZ1NGLUZhcDckJCIlNVRGLUZhcDckJCIlXVRGLUZhcDckJCIlIT4lRi1GYXA3JCQiJVNVRi1GYXA3JCQiJSFHJUYtRmFwNyQkIiVJVkYtRmFwNyRGW1xtRmFwNyRGZFxtRmFcbTckJCIlU1ZGLUZbW243JCQiJTVWRi0kIiVTQ0YtNyQkIiVxVUYtRl5bbjckJCIlZ1VGLUZeXGwtRj42JUZARkktRk82KDckRltcbUZmXG03JCQiJV1WRi1GZHA3JEZeY25GXGNsNyQkIiUhSCVGLUZmam03JEZbZG5GanA3JCQiJV1VRi1GanAtRj42JUZARkktRk82STckRmRbbUZiczckRl5jbkZiczckRmpkbkZiczckRmhibkZiczckJCIlIT0lRi1GYnM3JCQiJUlURi1GYnM3JCQiJSEzJUYtRmJzNyQkIiU/U0YtRmJzNyQkIiUhKVJGLUZiczckJCIlSVJGLUZiczckJCIlISpRRi1GYnM3JCQiJV1RRi1GYnM3JCQiJStRRi1GYnM3JCQiJV1QRi0kIiUhPiRGLTckJCIlNVBGLUZmZ243JCQiJXFPRi1GZmduNyRGanVGZmduNyRGaF9tRmZnbjckRmNebEZmZ243JEZoXmxGZmduNyRGY3RGZmduNyQkIiUhUiRGLSQiJSE9JEYtNyRGZV5tRmZobjckRltgbEZmaG43JEZgYGxGZmhuNyRGZWBsRmZobjckRl9zRmZobjckRlthbEZmaG43JEZnckZmaG43JEZkckZmaG43JEZpYGxGZmduNyRGZmFsRmZobjckRmZgbkZfczckRmNpbSQiJV1KRi03JCQiJVNKRi1GZ2luNyRGW2ltRmBpbTckRmZobkZoXW03JEZgaW1GaF1tRltqbi1GPjYlRkBGSS1GTzYqNyRGYXJGZmduNyRGY2FsRlxhbjckRmZgbkZgYGw3JEZpYG5GX15tNyRGW2FsJCIlIUckRi03JEZkaW5GYl5tNyRGZmduRmZfbDckRmVgbCQiJT9MRi0tRj42JUZARkktRk82NTckRl5pbSQiJUlHRi03JEZbYG5GYVtvNyQkIiVnbEYtRmFbbzckJCIlK21GLUZhW283JCQiJV1tRi1GYVtvNyQkIiU1bkYtRmFbbzckJCIlcW5GLUZhW283JCQiJT9vRi1GYVtvNyQkIiVncEYtRmFbbzckJCIlK3JGLUZlXWw3JCQiJVNzRi1GZV1sNyQkIiUhRyhGLUZlXWw3JCQiJTV0Ri1GX2psNyQkIiUhSChGLUZgYmw3JCQiJV1zRi1GYHE3JCQiJT9zRi0kIiU/RkYtNyQkIiUhPihGLUZdcTckJCIlcXJGLUZeXWxGY15vLUY+NiVGQEZJLUZPNik3JCQiJSt0Ri1GZV1sNyRGZl1vRmZxNyQkIiVnc0YtRmZebDckRl1db0ZpcTckJCIlNXNGLSQiJVNIRi03JCQiJSE9KEYtRmlgbDckRmReb0ZpYGwtRj42JUZARkktRk82NjckJCIkXSJGLUZgcTckJCIkIT5GLUZgcTckJCIkUyNGLUZgcTckJCIkIUdGLUZgcTckJCIkUyRGLUZgcTckJCIkKyVGLUZgcTckJCIkXSVGLUZgcTckJCIkKyZGLUZgcTckJCIkUyZGLUZgcTckJCIkIWVGLUZgcTckJCIkPydGLUZgcTckJCIkXSdGLUZiam03JCQiJCsoRi1GYmptNyQkIiRJKEYtRmJdbDckRl1vRmJdbDckRmdib0ZjYmw3JCQiJCFwRi1GXl5vNyRGYWJvRl1xNyRGXmJvJCIlIW8jRi03JCQiJDUnRi0kIiVnRUYtLUY+NiVGQEZJLUZPNio3JEZfckZiam03JEZqckZfamw3JEZnYm9GXWJsNyQkIiQ1KEYtRl5ebDckJCIkIW9GLUZjX2w3JEZhYm9GaXE3JEZeYm9GXmBsNyRGY2NvRmVfby1GPjYlRkBGSS1GTzZnbzckRmJkbyQiJWc4Ri03JEZcY28kIiUhUiJGLTckRmRib0ZjeDckJCIkPyhGLSQiJSFbIkYtNyRGZ2JvRlt5NyRGXnQkIiVTOkYtNyRGX3IkIiUhZiJGLTckRl1vJCIlPztGLTckRlJGZHk3JCQiJCspRi0kIiUhbyJGLTckRmh1JCIlNTxGLTckJCIkSSlGLSQiJVM8Ri03JCQiJF0pRi0kIiUheiJGLTckJCIkISkpRi0kIiVTPUYtNyQkIiQhKilGLUZgejckJCIkNSpGLUZjejckJCIkPypGLSQiJVM+Ri03JEZjZWxGZno3JEZldkZpZ2w3JCQiJGcqRi1GVzckRmZlbEZfZGw3JCQiJCEqKkYtRltcbjckJCIlNTVGLUZiaGw3JEZpZWwkIiUhPSNGLTckJCIlXTVGLUZmZ2w3JEZcZmxGYm83JCQiJSs2Ri1GZWNsNyQkIiVJNkYtRmJjbDckJCIlZzZGLUZfaWw3JEZhd0ZoY243JCQiJTU3Ri1GXnA3JCQiJUk3Ri1GZmlsNyQkIiVdN0YtRmVcbDckJCIlZzdGLUZpXG03JEZnd0ZncDckJCIlIUgiRi1GaWJsNyQkIiU1OEYtRl5dbDckJCIlXThGLUZdcTckRmB4Rl5ebzckRmN4RmNibDckJCIlZzlGLUZgcTckJCIlKzpGLUZiXWw3JEZpZW9GYl1sNyQkIiVxOkYtRmJqbTckJCIlIWUiRi1GXl1tNyRGXGZvRl1xNyRGXGZvRmVjbzckRmFccEZbXWw3JEZpZW9GZmptNyRGXVxwRmZqbTckRmpbcEZmam03JCQiJUk5Ri1GW11sNyQkIiUrOUYtRmpwNyQkIiVxOEYtRmVjbzckRmVbcEZmYmw3JCQiJT84Ri1GXl5vNyRGYltwRmBxNyRGX1twRmBibDckJCIlIUciRi1GX2psNyRGW1twRl1ibDckRmVqb0ZeXmw3JCQiJSs3Ri1GZl5sNyQkIiVxNkYtRltfbDckRl9mbEZmXmw3JEZbam9GZnE3JEZbam9GX2psNyRGX2ZsRmBibDckRmRecEZiam03JEZiam9GYHE3JCQiJVM3Ri1GY2JsNyRGXF5wRl5dbTckRmddcEZeXW03JEZlW3BGXl5vNyRGZWZsRmZibEZiX3AtRj42JUZARkktRk82Ym83JEZkXG0kIiU1OUYtNyRGZFxtJCIlXTlGLTckRltcbSQiJT86Ri03JEZbXG1GYVxwNyRGZFttRl1nbDckRmFbbSQiJWc7Ri03JEZcW20kIiUrPEYtNyQkIiVJV0YtJCIlSTxGLTckJCIlU1dGLSQiJXE8Ri03JEZjXG4kIiUrPUYtNyRGZl9tJCIlST1GLTckJCIlcVdGLUZgejckJCIlIVwlRi0kIiUrPkYtNyRGZlxuJCIlST5GLTckRltgbUZmejckJCIlSVhGLUZpZ2w3JCQiJV1YRi1GaXo3JEZeYG1GX2hsNyQkIiUhZSVGLUZnbjckJCIlK1lGLUZqbjckJCIlNVlGLUZiW2w3JCQiJVNZRi1GamZsNyQkIiVdWUYtRmhbbDckRmRgbUZbXGw3JCQiJSFwJUYtRl1cbTckJCIlNVpGLUZbcDckJCIlSVpGLUZbW243JCQiJV1aRi0kIiUrREYtNyQkIiVnWkYtRmFwNyQkIiUheSVGLUZkcDckJCIlK1tGLUZcY2w3JCQiJTVbRi1GW11sNyQkIiVJW0YtRmVjbzckJCIlU1tGLUZmYmw3JEZgYW1GY2JsNyQkIiUhKltGLUZiXWw3JCQiJTVcRi1GX2psNyQkIiVJXEYtRmZxNyQkIiVdXEYtRmZebDckJCIlISpcRi1GaXE3JEZpYW1GaV9sNyQkIiVdXUYtRmlxNyQkIiUhMyZGLUZjX2w3JCQiJSteRi1GXl5sNyRGY2dwRmFbbzckRlxibUZfamw3JCQiJVNdRi0kIiUheiNGLTckJCIlK11GLUZqZ3A3JCQiJWdcRi1GamdwNyRGY2ZwRmVdbDckRmBmcEZmcTckRl1mcEZmXmw3JEZgYW1GaXE3JEZpZXBGZV9vNyRGXWFtRmlgbDckRmNlcEZhYWw3JCQiJSF6JUYtRmRyNyRGamRwRmZgbjckRmdgbUZmYG43JCQiJStaRi1GZmFsNyRGXGRwRmhhbDckRl9pcEZlX283JEZnYG1GaXE3JEZlZHBGW19sNyRGXWVwRmZebDckRmNlcEZeXmw3JEZpZXBGaF1sNyRGYGFtRmZxNyQkIiUrXEYtRl1ibDckJCIlISlbRi1GYVtvLUY+NiVGQEZJLUZPNis3JCQiJWc6Ri1GY2BsNyRGZGpwRmFhbDckRmFccEZmYG43JEZkXHBGZ3I3JEZcZm9GaF1tNyRGYXlGYGltNyQkIiVJO0YtRl9zNyQkIiVTO0YtRmZnbjckRmNgcEZlYGwtRj42JUZARkktRk82KzckRmVlb0ZlYGw3JEZeYHBGW2ltNyRGaWVvRmdpbjckRmFccEZcczckRmF5RmhdbTckRl9bcUZpYG43JCQiJXE7Ri1GZ3I3JEZmYHBGaGltNyRGZmBwRmZgbi1GPjYlRkBGSS1GTzYsNyRGZ2dvRmFyNyQkIiU/PUYtRmNhbDckRmJnb0ZoaW03JCQiJWc8Ri1GaWBuNyRGanlGW2FsNyRGW2FwRmdpbjckRlthcEZmaG43JEZbYXBGXGFuNyRGW2FwJCIlZ0tGLTckRlthcEZfXm0tRj42JUZARkktRk82KzckRl5icEZjYWw3JEZeYnBGXmFsNyRGXmJwRmhdbTckRl5icEZnaW43JEZhZ2xGW2ltNyRGYWdsRmVgbDckRl5icEZgYGw3JEZiaG9GYGBsNyRGYmhvRmVzLUY+NiVGQEZJLUZPNiY3JCQiJWc9Ri1GYGltNyQkIiU/PkYtRlxzNyRGZGdsRlthbEZcX3EtRj42JUZARkktRk82LzckJCIlNT9GLUZhcjckRmJkbEZjYWw3JEZXRmhpbTckRml6RmlgbjckRl5bbUZbYWw3JEZfZGxGZ2luNyRGX2RsRmZobjckRlpGZWBsNyRGWkZgYGw3JEZaRmhzNyRGX2RsRmZfbDckRml6Rl9hbkZeYHEtRj42JUZARkktRk82NjckRmhibUZmYG43JCQiJUlfRi1GaWBuNyRGZWBxRlxzNyRGaGJtRmRpbjckJCIlIUcmRi1GX3M3JCQiJTVgRi1GZGluNyQkIiU/YEYtRlxzNyQkIiVJYEYtRmNpbTckJCIlU2BGLUZcczckRmZhcUZfczckRmZhcUZiczckRmZhcUZlczckRmNhcUZfXm03JEZjYXFGYl5tNyRGXWFxRlt0NyRGamBxRl9hbjckJCIlXV9GLUZgdDckJCIlNV9GLUZgdDckJCIlcV5GLUZfYW43JEZmYnFGZV5tLUY+NiVGQEZJLUZPNi43JEZkY21GYWFsNyQkIiVxYUYtRmRyNyQkIiVdYUYtRl5hbDckRmFjbUZjaW03JEZlXW5GXHM3JCQiJTVhRi1GZGluNyQkIiUrYUYtRmZobjckRmpjcUZcYW43JEZqY3FGYl1xNyRGZ2NxRmhzNyRGZV1uRmpqbkZfZHEtRj42JUZARkktRk82LDckRmdjbUZoaW03JCQiJVNiRi1GY2ltNyQkIiVJYkYtRmBpbTckJCIlP2JGLUZfczckRlxlcUZiczckJCIlNWJGLUZcXm03JCQiJStiRi1GYl1xNyRGY2VxRltgbDckRmlkcUZqam43JEZpZHFGYl5tLUY+NiVGQEZJLUZPNiY3JEZfY3FGZGluNyRGXGVxRmRpbjckRmhdbkZkaW43JCQiJSFmJkYtRmRpbi1GPjYlRkBGSS1GTzYvNyRGXWRtRmZhbDckJCIlXWNGLUZmYG43JCQiJWdjRi1GY2ltNyRGW2dxRmBpbTckJCIlcWNGLUZfczckRmBkbUZmZ243JEZgZG1GXF5tNyRGYGRtRl9ebTckRl9ncUZbYGw3JEZbZ3FGZl9sNyRGXWRtRl9hbjckJCIlNWNGLUZgdDckRmhncSQiJSFRJEYtLUY+NiVGQEZJLUZPNio3JCQiJD8iRi1GaWhsNyRGYmhxRmVvNyRGYmhxRmhvNyRGYmhxRltwNyRGYmhxRltbbjckJCIkNSJGLUZhXG03JEZpaHFGYXA3JEZpaHEkIiVnREYtLUY+NiVGQEZJLUZPNidGYWhxNyQkIiRnIkYtRmlobDckJCIkKyNGLUZpaGw3JEZmYG9GaWhsNyQkIiRdI0YtRl1lbC1GPjYlRkBGSS1GTzYmRmZocTckRmRpcUZbcDckRmdpcUZbcDckRmdpcUZeW24tRj42JUZARkktRk82KzckJCIkcSJGLUZmaWw3JCQiJDUjRi1GZmlsNyQkIiQ/I0YtRmVcbDckRl9bckZoXGw3JEZnaXFGZmptNyRGZGlxRmlibDckRmNgb0ZqcDckJCIkSSNGLUZqcDckRmZgb0ZqcC1GPjYlRkBGSS1GTzYqNyQkIiVJb0YtRl9pbDckRl5cckZoY243JEZeXHJGXnA3JEZkXG8kIiU1REYtNyRGZFxvRmVcbDckRmRcb0ZoXGw3JEZkXG9GW11sNyQkIiUrb0YtRmlibC1GPjYlRkBGSS1GTzYnNyRGXlxyRl5cbDckJCIlcW9GLUZeXGw3JCQiJStwRi1GX2lsNyQkIiVTcEYtRl9pbDckJCIlP3BGLUZfaWwtRj42JUZARkktRk82JTckRmRcb0ZmaWw3JCQiJWdvRi1GZmlsNyQkIiUhKm9GLUZmaWwtRj42JUZARkktRk82KzckRmVeckZcY2w3JEZqXXJGZ3A3JCQiJUlwRi1GaWJsNyRGXl9yRmBjbzckRmpdciQiJTVGRi03JEZlXnJGXl1tNyRGal1yRmNibDckRmdcb0ZjYmw3JEZnXXJGY2JsLUY+NiVGQEZJLUZPNik3JEZcX21GX283JEZcX21GYm83JEZmdCQiJStCRi03JEZmdEZcaWw3JEZmdEZeXGw3JEZmdEZhXGw3JEZmdEZoY24tRj42JUZARkktRk82JjckRlxfbUZiW2w3JEZpdEZcaGw3JEZoX21GXGhsNyRGW15sRlxobC1GPjYlRkBGSS1GTzYlNyQkIiUrTkYtRmVjbDckRltobUZlY2w3JEZgXmxGZWNsLUY+NiVGQEZJLUZPNiY3JEZgX21GXltuNyRGYF5sRmhjbjckRmBebEZecDckRmBebEZnZHAtRj42JUZARkktRk82JzckRmhfbUZpYGw3JEZgXmxGYWFsNyRGW2htRmhpbTckRltobUZjaW03JEZgX21GY2ltLUY+NiVGQEZJLUZPNig3JEZgX21GZV9vNyRGW15sRmVfbzckRmV4RmVfbzckJCIlIW8kRi1GZV9vNyQkIiVJUEYtRmVfbzckJCIlU1BGLUZeYGwtRj42JUZARkktRk82JjckRmhfbUZjYWw3JEZfdUZjYWw3JEZldUZkckZdZHItRj42JUZARkktRk82JTckRl91RmlgbjckRl91RlxzNyRGX3VGZGluLUY+NiVGQEZJLUZPNjg3JEZgeCQiJVNRRi03JEZgXXAkIiVxUUYtNyRGaF9wJCIlK1JGLTckRmN4RmhmbjckRl1dcEZlZm43JCQiJVM5Ri0kIiU1U0YtNyRGaltwRl9mbjckRmh4JCIlP1RGLTckRmVlb0ZiYm43JCQiJSFcIkYtRmllbjckRl1ccCQiJTVVRi03JEZpZW8kIiUrVUYtNyRGYVxwRmllbjckRl9mb0ZcZm43JEZjYHBGX2ZuNyRGZWZvJCIlXVNGLTckRmZgcEZnZXI3JCQiJT88Ri0kIiVxUkYtNyRGankkIiVTUkYtNyQkIiUheSJGLUZgZXI3JEZjYXBGXWVyNyRGY2FwRmZhbi1GPjYlRkBGSS1GTzYpNyRGaGZvRmJibjckRmZgcEZpZW43JEZneUZiZnI3JEZlZm9GaGJuNyRGXVxxRltkbjckRmNgcCQiJStWRi1GZGhyLUY+NiVGQEZJLUZPNjo3JEZlYm1GZmFuNyRGZWJtJCIlNVJGLTckRmVibSQiJV1SRi03JEZjYnEkIiUhKlJGLTckJCIlP19GLSQiJUlTRi03JEZlYHFGX2ZuNyRGaGJtRltmcjckRmBicUZiYm43JEZiXW5GaWVuNyRGamBxJCIlP1VGLTckJCIlIUgmRi1GXmVuNyRGYGFxRl9qcjckJCIlZ2BGLUZpZW43JEZnY3EkIiVTVEYtNyRGYWNtRltmcjckRmRjbSQiJSE0JUYtNyRGY2VxRltncjckRmBlcUZiZm43JEZcZXFGZWZuNyRGaWRxRmBpcjckRmZkcSQiJT9SRi03JEZnY21GW2duNyRGaF1uRmZhbkZnW3MtRj42JUZARkktRk82MDckRlxlcUZiZnI3JEZnY21GZWZyNyQkIiUhZSZGLUZlYm43JEZgZnFGX2pyNyRGamNtRl5lbjckRmpjbUZqZG43JCQiJXFiRi1GZmNuNyRGZmRxRl5jbjckRmBlcUZjY243JEZcZXFGW1xtNyRGaF1uRltcbTckRmpjbUZbXG03JEZdZG1GZFxtNyQkIiVJY0YtRmZkbg==

Massebevarelse giver for tank 1:

Diff(V[1]*x(t),t)=F[1]*y(t)-F[1]*x(t);

Massebevarelse giver for tank 2:

Diff(V[2]*y(t),t)=-F[1]*y(t)+F[1]*x(t)-F[2]*y(t);

Med LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYpLUkjbWlHRiQ2I1EhRictRiM2JS1JJW1zdWJHRiQ2JS1GLDYlUSJyRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSNtbkdGJDYkUSIxRicvRjtRJ25vcm1hbEYnLyUvc3Vic2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUSI9RidGQy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGTi8lKXN0cmV0Y2h5R0ZOLyUqc3ltbWV0cmljR0ZOLyUobGFyZ2VvcEdGTi8lLm1vdmFibGVsaW1pdHNHRk4vJSdhY2NlbnRHRk4vJSdsc3BhY2VHUSwwLjI3Nzc3NzhlbUYnLyUncnNwYWNlR0Znbi1JJm1mcmFjR0YkNigtRiM2Iy1GMjYlLUYsNiVRIkZGJ0Y3RjpGPUZFLUYjNiMtRjI2JS1GLDYlUSJWRidGN0Y6Rj1GRS8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGYHAvJSliZXZlbGxlZEdGTi1GSTYtUSIsRidGQ0ZML0ZQRjlGUUZTRlVGV0ZZL0ZmblEmMC4wZW1GJy9GaW5RLDAuMzMzMzMzM2VtRictRiM2JS1GMjYlRjQtRiM2Iy1GQDYkUSIyRidGQ0ZFRkgtRltvNihGXW8tRiM2Iy1GMjYlRmhvRmFxRkVGW3BGXnBGYXBGY3BGZXAtRiM2JS1GMjYlRjQtRiM2Iy1GQDYkUSIzRidGQ0ZFRkgtRltvNigtRiM2Iy1GMjYlRmFvRmFxRkVGaHFGW3BGXnBGYXBGY3BGKw== f\303\245s s\303\245 systemet

sys:={diff(x(t),t)=r1*y(t)-r1*x(t),diff(y(t),t)=r2*x(t)-(r2+r3)*y(t)};

Med A defineret som

A:=Matrix([[-r1,r1],[r2,-r2-r3]]);

P\303\245 vektorform kan systemet skrives som f\303\270lger (f\303\270rst kun til udstillingsbrug):

<diff(x(t),t),diff(y(t),t)> = A &. <x(t),y(t)>;

eller ganget ud:

<diff(x(t),t),diff(y(t),t)> = A.<x(t),y(t)>;

Taleksempler. r-v\303\246rdierne er valgt simple for overskuelighedens skyld:

r1:=4: r2:=1: r3:=3:

sys;

Fuldst\303\246ndig l\303\270sning ved dsolve:

dsolve(sys);

L\303\270sningen opskrevet som en vektor:

FULD1:=subs(%,<x(t),y(t)>);

"H\303\245ndkraft":

EV,P:=Eigenvectors(A);

Den fuldst\303\246ndige l\303\270sning

L:=c1*exp(EV[1]*t)&*Column(P,1)+c2*exp(EV[2]*t)&*Column(P,2):

L=eval(L,`&*`=`*`);

Det er det samme som

P.<c1*exp(EV[1]*t),c2*exp(EV[2]*t)>;

De v\303\246rdier af konstanterne c1 og c2, der svarer til l\303\270sningen, der opfylder begyndelsesbetingelsen x(0) = 8, y(0) = 2 findes ved inds\303\246ttelse af t = 0 i den fuldst\303\246ndige l\303\270sning. Dvs. at c1 og c2 skal opfylde ligningen

P.<c1,c2> = <8,2>;

Vi l\303\270ser dette ligningssystem:

LinearSolve(P,<8,2>);

Dermed er l\303\270sningen

subs(c1=%[1],c2=%[2],L): %=eval(%,`&*`=`*`);

Den l\303\270sning, der opfylder begyndelsesbetingelsen x(0) = 8, y(0) = 2 kan findes direkte i Maple s\303\245ledes:

dsolve(sys union {x(0)=8,y(0)=2});

L\303\270sningen som vektor:

XY:=subs(%,<x(t),y(t)>);

plot(XY,t=0..3,color=[red,blue],legend=["Koncentrationen i tank 1, x","Koncentrationen i tank 2, y"]);

FASEPLANSBILLEDE animeret:

animate(plot,[[XY[1],XY[2],t=0..T],labels=[x,y],thickness=3],'T'=0..3,frames=50);

Vi tegner mange l\303\270sninger svarende til forskellige startv\303\246rdier (ogs\303\245 urealistiske - negative koncentrationer!):

with(DEtools):

DEplot(sys,[x(t),y(t)],t=0..5,[seq([x(0)=0,y(0)=k],k=-5..5),seq([x(0)=k,y(0)=0],k=-5..5)],stepsize=.05,arrows=none,linecolor=red);

En animeret version:

animate(DEplot,[sys,[x(t),y(t)],t=0..T,[seq([x(0)=0,y(0)=k],k=1..5),seq([x(0)=k,y(0)=0],k=1..5)],stepsize=.05,arrows=none,linecolor=red],T=[.01,.05,.1,.2,.3,.4,.5,.7,1,2,5]);