No Title    |POLVER05_1 |2
f(Ca) = v*Cao-v*Ca+ra*V
f(Cb) = v*Cbo-v*Cb+rb*V
f(Cc) = -v*Cc+rc*V
f(Cd) = -v*Cd+rd*V
r1a = -k1a*Ca*Cb^2
r2c = -k2c*Ca^2*Cc^3
ra = r1a+r2a
rb = r1b
rc = r1c+r2c
rd = r2d
r1b = 2*r1a
r2a = 2/3*r2c
r1c = -r1a
r2d = -1/3*r2c
Cao = 2
Cbo = 2
k1a = 10
k2c = 15
V = 2500
v = 100
Scd = Cc/(Cd+0.001)
Ca(0)=2
Cb(0)=1
Cc(0)=1
Cd(0)=0.1