# Model Mathematics

### Component: cell

$i_Stim=-stim_amplitudeiftime-⌊timestim_period⌋⁢stim_period≥stim_start∧time-⌊timestim_period⌋⁢stim_period≦stim_start+stim_duration0otherwiseddtimeV=-i_Na+i_Nab+i_NaK+i_Kr+i_Ks+i_tos+i_tof+i_K1+i_Kp+i_NaCa+i_Cl_Ca+i_Clb+i_CaL+i_Cab+i_Cap+i_Stim$

### Component: model_parameters

$Cm=Cm_per_area⁢2⁢cell_radius10000⁢π⁢cell_length10000Vol_Cell=3.141592654⁢cell_radius10002⁢cell_length10003Vol_SR=0.035⁢Vol_CellVol_SL=0.02⁢Vol_CellVol_jct=0.00051⁢Vol_CellVol_cytosol=0.65⁢Vol_CellVol_mito=Vol_Cell-Vol_SR+Vol_SL+Vol_jct+Vol_cytosol$

### Component: reversal_potentials

$E_Na_jct=R⁢TF⁢ln⁡NaoNa_jctE_Na_SL=R⁢TF⁢ln⁡NaoNa_SLE_Ca_jct=R⁢T2⁢F⁢ln⁡CaoCa_jctE_Ca_SL=R⁢T2⁢F⁢ln⁡CaoCa_SLE_K=R⁢TF⁢ln⁡KoKiE_Cl=-R⁢TF⁢ln⁡CloCli$

### Component: INa

$openProb=m3⁢h⁢ji_Na_jct=Fx_Na_jct⁢G_INa⁢openProb⁢V-E_Na_jcti_Na_SL=Fx_Na_SL⁢G_INa⁢openProb⁢V-E_Na_SLi_Na=i_Na_jct+i_Na_SL$

### Component: INa_h_gate

$alpha_h=0.135⁢ⅇ80+V-6.8ifV<-400otherwisebeta_h=3.56⁢ⅇ0.079⁢V+3.1e5⁢ⅇ0.35⁢VifV<-4010.13⁢1+ⅇV+10.66-11.1otherwisetau_h=1alpha_h+beta_hh_infinity=alpha_halpha_h+beta_hddtimeh=h_infinity-htau_h$

### Component: INa_j_gate

$alpha_j=-1.2714e5⁢ⅇ0.2444⁢V-3.474e-5⁢ⅇ-0.04391⁢V⁢V+37.7811+ⅇ0.311⁢V+79.23ifV<-400otherwisebeta_j=0.1212⁢ⅇ-0.01052⁢V1+ⅇ-0.1378⁢V+40.14ifV<-400.3⁢ⅇ-2.535e-7⁢V1+ⅇ-0.1⁢V+32otherwisetau_j=1alpha_j+beta_jj_infinity=alpha_jalpha_j+beta_jddtimej=j_infinity-jtau_j$

### Component: INa_m_gate

$alpha_m=0.32⁢V+47.1311-ⅇ-0.1⁢V+47.13beta_m=0.08⁢ⅇ-V11tau_m=1alpha_m+beta_mm_infinity=alpha_malpha_m+beta_mddtimem=m_infinity-mtau_m$

### Component: INab

$i_Nab_jct=Fx_NaBk_jct⁢G_NaBk⁢V-E_Na_jcti_Nab_SL=Fx_NaBk_SL⁢G_NaBk⁢V-E_Na_SLi_Nab=i_Nab_jct+i_Nab_SL$

### Component: INaK

$Q_NaK=Q10_NaKT-31010Q_Km_Nai=Q10_Km_NaiT-31010sigma=ⅇNao67.3-17f_NaK=11+0.1245⁢ⅇ-0.1⁢V⁢FR⁢T+0.0365⁢sigma⁢ⅇ-V⁢FR⁢Ti_NaK_jct=Fx_NaK_jct⁢Q_NaK⁢I_NaK_max⁢f_NaK1+Q_Km_Nai⁢Km_NaiNa_jctH_NaK⁢KoKo+Km_Koi_NaK_SL=Fx_NaK_SL⁢Q_NaK⁢I_NaK_max⁢f_NaK1+Q_Km_Nai⁢Km_NaiNa_SLH_NaK⁢KoKo+Km_Koi_NaK=i_NaK_jct+i_NaK_SL$

