Team:Uppsala-SwedenModeling
From 2010.igem.org
(→Results) |
(→Code) |
||
Line 20: | Line 20: | ||
== Code == | == Code == | ||
+ | |||
+ | == odeIGEM03.m == | ||
+ | |||
+ | %lambda1 function file | ||
+ | |||
+ | function dX =odeIGEM02(t,X,P) | ||
+ | |||
+ | k1=P(1); | ||
+ | k2=P(2); | ||
+ | k3=P(3); | ||
+ | k4=P(4); | ||
+ | k5=P(5); | ||
+ | k6=P(6); | ||
+ | k7=P(7); | ||
+ | k8=P(8); | ||
+ | k9=P(9); | ||
+ | k10=P(10); | ||
+ | k11=P(11); | ||
+ | |||
+ | km1=P(12); | ||
+ | km2=P(13); | ||
+ | km3=P(14); | ||
+ | km4=P(15); | ||
+ | km5=P(16); | ||
+ | |||
+ | ahl=X(1); | ||
+ | lam=X(2); | ||
+ | ci=X(3); | ||
+ | laci=X(4); | ||
+ | gfp=X(5); | ||
+ | |||
+ | n1=1; | ||
+ | n2=1; | ||
+ | n3=1; | ||
+ | n4=1; | ||
+ | n5=2; | ||
+ | |||
+ | dlam=((k1*(ahl/km1)^n1)/(1+(ahl/km1)^n1))-k2*lam; | ||
+ | dahl=0;%k10*ahl-k11*ahl; | ||
+ | dci=((k3*(lam/km3)^n2)/(1+(lam/km3)^n2))-k4*ci; | ||
+ | dlaci=(k5/(1+(ci/km3)^n3))+((k6*(lam/km4)^n4)/(1+(lam/km4)^n4))-k7*laci; | ||
+ | dgfp=(k8/(1+(laci/km5)^n5))-k9*gfp; | ||
+ | |||
+ | dX=[dahl dlam dci dlaci dgfp]'; |
Revision as of 11:44, 24 October 2010
Modelling
Wiring Diagram
This is a outline of the regulatory topology of the concentration band-detect architecture.The different components and the way they interact with each other is represented similar to a electronic circuit.We use this as the starting point for building our model.First each of the components of our system was mathematically formulated. The components were then put together to form the whole system.
Model
Results
Code
odeIGEM03.m
%lambda1 function file
function dX =odeIGEM02(t,X,P)
k1=P(1); k2=P(2); k3=P(3); k4=P(4); k5=P(5); k6=P(6); k7=P(7); k8=P(8); k9=P(9); k10=P(10); k11=P(11);
km1=P(12); km2=P(13); km3=P(14); km4=P(15); km5=P(16);
ahl=X(1); lam=X(2); ci=X(3); laci=X(4); gfp=X(5);
n1=1; n2=1; n3=1; n4=1; n5=2;
dlam=((k1*(ahl/km1)^n1)/(1+(ahl/km1)^n1))-k2*lam; dahl=0;%k10*ahl-k11*ahl; dci=((k3*(lam/km3)^n2)/(1+(lam/km3)^n2))-k4*ci; dlaci=(k5/(1+(ci/km3)^n3))+((k6*(lam/km4)^n4)/(1+(lam/km4)^n4))-k7*laci; dgfp=(k8/(1+(laci/km5)^n5))-k9*gfp;
dX=[dahl dlam dci dlaci dgfp]';