No Title   |POLVER05_3
d(Cc)/d(t) =  (mu1*v1+mu2*v2)*Cc  #
d(Cs1)/d(t) =  -mu1*v1*Cc/Ycs1  #
d(Cs2)/d(t) =  -mu2*v2*Cc/Ycs2  #
d(e1)/d(t) =  alpha1*Cs1/(K1+Cs1)*u1 - beta1*e1 - (mu1*v1+mu2*v2)*e1  #
d(e2)/d(t) =  alpha2*Cs2/(K2+Cs2)*u2 - beta2*e2 - (mu1*v1+mu2*v2)*e2  #
mu1 =  mu1max*E1*Cs1/(K1+Cs1)  #
mu2 =  mu2max*E2*Cs2/(K2+Cs2)  #
u1 =  mu1/(mu1+mu2)  #
u2 =  mu2/(mu1+mu2)  #
v1 =  mu1/ mumax  #
v2 =  mu2/mumax  #
E1 =  e1/e1max  #
E2 =  e2/e2max  #
e1max =  alpha1/( mu1max+ beta1)  #
e2max =  alpha2/( mu2max+ beta2)  #
mu1max =  .9  #
mu2max =  .6  #
K1 =  .1  #
K2 =  .5  #
alpha1 =  .0001  #
alpha2 =  .0001  #
beta1 =  .05  #
beta2 =  .05  #
Ycs1 =  .4  #
Ycs2 =  .4  #
mumax =  if (mu1>mu2) then (mu1) else (mu2)  #
t(0)=0
Cc(0)=0.1
Cs1(0)=4
Cs2(0)=20
e1(0)=8.3e-5
e2(0)=1e-6
t(f)=10