Model Mathematics

Component: environment

Component: membrane

i_Stim=stim_amplitudeiftimestim_starttimestim_endtime-stim_start-time-stim_startstim_periodstim_periodstim_duration0otherwise RTONF=RTF ddtimeV=-i_Na+i_CaL+i_to+i_Kr+i_f+i_st+i_K1+i_NaCa+i_p+i_b+i_ACh+i_StimC

Component: hyperpolarising_activated_current

i_f=yg_fV--30

Component: hyperpolarising_activated_current_y_gate

y_inf=11+V+83.19--7.2ACh0.691.26 e -50.69+ACh0.6913.56 tau_y=0.25+2-V+702500 ddtimey=y_inf-ytau_y

Component: rapid_delayed_rectifier_potassium_current

E_K=RTONFlnKcKi i_Kr=g_Kr0.9paf+0.1paspikV-E_K

Component: rapid_delayed_rectifier_potassium_current_paf_gate

paf_infinity=11+V+10.22-8.5 tau_paf=1170.0398V+0.211-0.051V ddtimepaf=paf_infinity-paftau_paf

Component: rapid_delayed_rectifier_potassium_current_pas_gate

pas_infinity=11+V+10.22-8.5 tau_pas=0.33581+0.90673-V+102988.05 ddtimepas=pas_infinity-pastau_pas

Component: rapid_delayed_rectifier_potassium_current_pik_gate

pik_infinity=11+V+4.915.141-0.3-V2500 alpha_pik=92.01-0.0183V beta_pik=603.60.00942V tau_pik=1alpha_pik+beta_pik ddtimepik=pik_infinity-piktau_pik

Component: time_independent_potassium_current

g_K1_prime=g_K10.5+0.51+V+305 i_K1=g_K1_primeKcKc+0.593V+81.91+1.393V+81.9+3.6RTONF

Component: background_current

i_b=g_bV-E_b

Component: sodium_potassium_pump

i_p=I_pNai5.64+Nai3Kc0.621+Kc21.61.5+-V+6040

Component: sodium_calcium_exchange_current

do=1+CaoKco1+QcoVRTONF+NaoK1no+Nao2K1noK2no+Nao3K1noK2noK3no k32=QnV2RTONF k23=Nao2K1noK2no+Nao3K1noK2noK3no-QnV2RTONFdo k21=CaoKco-QcoVRTONFdo k34=NaoK3no+Nao di=1+CasubKci1+-QciVRTONF+NaiKcni+NaiK1ni+Nai2K1niK2ni+Nai3K1niK2niK3ni k41=-QnV2RTONF k14=Nai2K1niK2ni+Nai3K1niK2niK3niQnV2RTONFdi k12=CasubKci-QciVRTONFdi k43=NaiK3ni+Nai x1=k41k34k23+k21+k21k32k43+k41 x2=k32k43k14+k12+k41k12k34+k32 x3=k14k43k23+k21+k12k23k43+k41 x4=k23k34k14+k12+k14k21k34+k32 i_NaCa=kNaCax2k21-x1k12x1+x2+x3+x4

Component: fast_sodium_current

E_Na=RTONFlnNaoNai i_Na=g_Nam30.635h1+0.365h2NaoVFRTONFV-E_NaRTONF-1VRTONF-1

Component: fast_sodium_current_m_gate

E0_m=V+44.4 alpha_m=-460-12.673E0_m-12.673if|E0_m|<delta_m-460E0_mE0_m-12.673-1otherwise beta_m=18400E0_m-12.673 ddtimem=alpha_m1-m-beta_mm

Component: fast_sodium_current_h1_gate

alpha_h1=44.9V+66.9-5.57 beta_h1=14911+323.3V+94.6-12.9 h1_inf=alpha_h1alpha_h1+beta_h1 tau_h1=0.031+V+406+0.00035 ddtimeh1=h1_inf-h1tau_h1

Component: fast_sodium_current_h2_gate

alpha_h2=44.9V+66.9-5.57 beta_h2=14911+323.3V+94.6-12.9 h2_inf=alpha_h2alpha_h2+beta_h2 tau_h2=0.121+V+602+0.00295 ddtimeh2=h2_inf-h2tau_h2

Component: L_type_calcium_current

i_CaL=g_CaLd0.675f+0.325f2V-E_CaL1-i_AChACh0.9 e -4+ACh1

Component: L_type_calcium_current_d_gate

alpha_d=-26.12V+35V+35-2.5-1+-78.11V-0.208V-1 beta_d=10.52V-50.4V-5-1 d_inf=11+V--3.2+act_shiftslope_factor_act tau_d=1alpha_d+beta_d ddtimed=d_inf-dtau_d

