Trp.m

<pre>                                             
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   This file is part of Stanford iGEM team product, 2009               % 
%                                                                       % 
%   Copyright (C) YANG He     yanghemary@gmail.com                      % 
%   All rights reserved.                                                % 
%                                                                       % 
%   $File: SiGEM_enzyme.ido.m$     $Revision:  1$                       %   
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Editted by Chris VanLang 7/13/09


clear all;

% Constants Definition
km1=0.009;
km2=0.1;
e0=0.01;
i0=1;
s0=1;
k2=0.01;
k4=0.001;

% declare initial conditions
y0=[i0 s0 0 0];

% Timing
length=60;
timespan=0.01;
n=length/timespan;

tspan=[0,200000];

[t,y]=ode45(@IDO,tspan,y0,[]);


% Calculation
%s(1)=s0;
%i(1)=i0;
%t=timespan:timespan:length;



%for k=2:n
%    es(k)=e0*s(k-1)/(s(k-1)+km1/km2*i(k-1)+km1);
%    ds=-(es(k)-es(k-1))-k2*es(k);
%    s(k)=s(k-1)+ds*timespan;
%    ei(k)=es(k)*km1/km2*s(k-1)/i(k-1);
%    di=-(ei(k)-ei(k-1))-k4*ei(k);
%    i(k)=i(k-1)+ds*timespan;
%end;

% Graph
%plot(t,es,'r',t,ei,'b',t,s,'k',t,i,'y');
%plot(t,y(1),'r',t,y(2),'b',t,y(3),'k',t,y(4),'y');
plot(t,y)
</pre>