## Set file for Fridlyand.ode 13 Number of equations and auxiliaries 74 Number of parameters # Numerical stuff 1 nout 40 nullcline mesh 10 CVode 1800000 total 50 DeltaT 0 T0 0 Transient 100000000 Bound 1e-12 DtMin 1 DtMax 1e-07 Tolerance 1e-07 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 1800000 Current time 1800000 Last Time 0 MyStart 1 INFLAG # Delays 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 # Old ICs -59.23737317695019 V 0.2940743635025411 CAI 23.69658362865903 ADPF 5995.577830221578 RE 3.805286572034495 IP3 6603.65536107568 NI 48.30133793695131 CAER 0.001557809673010048 N 240.9902445689719 ADPB 3734.513171802411 ATP # Ending ICs -74.65058290194277 V 0.03251968698081115 CAI 3.5170113830354 ADPF 145018.5249051946 RE 0.04924649729026075 IP3 6850.202214093814 NI 3.411997342163442 CAER 0.0001725986501678411 N 35.91222052782078 ADPB 3959.77076808946 ATP # Parameters 24000 gmKATP 9.5 KCah 8.000000000000001e-05 kATPCa 0 VCRAN 5e-05 kATP 0.7 gmCRAN 0.002 mf5 670 gmVCa 17 kdd 1.5e-08 f44 8000 Ko 1200 gmNa 0.4 KIPCa 6158 Cm 0.0001 Pleak 0.0001 ksg 500 KmATP 0.1 KKCa 140000 No -19 VCah 0.0003 kIP 1.7 hgl 0.03 mb5 4000 Ao 0.5 KCarp 0.01 fi 3.2 KIP3 1.72e-17 b33 6e-07 b66 132400 Ki 0.0002 b44 4950 P 20 Kop 4e-05 Kdip 0.0002 kADPf 2e-05 kADPb 0.0012 PIP3 600 PNaK 3000 gmKDr 0.03 fer 11.5 f6 0.172 f3 10 f2 2.5e-10 f11 8 Gluo 7 Kg 0.75 KNaCa 130 gmKCa 26.73 RTdF -75 Vt 20 c 0.28 Ver 20 b 2600 Cao 65 a -14 Vn 0.764 Vi 0.105 PCaER 271 gNaCa 2000 PmCap 96480 F 0.077 KRCa 50 ktt 0.025 Pglu 26 ktd 31 KRe 0.0001 b22 100 b1 0.0002 Pop 200 KCar 7 Sn 0.1 KCap 600000 Dzstart 1 Dzamp # Graphics 0.7071097251469876 rm -0.4999999999826661 rm 0.4999958366025516 rm 0.7071038372138505 rm 0.5000041633974484 rm -0.4999999999826661 rm 0 rm 0.7071038372138505 rm 0.7071097251469876 rm 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 12 2 2 1 0 -1000 1000 0 0 3DFlag 0 Timeflag 0 Colorflag 0 Type 1 color scale 0 minscale 30 xmax 0 xmin 0.8 ymax 0 ymin 12 zmax -12 zmin 15 0.06666666666666667 0.4 2.5 0 0.08333333333333333 45 Theta 45 Phi 0 xshft 0 yshft 0 zshft 0 xlo 0 ylo 0 0 30 xhi 0.8 yhi 30 0.8 # 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 # 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 gmKATP -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 gmKATP 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=-((IVCA+ICAPUMP+INACA+ICRAN+INA+INAK+IKDR+IKCA+IKATP)/CM) dCAI/dT=((FI*(((-IVCA+(2.0*INACA)-(2.0*ICAPUMP))/(2.0*F*VI))-JERP+(JOUT/VI)))-(KSG*CAI)) dADPF/dT=((KADPB*ADPB)-JOP-(KADPF*ADPF)+((INAK+ICAPUMP)/(F*VI))+(0.5*JERP)+((KATP+(KATPCA*CAI))*ATP)) dRE/dT=((KRE*JGLU)-JOP) dIP3/dT=((KIP*CAI*CAI/((CAI*CAI)+(KIPCA*KIPCA)))-(KDIP*IP3)) dNI/dT=((-(3.0*INACA)-(3.0*INAK)-INA-ICRAN)/(F*VI)) dCAER/dT=(FER*((JERP*VI)-JOUT)/VER) dN/dT=((NB-N)/TN) dADPB/dT=((KADPF*ADPF)-(KADPB*ADPB)) dATP/dT=(-(0.5*JERP)-((KATP+(KATPCA*CAI))*ATP)+JOP-((INAK+ICAPUMP)/(F*VI))) TSEC=T/1000.0 TMIN=T/60000.0 OKATP=OKATP where ... DZ = DZAMP*HEAV(T - DZSTART) F4 = (F44 * KO * KO) B6 = (B66 * KI * KI) B4 = (B44 * P) B3 = (B33 * NO * NO * NO) F1 = (F11 * NI * NI * NI) F55 = (MF5 * EXP((V / (2.0 * RTDF)))) F5 = (F55 * ATP) B2 = (B22 * (AO - ATP)) B5 = (MB5 * EXP( - (V / (2.0 * RTDF)))) D = ((F2 * F3 * F4 * F5 * F6) + (B1 * F3 * F4 * F5 * F6) + (B1 * B2 * F4 * F5 * F6) + (B1 * B2 * B3 * F5 * F6) + (B1 * B2 * B3 * B4 * F6) + (B1 * B2 * B3 * B4 * B5)) INAK = (PNAK * ((F1 * F2 * F3 * F4 * F5 * F6) - (B1 * B2 * B3 * B4 * B5 * B6)) / D) JERP = (PCAER * CAI * CAI / ((CAI * CAI) + (KCARP * KCARP))) VNACA = (RTDF * ((3.0 * LOG((NO / NI))) - LOG((CAO / CAI)))) INACA = (GNACA * (V - VNACA) * CAI * CAI * CAI * CAI * CAI / ((KNACA * KNACA * KNACA * KNACA * KNACA) + (CAI * CAI * CAI * CAI * CAI))) VNA = (RTDF * LOG((NO / NI))) INA = (GMNA * (V - VNA) / (1.0 + EXP((0.125 * (104.0 + V))))) TN = (C / (EXP(((V - VT) / A)) + EXP(((VT - V) / B)))) PVCA = (1.0 / (1.0 + EXP(((VCAH - V) / KCAH)))) VCA = (0.5 * LOG((CAO / CAI)) * RTDF) IVCA = (GMVCA * PVCA * (V - VCA)) MGADP = (0.55 * ADPF) OKATP = (((0.08 * (1.0 + (2.0 * MGADP / KDD))) + (0.89 * MGADP * MGADP / KDD / KDD)) / ((1.0 + (MGADP / KDD)) * (1.0 + (MGADP / KDD)) * (1.0 + (0.45 * MGADP / KTD) + (ATP / KTT)))) VK = (RTDF * LOG((KO / KI))) IKATP = ((1 - DZ)*GMKATP*OKATP + DZ*GMKATP) * (V - VK) ICRAN = - ((V - VNA) * (V - VCRAN) * GMCRAN / (1.0 + EXP((0.3333333333333333 * (CAER - KCAR))))) JOP = (POP * RE * MGADP * MGADP / ((MGADP * MGADP) + (KOP * KOP))) JOUT = ((PLEAK + (PIP3 * CAI * IP3 * IP3 * IP3 / (CAI + KRCA) / ((IP3 * IP3 * IP3) + (KIP3 * KIP3 * KIP3)))) * (CAER - CAI)) IKDR = (GMKDR * N * (V - VK)) ICAPUMP = (PMCAP * CAI * CAI / ((CAI * CAI) + (KCAP * KCAP))) NB = (1.0 / (1.0 + EXP(((VN - V) / SN)))) JGLU = (PGLU * ATP * (GLUO ^ HGL) / (KMATP + ATP) / ((KG ^ HGL) + (GLUO ^ HGL))) IKCA = (GMKCA * (V - VK) * CAI * CAI * CAI * CAI / ((KKCA * KKCA * KKCA * KKCA) + (CAI * CAI * CAI * CAI)))