Model Mathematics

Component: environment

Component: intracellular_ion_concentrations

mtime=time-167time1671 Ca_i=11.8535847759 e -5mtime3+-0.001590341066858mtime2+0.043645898640826mtime1+0.167079436392938ifmtime1.17mtime<30.841-5.7458491 e -8mtime3+3.1122215647 e -5mtime2-0.006618494920814mtime1+0.720441888689583ifmtime30.840.216otherwise

Component: crossbridges

permtot=TRPNperm50nperm ddtimexb=k_xbpermtot1-xb-1permtotxb

Component: troponin

Ca_50=Ca_50ref1+beta_1lambda_m-1 ddtimeTRPN=k_offCa_iCa_50n_TRPN1-TRPN-TRPN

Component: Myofilaments

lambda=1 dlambdadt=0

Component: filament_overlap

lambda_m=1.2iflambda>1.2lambdaotherwise lambda_s=0.87iflambda_m0.87lambda_motherwise overlap=1+beta_0lambda_m+lambda_s-1.87

Component: isometric_tension

T_0=T_refxboverlap

Component: dynamic_stiffness

Q=Q_1+Q_2 Tension=T_0aQ+11-QifQ<0T_01+a+2Q1+Qotherwise ddtimeQ_1=A_1dlambdadt-alpha_1Q_1 ddtimeQ_2=A_2dlambdadt-alpha_2Q_2