bioreactor_9=function(t,y) { # # Assign state variables: xylose =y[1]; xylitol =y[2]; xylulose =y[3]; acetaldehyde=y[4]; ethanol =y[5]; acetate =y[6]; glycerol =y[7]; # # Compute fluxes J=rep(0,n) J[1]=k1*xylose; J[2]=k2*xylitol-km2*xylulose*ethanol; J[3]=k3*xylulose-km3*acetaldehyde*ethanol; J[4]=k4*acetaldehyde; J[5]=k5*acetaldehyde; J[6]=k6*xylulose; # # Time derivatives f=rep(0,n) f[1]=-J[1]; f[2]=J[1]-J[2]; f[3]=J[2]-2*J[3]-2*J[6]; f[4]=3*J[3]-J[4]-J[5]; f[5]=J[4]; f[6]=J[5]; f[7]=3*J[6]; # # Calls to bioreactor_9 ncall <<- ncall+1 # # Return derivative vector return(c(f)) }