Model Mathematics

Component: Nai

ddtimeNai=Nai_rate

Component: Vstim_para

Component: control_para

Component: outputs

Component: Unit_uSMC

Component: K_1

K_1=Cain_MCa_halfMLCKn_M+Cain_M⁢1

Component: E_Na

Component: E_Ca

Component: E_K

Component: Nernst_potential

E=R⁢Tz⁢F⁢ln⁡X_oX_i

Component: Nernst_potential

E=R⁢Tz⁢F⁢ln⁡X_oX_i

Component: Nernst_potential

E=R⁢Tz⁢F⁢ln⁡X_oX_i

Component: CB4HM

K_4=K_34 K_5=K_2 K_6=K_1 norm=M+Mp+AM+AMp ddtimeM=-K_1⁢Mnorm+K_2⁢Mpnorm+K_7⁢AMnorm ddtimeMp=K_4⁢AMpnorm+K_1⁢Mnorm-K_2+K_3⁢Mpnorm ddtimeAM=K_5⁢AMpnorm-K_6+K_7⁢AMnorm ddtimeAMp=K_3⁢Mpnorm+K_6⁢AMnorm-K_4+K_5⁢AMpnorm stress=AMp+AM phosphorylation=AMp+Mp

Component: J_VOCC

rho_vCa=11+ⅇV_Cahalf-VK_Cahalf J_VOCC=-I_VOCC2⁢V_cell⁢F

Component: J_CaPump

J_CaPump=-V_pmax⁢CainK_phn+Cain

Component: J_NaCa

V_mNaCa=3⁢E_Na-2⁢E_Ca J_NaCa=G_NaCa⁢CaiCai+K_NaCa⁢V-V_mNaCa

Component: Cai

ddtimeCai=J_VOCC+inhPump⁢J_CaPump+J_NaCa

Component: Ionic_currents

I=g_max⁢PO⁢V-E

Component: sPulse_protocol_ms

V=V_actHoldingiftime<t_ssV_actTestiftime≥t_ss∧time<t_act+t_ssV_actHoldingotherwise

Component: sPulse_protocol_s

V=V_actHoldingiftime<t_ssV_actTestiftime≥t_ss∧time<t_act+t_ssV_actHoldingotherwise

Component: mPulse_protocol_ms

stimPeriod=t_act+t_intp V=V_actHoldingiftime<t_ssV_actTestiftime≥t_ss∧time-t_ssmodstimPeriod<t_act∧time≦t_ss+stimPeriod⁢npV_actHoldingotherwise

Component: mPulse_protocol_s

stimPeriod=t_act+t_intp t=time-t_ss-stimPeriod⁢⌊time-t_ssstimPeriod⌋ V=V_actHoldingiftime≦t_ssV_actTestiftime>t_ss∧t≦t_act∧time≦t_ss+stimPeriod⁢npV_actHoldingotherwise

Component: time_ms

Component: time_s

Component: constants

Component: model_para

Component: initials