MATLAB: I have been asked to Stimulate V-S plot of enzyme activity from the equation :- dS/dt =-Vm*S/(Km+S),

bhavya

this is my program function dy = react(t,y) dy = zeros(3,1); n=1; Ki=0.1; Vs=5; S= y(1); dS = [-Vs*(S)^n/((Ki)^n + (S)^n)]; dy = [dS]; and in command window, i wrote [T,Y]=ode45('react',[0 15],[0.2]) i am getting this T =
         0
    0.0030
    0.0060
    0.0090
    0.0121
    0.0225
    0.0330
    0.0434
    0.0539
    0.0593.......till 15.000
 i have to plot dS/dt v/s S. please help

Best Answer

Use
function main
n=1; 
Ki=0.1; 
Vs=5;
[T,Y]=ode45(@(t,y)react(t,y,n,Ki,Vs),[0 15],[0.2]);
DY(:,1)=-Vs*Y(:,1).^n./(Ki^n + Y(:,1).^n);
plot(Y(:,1),DY(:,1));
function dy = react(t,y,n,Ki,Vs) 
S= y(1); 
dS = -Vs*S^n/(Ki^n + S^n); 
dy = dS;
in a single file.
Best wishes
Torsten.
Related Question