Model Mathematics

Component: Istim_parameters

Component: outputs

Component: g_parameters

Component: Membrane_potential

ddtimeV=-I_stim+I_ionCm I_ion=I_CaL+I_CaT+I_Kv+I_BK+I_Na+I_NCX+I_NaK+I_NsNa+I_NsK

Component: Cai

Cai_total_rate=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai_total=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai=Cai_total_rate1+n_CRTCRT_totalK_D_CRTCain_CRT-1Cain_CRT+K_D_CRT2+n_CaMCaM_totalK_D_CaMCain_CaM-1Cain_CaM+K_D_CaM2

Component: Nai

ddtimeNai=-I_Na+3I_NaK+3I_NCX+I_NsNaCvtFV_cell

Component: Ki

ddtimeKi=-I_Kv+I_BK+I_stim-2I_NaK+I_NsKCvtFV_cell

Component: Cai

Cai_total_rate=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai_total=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai=Cai_total_rate1+n_CRTCRT_totalK_D_CRTCain_CRT-1Cain_CRT+K_D_CRT2+n_CaMCaM_totalK_D_CaMCain_CaM-1Cain_CaM+K_D_CaM2

Component: Nai

ddtimeNai=-I_Na+3I_NaK+3I_NCX+I_NsNaCvtFV_cell

Component: Ki

ddtimeKi=-I_Kv+I_BK+I_stim-2I_NaK+I_NsKCvtFV_cell

Component: Cai

Cai_total_rate=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai_total=-I_CaL+I_CaT-2I_NCXCvt2FV_cell ddtimeCai=Cai_total_rate1+n_CRTCRT_totalK_D_CRTCain_CRT-1Cain_CRT+K_D_CRT2+n_CaMCaM_totalK_D_CaMCain_CaM-1Cain_CaM+K_D_CaM2

Component: Nai

ddtimeNai=-I_Na+3I_NaK+3I_NCX+I_NsNaCvtFV_cell

Component: Ki

ddtimeKi=-I_Kv+I_BK+I_stim-2I_NaK+I_NsKCvtFV_cell

Component: ICaL

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: ICaL_channel_states

alpha=phi0.731V30 beta=phi0.2149-V40 alpha_0=4alpha alpha_1=3alpha alpha_2=2alpha alpha_3=1alpha beta_0=1beta beta_1=2beta beta_2=3beta beta_3=4beta phi_f=phi0.4742V10 phi_s=phi0.005956-V40 xi_f=phi0.01407-V300 xi_s=phi0.01213V500 psi_f=phi0.02197V500 psi_s=phi0.00232-V280 omega_f=beta_3xi_fphi_falpha_3psi_f omega_s=beta_3xi_sphi_salpha_3psi_s omega_sf=xi_spsi_fxi_f omega_fs=psi_s delta=phi0.01 theta=phi41+1Cai norm=C3+C2+C1+C0+C3Ca+C2Ca+C1Ca+C0Ca+O_CaL+ICa+IVs+IVf+IVsCa+IVfCa ddtimeC3=alpha_2C2norm-alpha_3+beta_2+phi_s+phi_f+thetaC3norm+beta_3O_CaLnorm+omega_fIVfnorm+omega_sIVsnorm+deltaC3Canorm ddtimeC2=alpha_1C1norm-alpha_2+beta_1+thetaC2norm+beta_2C3norm+deltaC2Canorm ddtimeC1=alpha_0C0norm-alpha_1+beta_0+thetaC1norm+beta_1C2norm+deltaC1Canorm ddtimeC0=-alpha_0+thetaC0norm+beta_0C1norm+deltaC0Canorm ddtimeC3Ca=thetaC3norm+alpha_2C2Canorm-alpha_3+beta_2+phi_f+phi_s+deltaC3Canorm+beta_3ICanorm+omega_fIVfCanorm+omega_sIVsCanorm ddtimeC2Ca=thetaC2norm+alpha_1C1Canorm-alpha_2+beta_1+deltaC2Canorm+beta_2C3Canorm ddtimeC1Ca=thetaC1norm+alpha_0C0Canorm-alpha_1+beta_0+deltaC1Canorm+beta_1C2Canorm ddtimeC0Ca=thetaC0norm-alpha_0+deltaC0Canorm+beta_0C1Canorm ddtimeO_CaL=alpha_3C3norm-beta_3+psi_f+psi_s+thetaO_CaLnorm+xi_fIVfnorm+xi_sIVsnorm+deltaICanorm ddtimeICa=thetaO_CaLnorm+alpha_3C3Canorm-beta_3+psi_f+psi_s+deltaICanorm+xi_fIVfCanorm+xi_sIVsCanorm ddtimeIVs=phi_sC3norm+psi_sO_CaLnorm+omega_fsIVfnorm-omega_s+xi_s+omega_sf+thetaIVsnorm+deltaIVsCanorm ddtimeIVf=phi_fC3norm+psi_fO_CaLnorm-omega_f+xi_f+omega_fs+thetaIVfnorm+omega_sfIVsnorm+deltaIVfCanorm ddtimeIVsCa=thetaIVsnorm+phi_sC3Canorm+psi_sICanorm+omega_fsIVfCanorm-omega_s+xi_s+omega_sf+deltaIVsCanorm ddtimeIVfCa=thetaIVfnorm+phi_fC3Canorm+psi_fICanorm-omega_f+xi_f+omega_fs+deltaIVfCanorm+omega_sfIVsCanorm

