Edinburgh/SBCode
From 2007.igem.org
(Difference between revisions)
Luca.Gerosa (Talk | contribs) |
Luca.Gerosa (Talk | contribs) |
||
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 of the Division PoPer Construct of Edinburgh iGEM 2007 Team | Model of the Division PoPer Construct of Edinburgh iGEM 2007 Team | ||
********** MODEL NOTES | ********** MODEL NOTES |
Revision as of 14:20, 24 October 2007
This is the 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=10; DivisionTime=40; temp=0; temp = mod(t,2*DivisionTime+2*FlippingTime); if temp < (DivisionTime) result=1; else result=0; end return function [result] = BackwardPhase(t) FlippingTime=10; DivisionTime=40; 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