#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