No Title  |POLVER05_0 |1
d(Fa)/d(V) =     ra     #
d(Fb)/d(V) =     rb     #
d(Fc)/d(V) =     rc  -Rc   #
d(Fd)/d(V) =     rd     #
Ft =     Fa+Fb+Fc+Fd     #
r1a =     -k1a*Ca*Cb^2     #
r2c =     -k2c*Ca^2*Cc^3     #
r1b =     2*r1a     #
r1c =     -r1a     #
r2a =     2/3*r2c     #
r2d =     -1/3*r2c     #
d(Fcsg) / d(V) = Rc*Cc 
Fcsg(0) = 0
ra =     r1a+r2a     #
rb =     r1b     #
rc =     r1c+r2c     #
rd =     r2d     #
Ca =   Cto*(Fa/Ft)*y   #
Cb =   Cto*(Fb/Ft)*y   #
Cc =   Cto*(Fc/Ft)*y   #
Cd =   Cto*(Fd/Ft)*y   #
k1a =   1000  #
k2c =   60000   #
Fto = 20 
alpharho = .041
kcc = 2 
Rc = kcc*Cc 
d(y) / d(V) = -alpharho/2/y*(Ft/Fto) 
y(0) = 1
Cto =   0.2   #
Scd =if(V>0.0001)then((Fc+Fcsg)/Fd)else(0)  
V(0)=0
Fa(0)=10
Fb(0)=10
Fc(0)=0
Fd(0)=0
V(f)=50