### Component: IKr

$G_IKr=0.03⁢Ko5.4i_Kr=G_IKr⁢Xr⁢Rr⁢V-E_K$

### Component: IKr_Xr_gate

$ddtimeXr=Xr_infinity-Xrtau_XrXr_infinity=11+ⅇ-50+V7.5tau_Xr=10.00138⁢V+71-ⅇ-0.123⁢V+7+0.00061⁢V+10ⅇ0.145⁢V+10-1$

### Component: IKr_Rr_gate

$Rr=11+ⅇ33+V22.4$

### Component: IKs

$pCa_jct=-log10⁡Ca_jct1+3pCa_SL=-log10⁡Ca_SL1+3G_Ks_jct=0.07⁢0.057+0.191+ⅇ-7.2+pCa_jct0.6G_Ks_SL=0.07⁢0.057+0.191+ⅇ-7.2+pCa_SL0.6E_Ks_jct=R⁢TF⁢ln⁡Ko+pKNa⁢NaoKi+pKNa⁢Na_jctE_Ks_SL=R⁢TF⁢ln⁡Ko+pKNa⁢NaoKi+pKNa⁢Na_SLE_Ks=R⁢TF⁢ln⁡Ko+pKNa⁢NaoKi+pKNa⁢Naii_Ks_jct=Fx_Ks_jct⁢G_Ks_jct⁢Xs2⁢V-E_Ksi_Ks_SL=Fx_Ks_SL⁢G_Ks_SL⁢Xs2⁢V-E_Ksi_Ks=i_Ks_jct+i_Ks_SL$

### Component: IKs_Xs_gate

$ddtimeXs=Xs_infinity-Xstau_XsXs_infinity=11+ⅇ-V-1.516.7tau_Xs=17.19e-5⁢V+301-ⅇ-0.148⁢V+30+1.31e-4⁢V+30-1+ⅇ0.0687⁢V+30$

### Component: Itos

$i_tos=G_tos⁢X_tos⁢Y_tos+0.5⁢R_tos⁢V-E_K$

### Component: Itos_X_gate

$ddtimeX_tos=X_tos_infinity-X_tostau_X_tosX_tos_infinity=11+ⅇ-V+315tau_X_tos=91+ⅇV+315+0.5$

### Component: Itos_Y_gate

$ddtimeY_tos=Y_tos_infinity-Y_tostau_Y_tosY_tos_infinity=11+ⅇV+33.510tau_Y_tos=30001+ⅇV+6010+30$

### Component: Itos_R_gate

$R_tos_infinity=11+ⅇV+33.510tau_R_tos=2.8e31+ⅇV+6010+220R_tos=R_tos_infinityddtimeR_tos_other=R_tos_infinity-R_tostau_R_tos$

### Component: Itof

$i_tof=G_tof⁢X_tof⁢Y_tof⁢V-E_K$

### Component: Itof_X_gate

$ddtimeX_tof=X_tof_infinity-X_toftau_X_tofX_tof_infinity=11+ⅇ-V+315tau_X_tof=3.5⁢ⅇ-V302+1.5$

### Component: Itof_Y_gate

$ddtimeY_tof=Y_tof_infinity-Y_toftau_Y_tofY_tof_infinity=11+ⅇV+33.510tau_Y_tof=201+ⅇV+33.510+20$

### Component: IK1

$G_K1=0.9⁢Ko5.4i_K1=G_K1⁢K1_infinity⁢V-E_K$

### Component: IK1_K1_gate

$K1_infinity=alpha_K1alpha_K1+beta_K1alpha_K1=1.021+ⅇ0.2385⁢V-E_K+59.215beta_K1=0.49124⁢ⅇ0.08032⁢V-E_K+5.476+1⁢ⅇ0.06175⁢V-E_K+594.311+ⅇ-0.5143⁢V-E_K+4.753$

### Component: IKp

$i_Kp=g_Kp⁢V-E_K1+ⅇ7.488-V5.98$

### Component: ICl_Ca

