Trp.m
From 2009.igem.org
(Difference between revisions)
Chrisvanlang (Talk | contribs) (New page: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % This file is part of Stanford iGEM team product, 2009 ...) |
Chrisvanlang (Talk | contribs) |
||
(2 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
- | + | <pre> | |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
% This file is part of Stanford iGEM team product, 2009 % | % This file is part of Stanford iGEM team product, 2009 % | ||
Line 56: | Line 56: | ||
%plot(t,y(1),'r',t,y(2),'b',t,y(3),'k',t,y(4),'y'); | %plot(t,y(1),'r',t,y(2),'b',t,y(3),'k',t,y(4),'y'); | ||
plot(t,y) | plot(t,y) | ||
+ | </pre> |
Latest revision as of 06:19, 21 October 2009
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 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)