MATLAB Answers

Kurt
0

Differential equations solver for a chemostat

Asked by Kurt
on 22 Apr 2019
mu=0.72; %Constant
Ks=0.4; %Constant
Yxs=0.44; %Constant
D(0)=0.6; %D @ t=0 (range 0.6 to 0.25)
D(t)=((mu*S(t))/(Ks+S(t))); %D @ t=t
X(0)=0.13 %X @ t=0 (range 0.13 to 1.67)
X(t)=(Yxs*(S(0)-((Ks*D(t))/(mu-D(t))))); %X @ t=t
S(0)=3.7; %S @ t=0 (range 3.7 to .2)
S(t)=((D(t)*Ks)/(mu-D(t))); %S @ t=t
ode1==dX/dt = (((mu*S)/(Ks+S))-D)*X);
ode2==dS/dt = (D*S(1)-(X/Yxs)*((mu*S)/(Ks+S))-D*S);
tspan = [0 10]; %time range, cond1
y0 = [0.2 ; 3.7]; %range of S, cond2
[T, Y] = ode45(odes,tspan,y0);
plot(t,X,t,S);
xSol(t) = dsolce(ode1,cond1,ode2,cond2);
I need help solving these differnital equations for a chemstat.
I also need to show D from 0.6 to 0.24.
Thanks!

  0 Comments

Sign in to comment.

0 Answers