Component: temperature_factor

phi=Q10T-T010

Component: ICaT

P_CaT=d_CaTf_CaT

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: gating_kinetics

ddtimeX=X_inf-Xtau_Xphi

Component: temperature_factor

phi=Q10T-T010

Component: gating_kinetics

ddtimeX=X_inf-Xtau_Xphi

Component: temperature_factor

phi=Q10T-T010

Component: d_inf

X_inf=11+-V+60.55.3

Component: f_inf

X_inf=11+V+75.54.0

Component: tau_f

tau=0.381178.6+14.7-V+5029001

Component: IKv

P_Kv=x_Kvy_Kv

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: gating_kinetics

ddtimeX=X_inf-Xtau_Xphi

Component: temperature_factor

phi=Q10T-T010

Component: gating_kinetics

ddtimeX=X_inf-Xtau_Xphi

Component: temperature_factor

phi=Q10T-T010

Component: x_inf

X_inf=11+-V+43.017.36

Component: y_inf

X_inf=11+V-44.912.0096

Component: IBK

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: IBK_channel_states

alpha=18.47188V1T beta=1-7.77556V1T k_C0O0=0.02162alpha k_C1O1=0.000869alpha k_C2O2=0.0000281alpha k_C3O3=0.000781alpha k_C4O4=0.044324alpha k_O0C0=318.1084beta k_O1C1=144.1736beta k_O2C2=32.6594beta k_O3C3=0.095312beta k_O4C4=0.000106beta k_C0C1=4k_onCai k_C1C2=3k_onCai k_C2C3=2k_onCai k_C3C4=1k_onCai k_C4C3=4k_off_C k_C3C2=3k_off_C k_C2C1=2k_off_C k_C1C0=1k_off_C k_O0O1=4k_onCai k_O1O2=3k_onCai k_O2O3=2k_onCai k_O3O4=1k_onCai k_O4O3=4k_off_O k_O3O2=3k_off_O k_O2O1=2k_off_O k_O1O0=1k_off_O norm=C0+C1+C2+C3+C4+O0+O1+O2+O3+O4 ddtimeC4=-k_C4C3+k_C4O4C4norm+k_C3C4C3norm+k_O4C4O4norm ddtimeC3=-k_C3C2+k_C3O3+k_C3C4C3norm+k_C2C3C2norm+k_O3C3O3norm+k_C4C3C4norm ddtimeC2=-k_C2C1+k_C2O2+k_C2C3C2norm+k_C1C2C1norm+k_O2C2O2norm+k_C3C2C3norm ddtimeC1=-k_C1C0+k_C1O1+k_C1C2C1norm+k_C0C1C0norm+k_O1C1O1norm+k_C2C1C2norm ddtimeC0=-k_C0C1+k_C0O0C0norm+k_C1C0C1norm+k_O0C0O0norm ddtimeO4=-k_O4O3+k_O4C4O4norm+k_O3O4O3norm+k_C4O4C4norm ddtimeO3=-k_O3O2+k_O3C3+k_O3O4O3norm+k_O2O3O2norm+k_C3O3C3norm+k_O4O3O4norm ddtimeO2=-k_O2O1+k_O2C2+k_O2O3O2norm+k_O1O2O1norm+k_C2O2C2norm+k_O3O2O3norm ddtimeO1=-k_O1O0+k_O1C1+k_O1O2O1norm+k_O0O1O0norm+k_C1O1C1norm+k_O2O1O2norm ddtimeO0=-k_O0O1+k_O0C0O0norm+k_O1O0O1norm+k_C0O0C0norm