$i_Cl_Ca=G_Cl⁢V-E_Cl⁢Fx_Cl_jct1+Kd_ClCaCa_jct+Fx_Cl_SL1+Kd_ClCaCa_SL$

### Component: IClb

$i_Clb=G_ClBk⁢V-E_Cl$

### Component: ICaL

$Q_CaL=Q10_CaLT-31010temp=d⁢f⁢Q_CaL⁢V⁢F2R⁢Ti_CaL_Ca_jct=temp⁢fCa_jct⁢Fx_ICaL_jct⁢PCa⁢4⁢gamma_Cai⁢Ca_jct⁢ⅇ2⁢V⁢FR⁢T-gamma_Cao⁢Caoⅇ2⁢V⁢FR⁢T-1i_CaL_Na_jct=temp⁢fCa_jct⁢Fx_ICaL_jct⁢PNa⁢gamma_Nai⁢Na_jct⁢ⅇV⁢FR⁢T-gamma_Nao⁢NaoⅇV⁢FR⁢T-1i_CaL_Ca_SL=temp⁢fCa_SL⁢Fx_ICaL_SL⁢PCa⁢4⁢gamma_Cai⁢Ca_SL⁢ⅇ2⁢V⁢FR⁢T-gamma_Cao⁢Caoⅇ2⁢V⁢FR⁢T-1i_CaL_Na_SL=temp⁢fCa_SL⁢Fx_ICaL_SL⁢PNa⁢gamma_Nai⁢Na_SL⁢ⅇV⁢FR⁢T-gamma_Nao⁢NaoⅇV⁢FR⁢T-1i_CaL_K=temp⁢fCa_SL⁢Fx_ICaL_SL+fCa_jct⁢Fx_ICaL_jct⁢PK⁢gamma_Ki⁢Ki⁢ⅇV⁢FR⁢T-gamma_Ko⁢KoⅇV⁢FR⁢T-1i_CaL=i_CaL_Ca_SL+i_CaL_Ca_jct+i_CaL_Na_SL+i_CaL_Na_jct+i_CaL_K$

### Component: ICaL_d_gate

$d_infinity=11+ⅇ-V+14.56tau_d=1⁢d_infinity⁢1-ⅇ-V+14.560.035⁢V+14.5ddtimed=d_infinity-dtau_d$

### Component: ICaL_f_gate

$f_infinity=11+ⅇV+35.063.6+0.61+ⅇ50-V20tau_f=10.0197⁢ⅇ-0.0337⁢V+14.52+0.02ddtimef=f_infinity-ftau_f$

### Component: ICaL_fCa_gate

$fCa_SL=1-fCaB_SLfCa_jct=1-fCaB_jctddtimefCaB_SL=1.7⁢Ca_SL⁢1-fCaB_SL-11.9⁢fCaB_SLddtimefCaB_jct=1.7⁢Ca_jct⁢1-fCaB_jct-11.9⁢fCaB_jct$

### Component: INaCa

$temp_jct=ⅇeta⁢V⁢FR⁢T⁢Na_jctHNa⁢Cao-ⅇeta-1⁢V⁢FR⁢T⁢NaoHNa⁢Ca_jct1+ksat⁢ⅇeta-1⁢V⁢FR⁢Ttemp_SL=ⅇeta⁢V⁢FR⁢T⁢Na_SLHNa⁢Cao-ⅇeta-1⁢V⁢FR⁢T⁢NaoHNa⁢Ca_SL1+ksat⁢ⅇeta-1⁢V⁢FR⁢TQ_NCX=Q10_NCXT-31010Ka_SL=11+Kd_actCa_SL3Ka_jct=11+Kd_actCa_jct3i_NaCa_jct=Fx_NCX_jct⁢V_max⁢Ka_jct⁢Q_NCX⁢temp_jctK_mCai⁢NaoHNa⁢1+Na_jctK_mNaiHNa+K_mNaoHNa⁢Ca_jct⁢1+Ca_jctK_mCai+K_mCao⁢Na_jctHNa+Na_jctHNa⁢Cao+NaoHNa⁢Ca_jcti_NaCa_SL=Fx_NCX_SL⁢V_max⁢Ka_SL⁢Q_NCX⁢temp_SLK_mCai⁢NaoHNa⁢1+Na_SLK_mNaiHNa+K_mNaoHNa⁢Ca_SL⁢1+Ca_SLK_mCai+K_mCao⁢Na_SLHNa+Na_SLHNa⁢Cao+NaoHNa⁢Ca_SLi_NaCa=i_NaCa_jct+i_NaCa_SL$

