## Set file for single_cells.ode on Wed Dec 13 10:25:52 2017
19 Number of equations and auxiliaries
72 Number of parameters
# Numerical stuff
1 nout
40 nullcline mesh
8 Qual RK
4000000 total
0.5 DeltaT
0 T0
0 Transient
500000 Bound
1e-12 DtMin
1 DtMax
1e-07 Tolerance
0.001 Abs. Tolerance
0 Max Delay
100 Eigenvector iterates
0.001 Eigenvector tolerance
0.001 Newton tolerance
0 Poincare plane
1e-05 Boundary value tolerance
1e-05 Boundary value epsilon
20 Boundary value iterates
0 Poincare None
1 Poincare variable
1 Poincare sign
0 Stop on Section
0 Delay flag
4000000 Current time
4000000 Last Time
1 MyStart
1 INFLAG
# Delays
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
# Bndry conds
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Old ICs
-58.81137764996264 V
0.432356959113168 N
0.8538036058454112 HP
0.002082416993486555 SG
0.08000000000046259 IMAXA
-39.02998140466573 VA
149.1319593988436 KI
16.41369332787475 KE
1.285281166395606 NAI
116.0901080839575 NAE
71.95968297063686 KIA
66.24314341750998 NAIA
12.52512787448607 CLI
111.3857579081404 CLE
0.004525191854168424 GLUT
0.0009387716976116532 CAI
0.05904737645069838 CAER
2.801142735550066 CAE
0.9981259767393746 HER
# Ending ICs
-62.52638004813665 V
0.3687463987794178 N
0.9166237311852015 HP
0.3333333333333334 SG
0.08000000000046259 IMAXA
-26.3566736314976 VA
138.9603521990604 KI
10.84005294617766 KE
7.222072086379609 NAI
129.0789158035363 NAE
29.67248554159209 KIA
108.5285057274632 NAIA
11.03269938780366 CLI
112.0102445460947 CLE
0.05796814244940238 GLUT
0.00176875116910499 CAI
0.1920086511206416 CAER
0.1877525063186725 CAE
0.9964749701503084 HER
# Parameters
8310 R
310 Temp
96485 F
1e-05 pna
2e-09 pnal
-34 thm
5 sigm
3e-08 pnap
10000 taubar
-40 thmp
6 sigmp
-48 thhp
-6 sighp
-49 vt
6 sig
0.05 phihp
6.999999999999999e-05 pk
4e-07 pkl
-55 thn
14 sgn
0.05 taun0
0.27 taun1
-40 thnt
-12 sn
0.8 phin
3e-06 pnmda
-10 thetat
2 trise
1 tdecay
0.5 alphag
3 pca
0.01 thg
0.001 sigmag
2e-07 pcl
5 imax
4.8e-06 pka
1.5e-08 pnaa
3e-05 eps
0.08 imin
0.5 hr
3e-05 gnagl
2160 voln
922 s
1600 Sa
2000 vola
0.05 alpha0
1e-05 gnab
2e-05 gkb
6e-06 gclb
140 na0
4 k0
110 cl0
1 glui0
0.001 glut0
60 minute
1000 second
0.3 ip3
0.5 pcytosol
75 dcytosol
0.1 per
1000 der
0.01 fi
3000 vip3
0.01 vleak
0.25 dip3
0.5 dinh
1 dact
4 tau
110 vserca
0.4 kserca
6e-06 gcab
3 cae0
# Graphics
0.7071097251469876 rm
-0.4999999999826661 rm
0.4999958366025516 rm
0.7071038372138505 rm
0.5000041633974482 rm
-0.4999999999826661 rm
0 rm
0.7071038372138505 rm
0.7071097251469876 rm
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
0
1
1
1
0
-1000
1000
0
0 3DFlag
1 Timeflag
0 Colorflag
0 Type
1 color scale
0 minscale
4000000 xmax
0 xmin
0 ymax
-90 ymin
12 zmax
-12 zmin
2000000
5e-07
-45
0.02222222222222222
0
0.08333333333333333
45 Theta
45 Phi
0 xshft
0 yshft
0 zshft
0 xlo
-90 ylo
0
-80
4000000 xhi
0 yhi
50000
0
# Transpose variables etc
V
2 n columns
1 n rows
1 row skip
1 col skip
1 row 0
# Coupling stuff for H funs
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
# Array plot stuff
1 NCols
0 Row 1
50 NRows
8 RowSkip
0 Zmin
1 Zmax
# Torus information
0 Torus flag 1=ON
6.283185307179586 Torus period
# Range information
R
-1 eq-range stab col
0 shoot flag 1=on
10 eq-range steps
0 eq_range low
1 eq_range high
V
V
20 Range steps
0 Cycle color 1=on
1 Reset data 1=on
1 Use old I.C.s 1=yes
0 Par1 low
0 Par2 low
1 Par1 high
1 Par2 high
R
0 BVP side
0 color cycle flag 1=on
10 BVP range steps
0 BVP range low
1 BVP range high
RHS etc ...
dV/dT=-(INA(V,N,NAE,NAI)+INAP(V,HP,NAE,NAI)+IK(V,N,KE,KI)+ICL(V,CLE,CLI)+IPUMP+INMDA)
dN/dT=PHIN*(NINF(V)-N)/TAUN(V)
dHP/dT=PHIHP*(HINFP(V)-HP)/TAUHP(V)
dSG/dT=-SG/TDECAY+ALPHAG*SINFG(GLUT)*(1-SG)
dIMAXA/dT=EPS*(IMIN-IMAXA)
dVA/dT=-(IKA(VA,KE,KIA)+INAA(VA,NAE,NAIA)+IPUMPA(KE,NAIA)+INAGL)
dKI/dT=-CFN*(IK(V,N,KE,KI)+IKPUMP+IKNMDA)
dKE/dT=CFNE*(IK(V,N,KE,KI)+IKPUMP+IKNMDA)+CFAE*(IKA(VA,KE,KIA)-2*IPUMPA(KE,NAIA)-INAGL)+GKB*(K0-KE)
dNAI/dT=-CFN*(INA(V,N,NAE,NAI)+INAP(V,HP,NAE,NAI)+INAPUMP+INANMDA)
dNAE/dT=CFNE*(INA(V,N,NAE,NAI)+INAP(V,HP,NAE,NAI)+INAPUMP+INANMDA)+CFAE*(INAA(VA,NAE,NAIA)+3*IPUMPA(KE,NAIA)+(3)*INAGL)+GNAB*(NA0-NAE)
dKIA/dT=-CFA*(IKA(VA,KE,KIA)-2*IPUMPA(KE,NAIA)-INAGL)
dNAIA/dT=-CFA*(INAA(VA,NAE,NAIA)+3*IPUMPA(KE,NAIA)+(3)*INAGL)
dCLI/dT=CFN*ICL(V,CLE,CLI)
dCLE/dT=-CFNE*ICL(V,CLE,CLI)+GCLB*(CL0-CLE)
dGLUT/dT=CFAE*INAGL
dCAI/dT=-FI*CFN*(ICANMDA)-FI*(PER*DER*(JSERCA-JEROUT))/(PCYTOSOL*MINUTE*SECOND)
dCAER/dT=FI*(PER*DER*(JSERCA-JEROUT))/(PER*MINUTE*SECOND)
dCAE/dT=CFNE*ICANMDA+GCAB*(CAE0-CAE)
dHER/dT=(DINH-(CAI+DINH)*HER)/TAU
where ...
FRT = R*TEMP/F
INANMDA = PNMDA*SG*BINF(V)*F*PHI(V)*(NAE*EXP(-PHI(V))-NAI)/(EXP(-PHI(V))-1)
IKNMDA = PNMDA*SG*BINF(V)*F*PHI(V)*(KE*EXP(-PHI(V))-KI)/(EXP(-PHI(V))-1)
ICANMDA = PCA*2*PNMDA*SG*BINF(V)*F*PHI(2*V)*(CAE*EXP(-PHI(2*V))-CAI)/(EXP(-PHI(2*V))-1)
INMDA = INANMDA+IKNMDA+ICANMDA
IPUMP = IMAX/(((1+2/KE)^2)*(1+7.7/NAI)^3)
INAPUMP = 3*IPUMP
IKPUMP = -2*IPUMP
ENAGL = (FRT/2)*LN(((NAE/NAIA)^3)*(KIA/KE)*(GLUT/GLUI)*HR)
INAGL = GNAGL*(VA-ENAGL)
VOLE = ALPHA0*VOLN
CFN = 10.0*S/(F*VOLN)
CFNE = 10.0*S/(F*VOLE)
CFA = 10*SA/(F*VOLA)
CFAE = 10*SA/(F*VOLE)
GLTOT = CFAE*GLUI0+CFA*GLUT0
GLUI = (1/CFAE)*(GLTOT-CFA*GLUT)
JEROUT = (VIP3*((IP3/(IP3+DIP3))^3)*((CAI/(CAI+DACT))^3)*(HER^3)+VLEAK)*(CAER-CAI)
JSERCA = VSERCA*CAI^2/(KSERCA^2+CAI^2)
User-defined functions:
PHI(V) = V/FRT
MINF(V) = 1./(1.+EXP(-(V-THM)/SIGM))
INA(V,N,NAE,NAI) = PNA*(MINF(V)^3*(1-N)+PNAL)*F*PHI(V)*(NAE*EXP(-PHI(V))-NAI)/(EXP(-PHI(V))-1)
MINFP(V) = 1./(1.+EXP(-(V-THMP)/SIGMP))
HINFP(V) = 1./(1.+EXP(-(V-THHP)/SIGHP))
TAUHP(V) = TAUBAR/COSH((V-VT)/(2*SIG))
INAP(V,HP,NAE,NAI) = PNAP*MINFP(V)*HP*F*PHI(V)*(NAE*EXP(-PHI(V))-NAI)/(EXP(-PHI(V))-1)
NINF(V) = 1./(1.+EXP(-(V-THN)/SGN))
TAUN(V) = TAUN0+TAUN1/(1+EXP(-(V-THNT)/SN))
IK(V,N,KE,KI) = PK*(N^4+PKL)*F*PHI(V)*(KE*EXP(-PHI(V))-KI)/(EXP(-PHI(V))-1)
BINF(V) = 1/(1+EXP(-(V-THETAT)/16.13))
SINFG(X) = 1./(1.+EXP(-(X-THG)/SIGMAG))
ICL(V,CLE,CLI) = -PCL*F*PHI(-V)*(CLE*EXP(-PHI(-V))-CLI)/(EXP(-PHI(-V))-1)
IKA(V,KE,KI) = PKA*F*PHI(V)*(KE*EXP(-PHI(V))-KI)/(EXP(-PHI(V))-1)
INAA(V,NAE,NAI) = PNAA*F*PHI(V)*(NAE*EXP(-PHI(V))-NAI)/(EXP(-PHI(V))-1)
IPUMPA(KE,NAI) = IMAXA/(((1+2/KE)^2)*(1+10/NAI)^3)
