Team:IIT Delhi/Modeling
From 2014.igem.org
Revision as of 23:26, 15 October 2014 by AbhishekBharti (Talk | contribs)
i have the equations with the data; let x=mRNA, y=protein, z=NOx; then
dx/dt = 4.56 - 0.00345*x x(0)=0
dy/dt = 0.1*x - 1.148*10^-3*y y(0)=0
now i facing problem with NOx rate equation; i want to apply ,"input_rate - output_rate + generation_rate - consumption_rate = accumulation_rate".
and the equation is; dz/dt= - 770*y*z/(0.00003+z), z(0)=0; assuming no inlet of 'z'; doing this matlab gives z=0, for all t. thats the problem.
also dz/dt=1000 - 770*y*z/0.00003+z,z(0)=0', assuming 1000 NOx molecules enter per second in the cell, this doesn't give any explicit solution.
now if remove 'y' and put any numeric value of protein concentration, solving for 'z' gives a result but it shows all molecules of NOx are converted in less than a second. another error i guess
graphs for mRNA and protein concentration with time are attached.
graphs for mRNA and Protein are ok, but NOx levels are giving problem.