### Component: ICap

$V_maxAF_old=V_max⁢Vol_cytosol⁢2⁢FCmQ_SLCaP=Q10_SLCaPT-31010i_Cap_jct=Q_SLCaP⁢V_maxAF⁢Fx_SLCaP_jct1+KmCa_jctHi_Cap_SL=Q_SLCaP⁢V_maxAF⁢Fx_SLCaP_SL1+KmCa_SLHi_Cap=i_Cap_jct+i_Cap_SL$

### Component: ICab

$i_Cab_jct=G_CaBk⁢Fx_CaBk_jct⁢V-E_Ca_jcti_Cab_SL=G_CaBk⁢Fx_CaBk_SL⁢V-E_Ca_SLi_Cab=i_Cab_SL+i_Cab_jct$

### Component: Jrel_SR

$kCaSR=Max_SR-Max_SR-Min_SR1+EC50_SRCa_SRHSRkoSRCa=koCakCaSRkiSRCa=kiCa⁢kCaSRRI=1-R-O-IddtimeR=kim⁢RI-kiSRCa⁢Ca_jct⁢R-koSRCa⁢Ca_jct2⁢R-kom⁢OddtimeO=koSRCa⁢Ca_jct2⁢R-kom⁢O-kiSRCa⁢Ca_jct⁢O-kim⁢IddtimeI=kiSRCa⁢Ca_jct⁢O-kim⁢I-kom⁢I-koSRCa⁢Ca_jct2⁢RIj_rel_SR=ks⁢O⁢Ca_SR-Ca_jct$

### Component: Jleak_SR

$j_leak_SR=KSRleak⁢Ca_SR-Ca_jct$

### Component: Jpump_SR

$Q_SRCaP=Q10_SRCaPT-31010j_pump_SR=Q_SRCaP⁢V_max⁢Vol_cytosolVol_SR⁢CaiKmfH-Ca_SRKmrH1+CaiKmfH+Ca_SRKmrH$

### Component: ion_diffusion

$j_Na_jct_SL=D_Na_jct_SL⁢A_jct_SLx_jct_SLj_Na_SL_cytosol=D_Na_SL_cytosol⁢A_SL_cytosolx_SL_cytosolj_Ca_jct_SL=D_Ca_jct_SL⁢A_jct_SLx_jct_SLj_Ca_SL_cytosol=D_Ca_SL_cytosol⁢A_SL_cytosolx_SL_cytosolj_Na_jct_SL_new=1.8313e-14j_Na_SL_cytosol_new=1.6386e-12j_Ca_jct_SL_new=8.2413e-13j_Ca_SL_cytosol_new=3.7243e-12J_Na_jct_SL=Na_jct-Na_SL⁢j_Na_jct_SLJ_Na_SL_cytosol=Na_SL-Nai⁢j_Na_SL_cytosolJ_Ca_jct_SL=Ca_jct-Ca_SL⁢j_Ca_jct_SLJ_Ca_SL_cytosol=Ca_SL-Cai⁢j_Ca_SL_cytosol$

### Component: Na_buffer

$dNa_jct_buf=kon⁢Na_jct⁢Bmax_jct-Na_jct_buf-koff⁢Na_jct_bufdNa_SL_buf=kon⁢Na_SL⁢Bmax_SL-Na_SL_buf-koff⁢Na_SL_bufddtimeNa_jct_buf=dNa_jct_bufddtimeNa_SL_buf=dNa_SL_bufddtimeNa_jct=-Cm⁢i_Na_jct+3⁢i_NaCa_jct+i_Nab_jct+3⁢i_NaK_jct+i_CaL_Na_jctVol_jct⁢F-J_Na_jct_SLVol_jct-dNa_jct_bufddtimeNa_SL=-Cm⁢i_Na_SL+3⁢i_NaCa_SL+i_Nab_SL+3⁢i_NaK_SL+i_CaL_Na_SLVol_SL⁢F+J_Na_jct_SL-J_Na_SL_cytosolVol_SL-dNa_SL_bufddtimeNai=J_Na_SL_cytosolVol_cytosol$

