$\frac{-1}{\mathrm{Cm}}(\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})$ $\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}$ $\frac{-1}{\mathrm{Cm}}(\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\_2}+\mathrm{beta\_13}\mathrm{P\_O\_Na}+\mathrm{beta\_2}\mathrm{P\_IF}+\mathrm{alpha\_13}\mathrm{P\_C1}$ $\mathrm{beta\_12}+\mathrm{alpha\_13}+\mathrm{beta\_3}\mathrm{P\_C1}+\mathrm{alpha\_12}\mathrm{P\_C2}+\mathrm{beta\_13}\mathrm{P\_O\_Na}+\mathrm{alpha\_3}\mathrm{P\_IF}$ $\mathrm{beta\_11}+\mathrm{alpha\_12}\mathrm{P\_C2}+\mathrm{alpha\_11}\mathrm{P\_C3}+\mathrm{beta\_12}\mathrm{P\_C1}$ $\mathrm{beta\_11}\mathrm{P\_C2}-\mathrm{alpha\_11}\mathrm{P\_C3}$ $\mathrm{beta\_2}+\mathrm{alpha\_3}+\mathrm{alpha\_4}\mathrm{P\_IF}+\mathrm{beta\_3}\mathrm{P\_C1}+\mathrm{beta\_4}\mathrm{P\_IS}+\mathrm{alpha\_2}\mathrm{P\_O\_Na}$ $\mathrm{alpha\_4}\mathrm{P\_IF}-\mathrm{beta\_4}\mathrm{P\_IS}$ $\frac{3802}{0.1027e^{\frac{-V}{17}}+0.2e^{\frac{-V}{150}}}$ $191.7e^{\frac{-V}{20.3}}$ $\frac{3802}{0.1027e^{\frac{-V}{15}}+0.23e^{\frac{-V}{150}}}$ $200e^{\frac{-(V-5)}{20.3}}$ $\frac{3802}{0.1027e^{\frac{-V}{12}}+0.25e^{\frac{-V}{150}}}$ $220e^{\frac{-(V-10)}{20.3}}$ $9178e^{\frac{V}{29.68}}$ $\frac{\mathrm{alpha\_13}\mathrm{alpha\_2}\mathrm{alpha\_3}}{\mathrm{beta\_13}\mathrm{beta\_3}}$ $3.7933E-7e^{\frac{-V}{5.2}}$ $8.4+0.02V$ $\frac{\mathrm{alpha\_2}}{100}$ $\mathrm{alpha\_3}$ $\mathrm{g\_Na}\mathrm{P\_O\_Na}(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{0.001}{0.035\times 6.24} & \text{if $\left|\mathrm{E0\_d}\right|< 1E-5$}\\ \frac{0.001\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+32}{8}}}+\frac{0.6}{1+e^{\frac{50-V}{20}}}$ $\frac{0.001}{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}}}$ $0.0037+\frac{0.0061}{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.000875V+0.012 & \text{if $V\le 0$}\\ 0.012 & \text{otherwise}\end{cases}$ $\mathrm{g\_CaT}bbg(V-\mathrm{E\_Ca})$ $\mathrm{beta\_beta}+\mathrm{beta\_i}\mathrm{P\_O}+\mathrm{alpha\_alpha}\mathrm{P\_C1}+\mathrm{alpha\_i}\mathrm{P\_I}$ $\mathrm{beta\_in}+\mathrm{alpha\_alpha}+\mathrm{alpha\_alpha}\mathrm{P\_C1}+\mathrm{alpha\_in}\mathrm{P\_C2}+\mathrm{beta\_beta}\mathrm{P\_O}+\mathrm{mu}\mathrm{P\_I}$ $\mathrm{beta}+\mathrm{alpha\_in}\mathrm{P\_C2}+\mathrm{alpha}\mathrm{P\_C3}+\mathrm{beta\_in}\mathrm{P\_C1}$ $\mathrm{beta}\mathrm{P\_C2}-\mathrm{alpha}\mathrm{P\_C3}$ $\mathrm{mu}+\mathrm{alpha\_i}\mathrm{P\_I}+\mathrm{alpha\_alpha}\mathrm{P\_C1}+\mathrm{beta\_i}\mathrm{P\_O}$ $55.5e^{0.05547153(V-12)}$ $2.357e^{-0.036588V}\times 10.5$ $65.5e^{0.05547153(V-36)}$ $2.9357e^{-0.02158V}\times 6.3$ $\frac{439e^{-0.02352(V+25)}\times 4.5}{\mathrm{Ko}}$ $656e^{0.000942V}\left(\frac{4.5}{\mathrm{Ko}}\right)^{0.3}$ $\frac{\mathrm{alpha\_i}\mathrm{beta\_beta}\mathrm{alpha\_alpha}}{\mathrm{alpha\_alpha}\mathrm{beta\_i}}$ $\mathrm{g\_Kr}\mathrm{P\_O}(V-\mathrm{E\_Kr})$ $2.7\times 0.0135\mathrm{Ko}^{0.59}$ $\frac{RT}{F}\ln \left(\frac{\mathrm{Ko}}{\mathrm{Ki}}\right)$ $\frac{\mathrm{xs1\_infinity}-\mathrm{xs1}}{\mathrm{tau\_xs1}}$ $\frac{1}{1+e^{\frac{-(V-1.5)}{16.7}}}$ $\frac{0.001}{\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\times 0.001}{\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})$ $0.433(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{1020}{1+e^{0.2385(V-\mathrm{E\_K}-59.215)}}$ $\frac{1000(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)$ $0.75\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})$ $0\times 0.12848$ $\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{10000e^{\frac{V-40}{25}}}{1+e^{\frac{V-40}{25}}}$ $\frac{10000e^{\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{15}{1+e^{\frac{V+60}{5}}}$ $\frac{100e^{\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}$ $0\times 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}}$ $\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{-\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}$ $\frac{0.0048}{0.68}\mathrm{V\_myo}$ $\frac{0.0552}{0.68}\mathrm{V\_myo}$ $\begin{cases}100000(1-\mathrm{APtrack})-500\mathrm{APtrack} & \text{if $\mathrm{dVdt}> 150000$}\\ -500\mathrm{APtrack} & \text{otherwise}\end{cases}$ $\begin{cases}100000(1-\mathrm{APtrack2})-500\mathrm{APtrack2} & \text{if $(\mathrm{APtrack}< 0.2)\land (\mathrm{APtrack}> 0.18)$}\\ -500\mathrm{APtrack2} & \text{otherwise}\end{cases}$ $\begin{cases}100000(1-\mathrm{APtrack3})-500\mathrm{APtrack3} & \text{if $(\mathrm{APtrack}< 0.2)\land (\mathrm{APtrack}> 0.18)$}\\ -10\mathrm{APtrack3} & \text{otherwise}\end{cases}$ $\begin{cases}\frac{-\mathrm{A\_cap}(\mathrm{i\_CaCa}+\mathrm{i\_Ca\_T}-\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)$}\\ -500\mathrm{Cainfluxtrack} & \text{otherwise}\end{cases}$ $\begin{cases}50000(1-\mathrm{OVRLDtrack}) & \text{if $(\frac{1}{1+\frac{\mathrm{K\_mCSQN}}{\mathrm{Ca\_NSR}}}> \mathrm{CSQNthresh})\land (\mathrm{OVRLDtrack3}< 0.37)\land (\mathrm{APtrack3}< 0.37)$}\\ -500\mathrm{OVRLDtrack} & \text{otherwise}\end{cases}$ $\begin{cases}50000(1-\mathrm{OVRLDtrack2}) & \text{if $(\mathrm{OVRLDtrack}> \mathrm{Logicthresh})\land (\mathrm{OVRLDtrack2}< \mathrm{Logicthresh})$}\\ -500\mathrm{OVRLDtrack2} & \text{otherwise}\end{cases}$ $\begin{cases}50000(1-\mathrm{OVRLDtrack3}) & \text{if $(\mathrm{OVRLDtrack}> \mathrm{Logicthresh})\land (\mathrm{OVRLDtrack3}< \mathrm{Logicthresh})$}\\ -10\mathrm{OVRLDtrack3} & \text{otherwise}\end{cases}$ $\frac{-(\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{-(\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}$