Homework+3+(genetic+toggle+switch)

The goal of this homework is to recreate Figure 2D from the 2000 Nature paper Construction of a genetic toggle switch in Escherichia coli



1st Attempt

ODE File from teacher


 * 1) Parameters

p k1=10 p k2=10 p n1=2 p n2=2


 * 1) Constants

u' = k1/(1 + v^n1) - u v' = k2/(1 + u^n2) - v

i u=0 i v=1

@ maxstore=100000, total=100, dt=0.01 @ xplot=u, yplot=v @ xlo=-1, xhi=10 ylo=-1, yhi=10
 * 1) Control Parameters

done

Here's some pointers from Ethan

axis: hi lo main param: V x: -20:40 y: -20:20

numerics: ds = 0.02 nmax = 400 npr = 200 dsmax = 2.5 / 3.5 par min = -220 par max = 50 norm max = 1000

run: steady state grab point 2

axes: 2 par x: -50:50 y: -40:60 numerics: par max: 100 run grab lowest point numerics: ds = -0.02 run again

When I tried to follow these instructions, nothing graphed and I didn't want to spend too much time on this.