Model Mathematics

Component: Environment

FonRT=FRT

Component: cell

past=timestim_periodstim_period i_Stim=stim_amplitudeiftime-paststim_offsettime-paststim_offset+stim_duration0otherwise Itotal=-xina+xik1+xikr+xiks+xito+xiNaCa+xica+xiNaK+i_Stim ddtimeV=Itotal

Component: INa

am=0.321V+47.131--0.1V+47.13if|V+47.13|>0.0013.2otherwise bm=0.08-V11 ah=0.13580+V-6.8ifV<-400otherwise bh=3.560.079V+3100000.35VifV<-4010.131+V+10.66-11.1otherwise aj=-1271400.2444V-0.00003474-0.04391V1V+37.781+0.311V+79.23ifV<-400otherwise bj=0.1212-0.01052V1+-0.1378V+40.14ifV<-400.3-0.0000002535V1+-0.1V+32otherwise ddtimexh=ah1-xh-bhxh ddtimexj=aj1-xj-bjxj ddtimexm=am1-xm-bmxm xina=gnaxhxjxmxmxmV-ena

Component: ICaL

za=V2FonRT rxa=4pcaFFonRTcsmza-0.341Ca_o2FonRTif|za|<0.0014pcaVFFonRTcsmza-0.341Ca_oza-1otherwise poinf=11+-V-vths6 alpha=poinftaupo beta=1-poinftaupo fca=11+catCa_dyad3 s1=0.0182688fca k1=0.024168fca s2=s1r1r2k2k1 s2t=s1tr1r2k2tk1t poi=11+-V-vxsx k3=1-poitau3 k3t=k3 Pr=1-11+-V-vysy recov=10+4954V15.6 tau_ca=tca1+Ca_dyadcpt4+0.1 tauca=recov-tau_caPr+tau_ca tauba=recov-450Pr+450 Ps=11+-V-vyrsyr k6=fcaPstauca k5=1-Pstauca k6t=Pstauba k5t=1-Pstauba k4=k3alphabetak1k2k5k6 k4t=k3talphabetak1tk2tk5tk6t po=1-xi1ca-xi2ca-xi1ba-xi2ba-c1-c2 ddtimec1=alphac2+k2xi1ca+k2txi1ba+r2po-beta+r1+k1t+k1c1 ddtimec2=betac1+k5xi2ca+k5txi2ba-k6+k6t+alphac2 ddtimexi1ca=k1c1+k4xi2ca+s1po-k3+k2+s2xi1ca ddtimexi1ba=k1tc1+k4txi2ba+s1tpo-k3t+k2t+s2txi1ba ddtimexi2ca=k3xi1ca+k6c2-k5+k4xi2ca ddtimexi2ba=k3txi1ba+k6tc2-k5t+k4txi2ba jca=gcaporxa xica=2wcajca

Component: IK1

aki=1.021+0.2385V-ek-59.215 bki=0.491240.08032V-ek+5.476+10.06175V-ek-594.311+-0.5143V-ek+4.753 xkin=akiaki+bki xik1=gkixK_o5.4xkinV-ek

Component: IKr

xkrv1=0.001381V+71--0.123V+7if|V+7|>0.0010.001380.123otherwise xkrv2=0.000611V+100.145V+10-1if|V+10|>0.0010.000610.145otherwise taukr=1xkrv1+xkrv2 xkrinf=11+-V+507.5 rg=11+V+3322.4 xikr=gkrK_o5.4xrrgV-ek ddtimexr=xkrinf-xrtaukr

Component: IKs

xs1ss=11+-V-1.516.7 xs2ss=xs1ss tauxs1=10.00007190.148+0.0001310.0687if|V+30|<0.0010.068710.0000719V+301--0.148V+30+0.000131V+300.0687V+30-1otherwise tauxs2=4tauxs1 gksx=1+0.81+0.5Ca_i3 xiks=gksgksxxs1xs2V-eks ddtimexs1=xs1ss-xs1tauxs1 ddtimexs2=xs2ss-xs2tauxs2

Component: Ito

