#tyson.ode # Equations Y'=k1-(k2p+k2pp*P)*Y P'=(1-P)*(k3p+k3pp*A)/(J3+1-P)-P*k4*m*Y/(J4+P) A'=k5p+k5pp* ((m*Y/J5)^n)/(1+(y*m/J5)^n)-k6*A # Parameters with units of 1/time: p k1=0.04 p k2p=0.04,k2pp=1 p k3p=1,k3pp=10 p k4=35 p k5p=0.005,k5pp=0.2,k6=0.1 # mass of cell (try m=0.6, m=0.3, m=1) p m=1 # Dimensionless parameters: p J3=0.04,J4=0.04,J5=0.3,n=4 # Numerics @ TOTAL=2000,DT=.1,xlo=0,xhi=2000,ylo=0,yhi=6 @ NPLOT=1,XP1=t,YP1=Y @ MAXSTOR=10000000 @ BOUNDS=100000,METH=stiff @ dsmin=1e-5,dsmax=.1,parmin=-.5,parmax=.5,autoxmin=-.5,autoxmax=.5 @ autoymax=.4,autoymin=-.5 # IC Y(0)=1 P(0)=0.5 A(0)=0.1 done