# Model of a genetic toggle switch - stochastic version # Lepzelter et al., J. Phys. Chem. B 111, 10239 (2007). param chA=0.0004, chB=0.0002, cfA=0.2, cfB=0.2 param bA=1, bB=1, cgA=150, cgB=100 param c=1 # Reaction propensities r1=chA*OonA*B*(B-1)/2 r2=cfA*OoffA r3=chB*OonB*A*(A-1)/2 r4=cfB*OoffB r5=cgA*OonA r6=cgB*OonB r7=c*A r8=c*B OonA(0)=0 OoffA(0)=1 OonB(0)=0 OoffB(0)=1 A(0)=0 B(0)=0 special z=gill(0,r{1-8}) time'=time+z(0) OonA'=OonA-z(1)+z(2) OoffA'=OoffA+z(1)-z(2) OonB'=OonB-z(3)+z(4) OoffB'=OoffB+z(3)-z(4) A'=A-2*z(3)+2*z(4)+bA*z(5)-z(7) B'=B-2*z(1)+2*z(2)+bB*z(6)-z(8) @ MAXSTOR=1000000, BOUND=1e6 @ TOTAL=200000, METH=discrete done