Component: L_type_calcium_current_f_gate

f_inf=11+V--24+inact_shift6.31 tau_f=0.01+0.1539-V+402185.67 ddtimef=f_inf-ftau_f

Component: L_type_calcium_current_f2_gate

f2_inf=11+V--24+inact_shift6.31 tau_f2=0.06+0.480762.25-V--402138.04 ddtimef2=f2_inf-f2tau_f2

Component: transient_outward_potassium_current

i_to=g_tor0.45q_fast+0.55q_slowV-E_K E_K=RTONFlnKcKi

Component: transient_outward_potassium_current_r_gate

ddtimer=r_infinity-rtau_r r_infinity=11+V-7.44-16.4 tau_r=0.596 e -3+3.118 e -31.0370.09V+30.61+0.396-0.12V+23.84

Component: transient_outward_potassium_current_qfast_gate

ddtimeq_fast=qfast_infinity-q_fasttau_qfast qfast_infinity=11+V+33.86.12 tau_qfast=12.66 e -3+4.727161+V+154.523.96

Component: transient_outward_potassium_current_qslow_gate

ddtimeq_slow=qslow_infinity-q_slowtau_qslow qslow_infinity=11+V+33.86.12 tau_qslow=0.1+4-V+652500

Component: sustained_outward_potassium_current

i_st=g_stqaqiV-E_st

Component: sustained_outward_potassium_current_qa_gate

ddtimeqa=qa_infinity-qatau_qa qa_infinity=11+V--49.1-8.98 alpha_qa=10.15-V11+0.2-V700 beta_qa=116V8+15V50 tau_qa=0.001alpha_qa+beta_qa

Component: sustained_outward_potassium_current_qi_gate

ddtimeqi=qi_infinity-qitau_qi alpha_qi=0.15043100V13+700V70 beta_qi=0.150495-V10+50-V700+0.0002291+-V5 qi_infinity=alpha_qialpha_qi+beta_qi tau_qi=0.001alpha_qi+beta_qi

Component: acetylcholine_sensitive_current

g_ACh=g_ACh_maxachfachsACh1.5K_ACh1.5+ACh1.5 i_ACh=g_AChKc10+KcV-E_K1+V-E_K-1402.5RTONF

Component: acetylcholine_sensitive_current_achf_gate

beta_achf=1201+-V+5015 ddtimeachf=alpha_achf1-achf-beta_achfachf

Component: acetylcholine_sensitive_current_achs_gate

beta_achs=5.821+-V+5015 ddtimeachs=alpha_achs1-achs-beta_achsachs

Component: extracellular_sodium_concentration

Component: intracellular_sodium_concentration

Component: extracellular_calcium_concentration

Component: intracellular_calcium_concentration

i_diff=Casub-Cai0.04 e -3 i_up=51+K_upCai i_tr=Ca_up-Ca_reltau_tr i_rel=P_relCa_rel-Casub1+0.0012Casub2 ddtimeCa_up=i_up-i_trV_relV_up V_up=0.0116V_cell ddtimeCa_rel=i_tr-i_rel-10diff_f_CQ V_rel=0.0012V_cell ddtimeCai=i_diffV_sub-i_upV_upVi-0.045diff_f_CMi+0.031diff_f_TC+0.062diff_f_TMC ddtimeCasub=-i_CaL-2i_NaCa2F+i_relV_relV_sub-i_diff-0.045diff_f_CMs-0.0311.2diff_f_CSL Vi=0.46V_cell-V_sub V_sub=0.01V_cell diff_f_TC=88.8 e 3Cai1-f_TC-0.446 e 3f_TC ddtimef_TC=diff_f_TC diff_f_TMC=227.7 e 3Cai1-f_TMC-f_TMM-0.00751 e 3f_TMC ddtimef_TMC=diff_f_TMC diff_f_TMM=2.277 e 32.51-f_TMC-f_TMM-0.751 e 3f_TMM ddtimef_TMM=diff_f_TMM diff_f_CMi=227.7 e 3Cai1-f_CMi-0.542 e 3f_CMi ddtimef_CMi=diff_f_CMi diff_f_CMs=227.7 e 3Casub1-f_CMs-0.542 e 3f_CMs ddtimef_CMs=diff_f_CMs diff_f_CQ=0.534 e 3Ca_rel1-f_CQ-0.445 e 3f_CQ ddtimef_CQ=diff_f_CQ diff_f_CSL=0.001115Casub1-f_CSL-1000f_CSL ddtimef_CSL=diff_f_CSL

Component: extracellular_potassium_concentration

Component: intracellular_potassium_concentration