Edinburgh/SBCode
From 2007.igem.org
< Edinburgh(Difference between revisions)
Luca.Gerosa (Talk | contribs) |
Luca.Gerosa (Talk | contribs) |
||
(2 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
+ | This is the [http://www.sbtoolbox.org/ SBToolBox] code we developed for simulation of the ODEs model. | ||
+ | |||
+ | == The SBToolBox Code == | ||
+ | |||
********** MODEL NAME | ********** MODEL NAME | ||
Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | ||
Line 20: | Line 24: | ||
********** MODEL MATLAB FUNCTIONS | ********** MODEL MATLAB FUNCTIONS | ||
function [result] = ForwardPhase(t) | function [result] = ForwardPhase(t) | ||
- | FlippingTime= | + | FlippingTime=15; |
- | DivisionTime= | + | DivisionTime=80; |
temp=0; | temp=0; | ||
temp = mod(t,2*DivisionTime+2*FlippingTime); | temp = mod(t,2*DivisionTime+2*FlippingTime); | ||
Line 31: | Line 35: | ||
function [result] = BackwardPhase(t) | function [result] = BackwardPhase(t) | ||
- | FlippingTime= | + | FlippingTime=15; |
- | DivisionTime= | + | DivisionTime=80; |
temp=0; | temp=0; | ||
temp = mod(t,2*DivisionTime+2*FlippingTime); | temp = mod(t,2*DivisionTime+2*FlippingTime); |
Latest revision as of 15:34, 24 October 2007
This is the [http://www.sbtoolbox.org/ SBToolBox] code we developed for simulation of the ODEs model.
The SBToolBox Code
********** MODEL NAME Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team ********** MODEL NOTES The model represents the Division PoPper construct with the addition of a reporter protein placed downstream ********** MODEL STATES d/dt(YFP) = YFPexp*ForwardPhase(time)-YFPdeg*YFP d/dt(GFP) = GFPexp*BackwardPhase(time)-GFPdeg*GFP YFP(0) = 0 GFP(0) = 0 ********** MODEL PARAMETERS YFPexp = 0.0001 GFPexp = 0.0001 GFPdeg = 0.05 YFPdeg = 0.05 ********** MODEL VARIABLES ********** MODEL REACTIONS ********** MODEL FUNCTIONS ********** MODEL EVENTS
********** MODEL MATLAB FUNCTIONS function [result] = ForwardPhase(t) FlippingTime=15; DivisionTime=80; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp < (DivisionTime) result=1; else result=0; end return function [result] = BackwardPhase(t) FlippingTime=15; DivisionTime=80; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp > (DivisionTime+FlippingTime) if temp < (2*DivisionTime+FlippingTime) result=1; else result=0; end else result=0; end return