Component: INa

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: INa_channel_states

k_I2I1=phi_T0.00392392.6793+0.0061468V k_I1O=phi_T0.12052-9.6028+0.083025V k_OC1=phi_T2.391-13.335-0.25289V k_C1C2=phi_T3.15660.36352+0.077193V k_C2C3=phi_T0.55432-0.099074+0.036441V k_C3C2=phi_T0.00052548-0.069102+0.0031945V k_C2C1=phi_T1.4496-0.1566+0.058353V k_C1O=phi_T1.53290.0093193+0.041075V k_OI1=phi_T1.61640.30763+0.0060535V k_I1I2=phi_T0.0277350.05149-0.046865V k_I1C1=phi_T1.9046-2.484+0.020406V k_C1I1=phi_T0.00021688-0.063438+0.0046683V norm=O_Na+C1+C2+C3+I1+I2 ddtimeC3=-k_C3C2C3norm+k_C2C3C2norm ddtimeC2=-k_C2C1+k_C2C3C2norm+k_C1C2C1norm+k_C3C2C3norm ddtimeC1=-k_C1C2+k_C1O+k_C1I1C1norm+k_OC1O_Nanorm+k_C2C1C2norm+k_I1C1I1norm ddtimeO_Na=-k_OC1+k_OI1O_Nanorm+k_C1OC1norm+k_I1OI1norm ddtimeI2=-k_I2I1I2norm+k_I1I2I1norm ddtimeI1=-k_I1O+k_I1I2+k_I1C1I1norm+k_I2I1I2norm+k_C1I1C1norm+k_OI1O_Nanorm

Component: temperature_factor

phi=Q10T-T010

Component: INCX

I_NCX=P_NCXgammaVFRTNai3Cao-2.5gamma-1VFRTNao3CaiK_mNai3+Nao3K_mCa+Cao1+k_satgamma-1VFRT

Component: INaK

I_NaK=P_NaKKoNaiKo+K_mKK_mNa+Nai1+0.1245-0.1VFRT+0.0353-VFRT

Component: INS_Na

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: INS_K

Component: Ionic_currents

I=g_maxPOV-E

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Nernst_potential

E=RTzFlnX_oX_i

Component: Periodic_IStim_protocol

I_stim=GcoupleV-V_ICC t=timemodt_period V_ICC=V_ICCrest+V_ICCamptf_2ift<t_ICCpeakV_ICCrest+V_ICCamp1+-f_12t_slope1+t-f_2-0.5f_1t_slopeotherwise

Component: time

Component: constants

Component: model_parameters

Component: initial_conditions