# Bullfrog Sympathetic Ganglion Neuron, closed cell, Ca2+ oscillations
#  Model as in book and my notes
#  
p CaT=1.8
#  Try also CaT=2.2, 3
p fi=.01
p vRyR=5
p vlk=.2
p vserca=100
p Kserca=.2
p sigma=.02
p Ka=.4
p Kb=.6
p Kc=.1
p kcm=.1
p jmin=0
p jmout=0

i Cai=.1
i w=.5

winf(Cai)=(1+(Ka/Cai)^4+(Cai/Kb)^3)/(1+1/Kc+(Ka/Cai)^4+(Cai/Kb)^3)
Po(Cai)=w*(1+(Cai/Kb)^3)/(1+(Ka/Cai)^4+(Cai/Kb)^3)
tau(Cai)=winf(Cai)/kcm
CaER(Cai)=(CaT-Cai)/sigma
jRyR(Cai)=vRyR*Po(Cai)*(CaER(Cai)-Cai)
jlk(Cai)=vlk*(CaER(Cai)-Cai)
jserca(Cai)=vserca*Cai^2/(Cai^2+Kserca^2)

Cai'=fi*(jRyR(Cai)+jlk(Cai)-jserca(Cai))
w'=(winf(Cai)-w)/tau(Cai)

@ xlo=0,xhi=4, ylo=0, yhi=1.5, xp=Cai, yp=w

done
#
#  For winpp set axes to display Cai vs w  with ranges  Cai: 0 to 4, w: 0 to 1.5
#  For locbif, set axes to display CaT vs Cai with ranges 0 to 4