### Component: Ca_buffer

$dCalsequestrin=kon_Calsequestrin⁢Ca_SR⁢Bmax_Calsequestrin⁢Vol_cytosolVol_SR-Ca_Calsequestrin-koff_Calsequestrin⁢Ca_CalsequestrinddtimeCa_Calsequestrin=dCalsequestrindCa_SLB_SL=kon_SL⁢Ca_SL⁢Bmax_SLB_SL⁢Vol_cytosolVol_SL-Ca_SLB_SL-koff_SLB⁢Ca_SLB_SLdCa_SLB_jct=kon_SL⁢Ca_jct⁢Bmax_SLB_jct⁢Vol_cytosolVol_jct-Ca_SLB_jct-koff_SLB⁢Ca_SLB_jctdCa_SLHigh_SL=kon_SL⁢Ca_SL⁢Bmax_SLHigh_SL⁢Vol_cytosolVol_SL-Ca_SLHigh_SL-koff_SLHigh⁢Ca_SLHigh_SLdCa_SLHigh_jct=kon_SL⁢Ca_jct⁢Bmax_SLHigh_jct⁢Vol_cytosolVol_jct-Ca_SLHigh_jct-koff_SLHigh⁢Ca_SLHigh_jctddtimeCa_SLB_SL=dCa_SLB_SLddtimeCa_SLB_jct=dCa_SLB_jctddtimeCa_SLHigh_SL=dCa_SLHigh_SLddtimeCa_SLHigh_jct=dCa_SLHigh_jctdCa_jct_tot_bound=dCa_SLB_jct+dCa_SLHigh_jctdCa_SL_tot_bound=dCa_SLB_SL+dCa_SLHigh_SLi_Ca_jct_tot=i_CaL_Ca_jct-2⁢i_NaCa_jct+i_Cab_jct+i_Cap_jcti_Ca_SL_tot=i_CaL_Ca_SL-2⁢i_NaCa_SL+i_Cab_SL+i_Cap_SLddtimeCa_SR=j_pump_SR-j_leak_SR⁢Vol_cytosolVol_SR+j_rel_SR-dCalsequestrinddtimeCa_jct=-i_Ca_jct_tot⁢CmVol_jct⁢2⁢F-J_Ca_jct_SLVol_jct+j_rel_SR⁢Vol_SRVol_jct+j_leak_SR⁢Vol_cytosolVol_jct-1⁢dCa_jct_tot_boundddtimeCa_SL=-i_Ca_SL_tot⁢CmVol_SL⁢2⁢F+J_Ca_jct_SL-J_Ca_SL_cytosolVol_SL-1⁢dCa_SL_tot_boundddtimeCai=-j_pump_SR⁢Vol_SRVol_cytosol+J_Ca_SL_cytosolVol_cytosol-1⁢dCa_cytosol_tot_bound$

### Component: cytosolic_Ca_buffer

