$\mathrm{dVdt}$ $\begin{cases}\mathrm{stim\_amplitude} & \text{if $(\mathrm{time\_}\ge \mathrm{stim\_start})\land (\mathrm{time\_}\le \mathrm{stim\_end})\land (\mathrm{time\_}-\mathrm{stim\_start}-\lfloor \frac{\mathrm{time\_}-\mathrm{stim\_start}}{\mathrm{stim\_period}}\rfloor \mathrm{stim\_period}\le \mathrm{stim\_duration})$}\\ 0 & \text{otherwise}\end{cases}$ $-(\mathrm{i\_Na}+\mathrm{i\_Ca\_L}+\mathrm{i\_Ca\_T}+\mathrm{i\_Kr}+\mathrm{i\_Ks}+\mathrm{i\_K\_Na}+\mathrm{i\_K\_ATP}+\mathrm{i\_to}+\mathrm{i\_K1}+\mathrm{i\_Kp}+\mathrm{i\_NaCa}+\mathrm{i\_p\_Ca}+\mathrm{i\_Na\_b}+\mathrm{i\_Ca\_b}+\mathrm{i\_NaK}+\mathrm{i\_ns\_Ca}+\mathrm{I\_st})$ $\mathrm{alpha\_m}(1-m)-\mathrm{beta\_m}m$ $\begin{cases}\frac{0.32\mathrm{E0\_m}}{1-e^{-0.1\mathrm{E0\_m}}} & \text{if $\left|\mathrm{E0\_m}\right|\ge \mathrm{delta\_m}$}\\ 3.2 & \text{otherwise}\end{cases}$ $0.08e^{\frac{-V}{11}}$ $V+47.13$ $\mathrm{alpha\_h}(1-h)-\mathrm{beta\_h}h$ $\begin{cases}0.135e^{\frac{80+V}{-6.8}} & \text{if $V< -40$}\\ 0 & \text{otherwise}\end{cases}$ $\begin{cases}3.56e^{0.079V}+310000e^{0.35V} & \text{if $V< -40$}\\ \frac{1}{0.13(1+e^{\frac{V+10.66}{-11.1}})} & \text{otherwise}\end{cases}$ $\mathrm{alpha\_j}(1-j)-\mathrm{beta\_j}j$ $\begin{cases}\frac{-(127140e^{0.2444V}+3.474E-5e^{-0.04391V})(V+37.78)}{1+e^{0.311(V+79.23)}} & \text{if $V< -40$}\\ 0 & \text{otherwise}\end{cases}$ $\begin{cases}\frac{0.1212e^{-0.01052V}}{1+e^{-0.1378(V+40.14)}} & \text{if $V< -40$}\\ \frac{0.3e^{-2.535E-7V}}{1+e^{-0.1(V+32)}} & \text{otherwise}\end{cases}$ $\mathrm{g\_Na}m^{3}hj(V-\mathrm{E\_Na})$ $\frac{RT}{F}\ln \left(\frac{\mathrm{Nao}}{\mathrm{Nai}}\right)$ $\mathrm{alpha\_d}(1-d)-\mathrm{beta\_d}d$ $\frac{\mathrm{d\_infinity}}{\mathrm{tau\_d}}$ $\frac{1-\mathrm{d\_infinity}}{\mathrm{tau\_d}}$ $\frac{1}{1+e^{\frac{-\mathrm{E0\_d}}{6.24}}}$ $\begin{cases}\frac{1}{0.035\times 6.24\times 2} & \text{if $\left|\mathrm{E0\_d}\right|< 1E-5$}\\ \frac{1\mathrm{d\_infinity}(1-e^{\frac{-\mathrm{E0\_d}}{6.24}})}{0.035\mathrm{E0\_d}} & \text{otherwise}\end{cases}$ $V+10$ $\mathrm{alpha\_f}(1-f)-\mathrm{beta\_f}f$ $\frac{\mathrm{f\_infinity}}{\mathrm{tau\_f}}$ $\frac{1-\mathrm{f\_infinity}}{\mathrm{tau\_f}}$ $\frac{1}{1+e^{\frac{V+35.06}{8.6}}}+\frac{0.6}{1+e^{\frac{50-V}{20}}}$ $\frac{1}{0.0197e^{-(0.0337(V+10))^2}+0.02}$ $\frac{1}{1+\frac{\mathrm{Cai}}{\mathrm{Km\_Ca}}}$ $\mathrm{i\_CaCa}+\mathrm{i\_CaK}+\mathrm{i\_CaNa}$ $df\mathrm{f\_Ca}\mathrm{I\_CaCa}$ $df\mathrm{f\_Ca}\mathrm{I\_CaK}$ $df\mathrm{f\_Ca}\mathrm{I\_CaNa}$ $\frac{\frac{\mathrm{P\_Ca}2^{2}VF^2}{RT}(\mathrm{gamma\_Cai}\mathrm{Cai}e^{\frac{2VF}{RT}}-\mathrm{gamma\_Cao}\mathrm{Cao})}{e^{\frac{2VF}{RT}}-1}$ $\frac{\frac{\mathrm{P\_K}1^{2}VF^2}{RT}(\mathrm{gamma\_Ki}\mathrm{Ki}e^{\frac{1VF}{RT}}-\mathrm{gamma\_Ko}\mathrm{Ko})}{e^{\frac{1VF}{RT}}-1}$ $\frac{\frac{\mathrm{P\_Na}1^{2}VF^2}{RT}(\mathrm{gamma\_Nai}\mathrm{Nai}e^{\frac{1VF}{RT}}-\mathrm{gamma\_Nao}\mathrm{Nao})}{e^{\frac{1VF}{RT}}-1}$ $\frac{\mathrm{b\_inf}-b}{\mathrm{tau\_b}}$ $\frac{1}{1+e^{\frac{-(V+14)}{10.8}}}$ $3.7+\frac{6.1}{1+e^{\frac{V+25}{4.5}}}$ $\frac{\mathrm{g\_inf}-g}{\mathrm{tau\_g}}$ $\frac{1}{1+e^{\frac{V+60}{5.6}}}$ $\begin{cases}0.875V+12 & \text{if $V\le 0$}\\ 12 & \text{otherwise}\end{cases}$ $\mathrm{g\_CaT}bbg(V-\mathrm{E\_Ca})$ $\frac{\mathrm{xr\_infinity}-\mathrm{xr}}{\mathrm{tau\_xr}}$ $\frac{1}{1+e^{\frac{-(V+21.5)}{7.5}}}$ $\frac{1}{\frac{0.00138(V+14.2)}{1-e^{-0.123(V+14.2)}}+\frac{0.00061(V+38.9)}{e^{0.145(V+38.9)}-1}}$ $\mathrm{g\_Kr}\mathrm{xr}\mathrm{Rect}(V-\mathrm{E\_K})$ $\mathrm{G\_Kr}\sqrt{\frac{\mathrm{Ko}}{5.4}}$ $\frac{1}{1+e^{\frac{V+9}{22.4}}}$ $\frac{\mathrm{xs1\_infinity}-\mathrm{xs1}}{\mathrm{tau\_xs1}}$ $\frac{1}{1+e^{\frac{-(V-1.5)}{16.7}}}$ $\frac{1}{\frac{7.19E-5(V+30)}{1-e^{-0.148(V+30)}}+\frac{0.000131(V+30)}{e^{0.0687(V+30)}-1}}$ $\frac{\mathrm{xs2\_infinity}-\mathrm{xs2}}{\mathrm{tau\_xs2}}$ $\frac{1}{1+e^{\frac{-(V-1.5)}{16.7}}}$ $\frac{4}{\frac{7.19E-5(V+30)}{1-e^{-0.148(V+30)}}+\frac{0.000131(V+30)}{e^{0.0687(V+30)}-1}}$ $\mathrm{g\_Ks}\mathrm{xs1}\mathrm{xs2}(V-\mathrm{E\_Ks})$ $\mathrm{G\_Ks}(1+\frac{0.6}{1+\left(\frac{3.8E-5}{\mathrm{Cai}}\right)^{1.4}})$ $\frac{RT}{F}\ln \left(\frac{\mathrm{Ko}+\mathrm{PNaK}\mathrm{Nao}}{\mathrm{Ki}+\mathrm{PNaK}\mathrm{Nai}}\right)$ $\frac{\mathrm{alpha\_K1}}{\mathrm{alpha\_K1}+\mathrm{beta\_K1}}$ $\frac{1.02}{1+e^{0.2385(V-\mathrm{E\_K}-59.215)}}$ $\frac{1(0.49124e^{0.08032(V-\mathrm{E\_K}+5.476)}+e^{0.06175(V-\mathrm{E\_K}-594.31)})}{1+e^{-0.5143(V-\mathrm{E\_K}+4.753)}}$ $\mathrm{g\_K1}\mathrm{K1\_infinity}(V-\mathrm{E\_K})$ $\frac{RT}{F}\ln \left(\frac{\mathrm{Ko}}{\mathrm{Ki}}\right)$ $\mathrm{G\_K1}\sqrt{\frac{\mathrm{Ko}}{5.4}}$ $\mathrm{g\_Kp}\mathrm{Kp}(V-\mathrm{E\_K})$ $\frac{1}{1+e^{\frac{7.488-V}{5.98}}}$ $\mathrm{g\_K\_Na}\mathrm{pona}\mathrm{pov}(V-\mathrm{E\_K})$ $\frac{0.85}{1+\left(\frac{\mathrm{kdKNa}}{\mathrm{Nai}}\right)^{\mathrm{nKNa}}}$ $0.8-\frac{0.65}{1+e^{\frac{V+125}{15}}}$ $\mathrm{GKbaraATP}(V-\mathrm{E\_K})$ $\frac{\mathrm{i\_K\_ATP\_on}\times 0.000193}{\mathrm{nicholsarea}}$ $\frac{1}{1+\left(\frac{\mathrm{ATPi}}{\mathrm{kATP}}\right)^{\mathrm{hATP}}}$ $\mathrm{g\_K\_ATP}\mathrm{pATP}\left(\frac{\mathrm{Ko}}{4}\right)^{\mathrm{nATP}}$ $\frac{\mathrm{zdv\_ss}-\mathrm{zdv}}{\mathrm{tau\_zdv}}$ $\frac{10e^{\frac{V-40}{25}}}{1+e^{\frac{V-40}{25}}}$ $\frac{10e^{\frac{-(V+90)}{25}}}{1+e^{\frac{-(V+90)}{25}}}$ $\frac{1}{\mathrm{alpha\_zdv}+\mathrm{beta\_zdv}}$ $\frac{\mathrm{alpha\_zdv}}{\mathrm{alpha\_zdv}+\mathrm{beta\_zdv}}$ $\frac{\mathrm{ydv\_ss}-\mathrm{ydv}}{\mathrm{tau\_ydv}}$ $\frac{0.015}{1+e^{\frac{V+60}{5}}}$ $\frac{0.1e^{\frac{V+25}{5}}}{1+e^{\frac{V+25}{5}}}$ $\frac{1}{\mathrm{alpha\_ydv}+\mathrm{beta\_ydv}}$ $\frac{\mathrm{alpha\_ydv}}{\mathrm{alpha\_ydv}+\mathrm{beta\_ydv}}$ $\mathrm{g\_to}\mathrm{zdv}^3\mathrm{ydv}\mathrm{rvdv}(V-\mathrm{E\_K})$ $0\times 0.5$ $e^{\frac{V}{100}}$ $\frac{\mathrm{I\_pCa}\mathrm{Cai}}{\mathrm{K\_mpCa}+\mathrm{Cai}}$ $\mathrm{g\_Nab}(V-\mathrm{E\_Na})$ $\mathrm{g\_Cab}(V-\mathrm{E\_Ca})$ $\frac{RT}{2F}\ln \left(\frac{\mathrm{Cao}}{\mathrm{Cai}}\right)$ $\frac{\frac{\mathrm{I\_NaK}\mathrm{f\_NaK}\times 1}{1+\left(\frac{\mathrm{K\_mNai}}{\mathrm{Nai}}\right)^2}\mathrm{Ko}}{\mathrm{Ko}+\mathrm{K\_mKo}}$ $\frac{1}{1+0.1245e^{\frac{-0.1VF}{RT}}+0.0365\mathrm{sigma}e^{\frac{-VF}{RT}}}$ $\frac{1}{7}(e^{\frac{\mathrm{Nao}}{67.3}}-1)$ $\mathrm{i\_ns\_Na}+\mathrm{i\_ns\_K}$ $\frac{\mathrm{I\_ns\_Na}\times 1}{1+\left(\frac{\mathrm{K\_m\_ns\_Ca}}{\mathrm{Cai}}\right)^3}$ $\frac{\mathrm{I\_ns\_K}\times 1}{1+\left(\frac{\mathrm{K\_m\_ns\_Ca}}{\mathrm{Cai}}\right)^3}$ $1.75E-7$ $\frac{\frac{\mathrm{P\_ns\_Ca}1^{2}VF^2}{RT}(\mathrm{gamma\_Nai}\mathrm{Nai}e^{\frac{1VF}{RT}}-\mathrm{gamma\_Nao}\mathrm{Nao})}{e^{\frac{1VF}{RT}}-1}$ $\frac{\frac{\mathrm{P\_ns\_Ca}1^{2}VF^2}{RT}(\mathrm{gamma\_Ki}\mathrm{Ki}e^{\frac{1VF}{RT}}-\mathrm{gamma\_Ko}\mathrm{Ko})}{e^{\frac{1VF}{RT}}-1}$ $\frac{\mathrm{c1}e^{\frac{(\mathrm{gamma}-1)VF}{RT}}(e^{\frac{VF}{RT}}\mathrm{Nai}^3\mathrm{Cao}-\mathrm{Nao}^3\mathrm{Cai})}{1+\mathrm{c2}e^{\frac{(\mathrm{gamma}-1)VF}{RT}}(e^{\frac{VF}{RT}}\mathrm{Nai}^3\mathrm{Cao}+\mathrm{Nao}^3\mathrm{Cai})}$ $\mathrm{G\_rel}(\mathrm{Ca\_JSR}-\mathrm{Cai})$ $\frac{\mathrm{I\_up}\mathrm{Cai}}{\mathrm{Cai}+\mathrm{K\_mup}}$ $\mathrm{K\_leak}\mathrm{Ca\_NSR}$ $\frac{\mathrm{Ca\_NSR}-\mathrm{Ca\_JSR}}{\mathrm{tau\_tr}}$ $\frac{\mathrm{RyRopen}\mathrm{RyRclose}\times 150}{1+e^{\frac{\mathrm{i\_CaCa}+\mathrm{i\_Ca\_T}-2\mathrm{i\_NaCa}+\mathrm{i\_p\_Ca}+\mathrm{i\_Ca\_b}+5}{0.9}}}$ $\begin{cases}\frac{\mathrm{G\_rel\_max}(\mathrm{Cainfluxtrack}-\mathrm{delta\_Ca\_ith})}{\mathrm{K\_mrel}+\mathrm{Cainfluxtrack}-\mathrm{delta\_Ca\_ith}}(1-\mathrm{APtrack2})\mathrm{APtrack2} & \text{if $\mathrm{Cainfluxtrack}> \mathrm{delta\_Ca\_ith}$}\\ \mathrm{G\_rel\_overload}(1-\mathrm{OVRLDtrack2})\mathrm{OVRLDtrack2} & \text{if $(\mathrm{Cainfluxtrack}\le \mathrm{delta\_Ca\_ith})\land (\mathrm{OVRLDtrack2}> 0)$}\\ 0 & \text{otherwise}\end{cases}$ $\frac{\mathrm{I\_up}}{\mathrm{Ca\_NSR\_max}}$ $\frac{1}{1+\frac{\mathrm{CMDN\_max}\mathrm{K\_mCMDN}}{(\mathrm{K\_mCMDN}+\mathrm{Cai})^2}+\frac{\mathrm{Tn\_max}\mathrm{K\_mTn}}{(\mathrm{K\_mTn}+\mathrm{Cai})^2}}(\frac{-1\mathrm{A\_cap}(\mathrm{i\_CaCa}+\mathrm{i\_Ca\_T}-2\mathrm{i\_NaCa}+\mathrm{i\_p\_Ca}+\mathrm{i\_Ca\_b})}{2\mathrm{V\_myo}F}+\frac{\mathrm{i\_rel}\mathrm{V\_JSR}}{\mathrm{V\_myo}}+\frac{(\mathrm{i\_leak}-\mathrm{i\_up})\mathrm{V\_NSR}}{\mathrm{V\_myo}})$ $\frac{1}{1+\frac{\mathrm{CSQN\_max}\mathrm{K\_mCSQN}}{(\mathrm{K\_mCSQN}+\mathrm{Ca\_JSR})^2}}(\mathrm{i\_tr}-\mathrm{i\_rel})$ $\frac{-\mathrm{i\_tr}\mathrm{V\_JSR}}{\mathrm{V\_NSR}}-\mathrm{i\_leak}+\mathrm{i\_up}$ $\begin{cases}100(1-\mathrm{APtrack})-0.5\mathrm{APtrack} & \text{if $\mathrm{dVdt}> 150$}\\ -0.5\mathrm{APtrack} & \text{otherwise}\end{cases}$ $\begin{cases}100(1-\mathrm{APtrack2})-0.5\mathrm{APtrack2} & \text{if $(\mathrm{APtrack}< 0.2)\land (\mathrm{APtrack}> 0.18)$}\\ -0.5\mathrm{APtrack2} & \text{otherwise}\end{cases}$ $\begin{cases}100(1-\mathrm{APtrack3})-0.5\mathrm{APtrack3} & \text{if $(\mathrm{APtrack}< 0.2)\land (\mathrm{APtrack}> 0.18)$}\\ -0.01\mathrm{APtrack3} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{-1\mathrm{A\_cap}(\mathrm{i\_CaCa}+\mathrm{i\_Ca\_T}-2\mathrm{i\_NaCa}+\mathrm{i\_p\_Ca}+\mathrm{i\_Ca\_b})}{2\mathrm{V\_myo}F} & \text{if $\mathrm{APtrack}> 0.2$}\\ 0 & \text{if $(\mathrm{APtrack2}> 0.01)\land (\mathrm{APtrack}\le 0.2)$}\\ -0.5\mathrm{Cainfluxtrack} & \text{otherwise}\end{cases}$ $\begin{cases}50(1-\mathrm{OVRLDtrack}) & \text{if $(\frac{1}{1+\frac{\mathrm{K\_mCSQN}}{\mathrm{Ca\_JSR}}}> \mathrm{CSQNthresh})\land (\mathrm{OVRLDtrack3}< 0.37)\land (\mathrm{APtrack3}< 0.37)$}\\ -0.5\mathrm{OVRLDtrack} & \text{otherwise}\end{cases}$ $\begin{cases}50(1-\mathrm{OVRLDtrack2}) & \text{if $(\mathrm{OVRLDtrack}> \mathrm{Logicthresh})\land (\mathrm{OVRLDtrack2}< \mathrm{Logicthresh})$}\\ -0.5\mathrm{OVRLDtrack2} & \text{otherwise}\end{cases}$ $\begin{cases}50(1-\mathrm{OVRLDtrack3}) & \text{if $(\mathrm{OVRLDtrack}> \mathrm{Logicthresh})\land (\mathrm{OVRLDtrack3}< \mathrm{Logicthresh})$}\\ -0.01\mathrm{OVRLDtrack3} & \text{otherwise}\end{cases}$ $\frac{1}{1+e^2\mathrm{APtrack2}}$ $1-\mathrm{RyRopen}$ $\frac{-1(\mathrm{i\_Na}+\mathrm{i\_CaNa}+\mathrm{i\_Na\_b}+\mathrm{i\_ns\_Na}+\mathrm{i\_NaCa}\times 3+\mathrm{i\_NaK}\times 3)\mathrm{A\_cap}}{\mathrm{V\_myo}F}$ $\frac{-1(\mathrm{I\_st}+\mathrm{i\_CaK}+\mathrm{i\_Kr}+\mathrm{i\_Ks}+\mathrm{i\_K1}+\mathrm{i\_Kp}+\mathrm{i\_K\_Na}+\mathrm{i\_K\_ATP}+\mathrm{i\_to}+\mathrm{i\_ns\_K}-\mathrm{i\_NaK}\times 2)\mathrm{A\_cap}}{\mathrm{V\_myo}F}$ $0.68\mathrm{volume}$ $0.0048\mathrm{volume}$ $0.0552\mathrm{volume}$