Scilab 簡単なシミュレーション

'初期値
Mmax=0.3 ; K=1 ; Dup=0.3
C0=[100 1 100]

'--------------

dC=[0 C0(1,2)*(Mmax*C0(1,1)/(K+C0(1,1))) 0;0 0 C0(1,2)*0.2;C0(1,3)*Dup 0 0]

fC=[[dC(3,1)-dC(1,2)] [dC(1,2)-dC(2,3)] [dC(2,3)-dC(3,1)]]

C1=C0+fC

C0=C1

'--------------

これを2行にまとめると、

Mmax=0.3 ; K=1 ; Dup=0.3; C0=[100 1 100]

dC=[0 C0(1,2)*(Mmax*C0(1,1)/(K+C0(1,1))) 0;0 0 C0(1,2)*0.2;C0(1,3)*Dup 0 0] ; fC=[[dC(3,1)-dC(1,2)] [dC(1,2)-dC(2,3)] [dC(2,3)-dC(3,1)]] ; C1=C0+fC ; C0=C1