# Bullfrog Sympathetic Ganglion Neuron, open cell, Ca2+ oscillations # Model as in book and my notes (see comments at bottom of file) # p jin=1.5 # 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 p vpmca=5 p Kpmca=.6 i Cai=.1 i w=.5 i CaT=2.2 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) jpmca(Cai)=vpmca*Cai^2/(Cai^2+Kpmca^2) Cai'=fi*(jRyR(Cai)+jlk(Cai)-jserca(Cai)+jin-jpmca(Cai)) w'=(winf(Cai)-w)/tau(Cai) CaT'=fi*(jin-jpmca(Cai)) done