Scilab 簡単なシミュレーション （修正案）

下の条件だと、
20. 108.6 72.4
で安定する。

-----------------------
プロット（案）
t=0 ; G(t,1)=C0(1,1) ; G(t,2)=C0(1,2) ; G(t,3)=C0(1,3)
t=t+1 ; G(t,1)=C0(t,1) ; G(t,2)=C0(t,2) ; G(t,3)=C0(t,3)
-----------------------

'初期値
Mmax=0.3 ; K=10 ; 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=10 ; 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 ;