#astro.ode #Li-Rinzel model for Ca-dynamics #This file as it is show correct q-nullcline but not a correct ca-nullcline - too high #parameters p ip3 p k3=.05 p tau_ip3=7.1429 p c0=2.0,c1=.185,v1=6.0,v2=.11,v3=.9 p d1=.13,d2=1.049,d3=.9434,d5=.08234 p a2=.2 ninf(ca)=ca/(ca+d5) jchan(ca,q)=c1*v1*(ip3/(ip3+d1))^3*ninf(ca)^3*q^3*(ca-caER(ca)) jpump(ca)=(v3*ca^2)/(k3^2+ca^2) jleak(ca)=c1*v2*(ca-caER(ca)) caER(ca)=(c0-ca)/c1 dca/dt=-jchan(ca,q)-jpump(ca)-jleak(ca) dq/dt=(a2*d2*(ip3+d1)/(ip3+d3))*(1-q)-a2*ca*q #initial conditions ca(0)=.1 q(0)=.01 #My part @ xp=ca, yp=q, xlo=0.0, ylo=0.4, xhi=1.5, yhi=1.2 @ total=1000, dt=.05, meth=backeul @ AUTOVAR=ca, AUTOXMIN=0, AUTOXMAX=1, AUTOYMIN=0, AUTOYMAX=1 @ DS=.001,DSMIN=.0001,DSMAX=.05, PARMIN=0, PARMAX=.9, NORMMAX=1000 done