rt1=-V+315 rt2=V+33.510 rt3=V+6010 rt4=-V30V30 rt5=V+33.510 xtos_inf=11+rt1 ytos_inf=11+rt2 xtof_inf=xtos_inf ytof_inf=ytos_inf rs_inf=11+rt2 txs=91+-rt1+0.5 tys=30001+rt3+30 txf=3.5rt4+1.5 tyf=201+rt5+20 ddtimextos=xtos_inf-xtostxs ddtimeytos=ytos_inf-ytostys ddtimextof=xtof_inf-xtoftxf ddtimeytof=ytof_inf-ytoftyf xitos=gtosxtosytos+0.5rs_infV-ek xitof=gtofxtofytofV-ek xito=xitos+xitof

Component: INaK

sigma=Na_o67.3-17 fNaK=11+0.1245-0.1VFonRT+0.0365sigma-VFonRT xiNaK=gNaKfNaKNa_iNa_i+xkmnaiK_oK_o+xkmko

Component: INaCa

zw3=Na_i3Ca_oV0.35FonRT-Na_o3csmV0.35-1FonRT zw4=1+0.2V0.35-1FonRT aloss=11+xkdnaCa_submem3 yz1=xmcaoNa_i3+xmnao3csm yz2=xmnai3Ca_o1+csmxmcai yz3=xmcaiNa_o31+Na_ixmnai3 yz4=Na_i3Ca_o+Na_o3csm zw8=yz1+yz2+yz3+yz4 jNaCa=gNaCaalosszw3zw4zw8 xiNaCa=wcajNaCa

Component: Irel

ddtimeCa_JSR=Ca_NSR-Ca_JSRtaua bv=1-avcstar-50 Qr0=Ca_JSR-501ifCa_JSR>50Ca_JSR<cstaravCa_JSR+bvifCa_JSRcstar0otherwise Qr=Ca_NSRQr0cstar sparkV=-ayV+301+-ayV+30 spark_rate=gryr1po|rxa|sparkV ddtimexir=spark_rateQr-xir1-taurdCa_JSRCa_NSRtaur xirp=poQr|rxa|gbarsr1-axV+301+-axV+30 xicap=pogdyad|rxa| xiryr=xirp+xicap

Component: Ileak_Iup_Ixfer

jup=vupCa_iCa_iCa_iCa_i+cupcup jleak=gleakCa_NSRCa_NSRCa_NSRCa_NSR+kjkjCa_NSR16.667-Ca_i

Component: Na

ddtimeNa_i=-xina+3xiNaK+3xiNaCawca1000

Component: Ca

bpxs=bcalxkcalxkcal+Ca_submemxkcal+Ca_submem spxs=srmaxsrkdsrkd+Ca_submemsrkd+Ca_submem mempxs=bmemkmemkmem+Ca_submemkmem+Ca_submem sarpxs=bsarksarksar+Ca_submemksar+Ca_submem dcsib=11+bpxs+spxs+mempxs+sarpxs bpxi=bcalxkcalxkcal+Ca_ixkcal+Ca_i spxi=srmaxsrkdsrkd+Ca_isrkd+Ca_i mempxi=bmemkmemkmem+Ca_ikmem+Ca_i sarpxi=bsarksarksar+Ca_iksar+Ca_i dciib=11+bpxi+spxi+mempxi+sarpxi jd=Ca_submem-Ca_itaud xbi=xkonCa_ibtrop-tropi-xkofftropi xbs=xkonCa_submembtrop-trops-xkofftrops dCa_JSR=-xir+jup-jleak ddtimeCa_dyad=xiryr-Ca_dyad-Ca_submemtaups ddtimeCa_submem=dcsib50xir-jd-jca+jNaCa-xbs ddtimeCa_i=dciibjd-jup+jleak-xbi ddtimeCa_NSR=dCa_JSR ddtimetropi=xbi ddtimetrops=xbs csm=Ca_submem1000

Component: reversal_potentials

ek=1FonRTlnK_oK_i eks=1FonRTlnK_o+prNaKNa_oK_i+prNaKNa_i ena=1FonRTlnNa_oNa_i