<p fontName="Arial"></p>

87bea091-33eb-4fa8-94c8-d9bd1472a308 53 4 5 false false 0 decimal

&[DATE] &[TIME] - &[FILENAME]

A.p 0.2 :

P.atm 1 :

L 0.04 :

V 2 :

W.s 0.300 :

O.2o 21.0 :

V.m1 22.71 :

K.m1 7.63 :

K.i1 14.42 :

P.O2 0.2 :

CO.2o 0.03 :

V.m2 17.64 :

K.m2 5.08 :

K.i2 11.99 :

P.CO2 0.2 :

Ode solver setup

oco.0 21.0 0.03 stack :

t.min 0 :

t.max 80 :

steps 80 :

Ensure that rx is not assigned.

r1 r2 Clear 1

t o D.1 1 V / A.p P.O2 * P.atm * L / 0.01 O.2o o - (* (* W.s r1 * - (* :

100 replaced with 1

t co D.2 1 V / A.p P.CO2 * P.atm * L / 0.01 CO.2o co - (* (* W.s r2 * + (* :

O.2i CO.2i V.m K.m K.i r V.m O.2i * K.m 1 CO.2i K.i / + (O.2i * + (1 -^* :

oco O.2i CO.2i 12 mat oco.0 transpose : t t.min : Δt t.max t.min - (steps 1 -^* : out t O.2i CO.2i 13 mat : n 1 steps range r1 O.2i CO.2i V.m1 K.m1 K.i1 r : r2 O.2i CO.2i V.m2 K.m2 K.i2 r : O.2i O.2i t t Δt + 1 D.1 rkfixed 22 el : CO.2i CO.2i t t Δt + 1 D.2 rkfixed 22 el : t t Δt + : out out t O.2i CO.2i 13 mat stack : 61 line for out 61 line :

t x D O.2i CO.2i 12 mat x transpose : r1 O.2i CO.2i V.m1 K.m1 K.i1 r : r2 O.2i CO.2i V.m2 K.m2 K.i2 r : t O.2i D.1 t CO.2i D.2 12 mat transpose 41 line :

oco' oco.0 t.min t.max steps D Rkadapt :

O.2i oco 1 col oco 2 col augment :

O'.2i oco' 1 col oco' 2 col augment :

CO.2i oco 1 col oco 3 col augment :

CO'.2i oco' 1 col oco' 3 col augment :

O.2i CO.2i O'.2i + augment CO'.2i o augment 4 1 sys

ODE system with 4 variables

Initials

O0.2i CO0.2i 12 mat oco.0 transpose :

r10 O0.2i CO0.2i V.m1 K.m1 K.i1 r :

r20 O0.2i CO0.2i V.m2 K.m2 K.i2 r :

oco'.0 O0.2i CO0.2i r10 r20 14 mat transpose :

t x D O.2i CO.2i r1 r2 14 mat x transpose : r1 O.2i CO.2i V.m1 K.m1 K.i1 r : r2 O.2i CO.2i V.m2 K.m2 K.i2 r : D1 1 V / A.p P.O2 * P.atm * L / 0.01 O.2o O.2i - (* (* W.s r1 * - (* : D2 1 V / A.p P.CO2 * P.atm * L / 0.01 CO.2o CO.2i - (* (* W.s r2 * + (* : D1 D2 r1 r2 14 mat transpose 61 line :

Define first r1 and r2

oco' oco'.0 t.min t.max steps D Rkadapt :

O'.2i oco' 1 col oco' 2 col augment :

CO'.2i oco' 1 col oco' 3 col augment :

O.2i CO.2i O'.2i + augment CO'.2i o augment 4 1 sys <p fontName="Arial"></p>

Defining first D1 and D2 fails. This code is the same, except the order of the updates.

t x D O.2i CO.2i r1 r2 14 mat x transpose : D1 1 V / A.p P.O2 * P.atm * L / 0.01 O.2o O.2i - (* (* W.s r1 * - (* : D2 1 V / A.p P.CO2 * P.atm * L / 0.01 CO.2o CO.2i - (* (* W.s r2 * + (* : r1 O.2i CO.2i V.m1 K.m1 K.i1 r : r2 O.2i CO.2i V.m2 K.m2 K.i2 r : D1 D2 r1 r2 14 mat transpose 61 line :

oco' oco'.0 t.min t.max steps D rkfixed :

Rkadapt fails, use rkfixed instead

O'.2i oco' 1 col oco' 2 col augment :

CO'.2i oco' 1 col oco' 3 col augment :

O.2i CO.2i O'.2i + augment CO'.2i o augment 4 1 sys

Using parameters into D

t x r1 r2 D.4 O.2i CO.2i r1 r2 14 mat x transpose : D1 1 V / A.p P.O2 * P.atm * L / 0.01 O.2o O.2i - (* (* W.s r1 * - (* : D2 1 V / A.p P.CO2 * P.atm * L / 0.01 CO.2o CO.2i - (* (* W.s r2 * + (* : r1 O.2i CO.2i V.m1 K.m1 K.i1 r : r2 O.2i CO.2i V.m2 K.m2 K.i2 r : D1 D2 r1 r2 14 mat transpose 61 line :

oco' oco'.0 t.min t.max steps D.4 rkfixed :

Rkadapt fails too

O'.2i oco' 1 col oco' 2 col augment :

CO'.2i oco' 1 col oco' 3 col augment :

O.2i CO.2i O'.2i + augment CO'.2i o augment 4 1 sys