$dCa_TroponinC=kon_TroponinC⁢Cai⁢Bmax_TroponinC-Ca_TroponinC-koff_TroponinC⁢Ca_TroponinCdCa_TroponinC_Ca_Mg=kon_TroponinC_Ca_Mg_Ca⁢Cai⁢Bmax_TroponinC_Ca_Mg_Ca-Ca_TroponinC_Ca_Mg+Mg_TroponinC_Ca_Mg-koff_TroponinC_Ca_Mg_Ca⁢Ca_TroponinC_Ca_MgdMg_TroponinC_Ca_Mg=kon_TroponinC_Ca_Mg_Mg⁢Mgi⁢Bmax_TroponinC_Ca_Mg_Mg-Ca_TroponinC_Ca_Mg+Mg_TroponinC_Ca_Mg-koff_TroponinC_Ca_Mg_Mg⁢Mg_TroponinC_Ca_MgdCa_Calmodulin=kon_Calmodulin⁢Cai⁢Bmax_Calmodulin-Ca_Calmodulin-koff_Calmodulin⁢Ca_CalmodulindCa_Myosin=kon_Myosin_Ca⁢Cai⁢Bmax_Myosin_Ca-Ca_Myosin+Mg_Myosin-koff_Myosin_Ca⁢Ca_MyosindMg_Myosin=kon_Myosin_Mg⁢Mgi⁢Bmax_Myosin_Mg-Ca_Myosin+Mg_Myosin-koff_Myosin_Mg⁢Mg_MyosindCa_SRB=kon_SRB⁢Cai⁢Bmax_SRB-Ca_SRB-koff_SRB⁢Ca_SRBdCa_cytosol_tot_bound=dCa_TroponinC+dCa_TroponinC_Ca_Mg+dMg_TroponinC_Ca_Mg+dCa_Calmodulin+dCa_Myosin+dMg_Myosin+dCa_SRBddtimeCa_TroponinC=dCa_TroponinCddtimeCa_TroponinC_Ca_Mg=dCa_TroponinC_Ca_MgddtimeMg_TroponinC_Ca_Mg=dMg_TroponinC_Ca_MgddtimeCa_Calmodulin=dCa_CalmodulinddtimeCa_Myosin=dCa_MyosinddtimeMg_Myosin=dMg_MyosinddtimeCa_SRB=dCa_SRB$

### Component: indo_fluo_Ca_buffer_not_connected

$Indo1Bound=Ca_Indo1_Cytosol+Ca_Indo1_jct⁢Vol_jctVol_cytosol+Ca_Indo1_SL⁢Vol_SLVol_cytosolFluo3Bound=Ca_Fluo3_Cytosol+Ca_Fluo3_jct⁢Vol_jctVol_cytosol+Ca_Fluo3_SL⁢Vol_SLVol_cytosoldCa_Indo1_Cytosol=kon_Indo1⁢Indo1-Indo1Bound⁢Cai⁢Bmax_Indo1_Cytosol-Ca_Indo1_Cytosol-koff_Indo1⁢Ca_Indo1_CytosoldCa_Indo1_jct=kon_Indo1⁢Indo1-Indo1Bound⁢Ca_jct⁢Bmax_Indo1_jct⁢Vol_cytosolVol_jct-Ca_Indo1_jct-koff_Indo1⁢Ca_Indo1_jctdCa_Indo1_SL=kon_Indo1⁢Indo1-Indo1Bound⁢Ca_SL⁢Bmax_Indo1_SL⁢Vol_cytosolVol_SL-Ca_Indo1_SL-koff_Indo1⁢Ca_Indo1_SLdCa_Fluo3_Cytosol=kon_Fluo3⁢Fluo3-Fluo3Bound⁢Cai⁢Bmax_Fluo3_Cytosol-Ca_Fluo3_Cytosol-koff_Fluo3⁢Ca_Fluo3_CytosoldCa_Fluo3_jct=kon_Fluo3⁢Fluo3-Fluo3Bound⁢Ca_jct⁢Bmax_Fluo3_jct⁢Vol_cytosolVol_jct-Ca_Fluo3_jct-koff_Fluo3⁢Ca_Fluo3_jctdCa_Fluo3_SL=kon_Fluo3⁢Fluo3-Fluo3Bound⁢Ca_SL⁢Bmax_Fluo3_SL⁢Vol_cytosolVol_SL-Ca_Fluo3_SL-koff_Fluo3⁢Ca_Fluo3_SLddtimeCa_Indo1_Cytosol=dCa_Indo1_CytosolddtimeCa_Indo1_SL=dCa_Indo1_SLddtimeCa_Indo1_jct=dCa_Indo1_jctddtimeCa_Fluo3_Cytosol=dCa_Fluo3_CytosolddtimeCa_Fluo3_SL=dCa_Fluo3_SLddtimeCa_Fluo3_jct=dCa_Fluo3_jct$

### Component: intracellular_currents

Source
Derived from workspace Shannon, Wang, Puglisi, Weber, Bers, 2004 at changeset bb85ada8052b.
Collaboration
To begin collaborating on this work, please use your git client and issue this command: