Location: Tran 2017 - Cross-bridge model of shortening heat @ eb49874a619c / IC_XBModel.m

Author:

Kenneth Tran <k.tran@auckland.ac.nz>

Date:

2017-08-31 14:58:36+12:00

Desc:

Adding MatLab code for users to reproduce simulation figures in the J Physiol publication: DOI: 10.1113/JP274680.

Permanent Source URI:

https://models.cellml.org/workspace/4a2/rawfile/eb49874a619c53bd0b587ea6854313751e76edd3/IC_XBModel.m

%% Initial State variable values

Nxb = 0.9997;  % Non-permissive fraction in overlap region
XBpreR = 0.0001;   % Fraction of XBs in pre-rotated state (Stiffness generation)
XBpostR = 0.0001;  % Fraction of XBS in post-rotated state (force generation)
x_XBpreR = 0;    % Strain of XBs in the pre-rotated state
x_XBpostR = 0.007;  % Strain of XBs in the post-rotated state
TropCaL = 0.0144; % Ca bound to the low affinity troponin site
TropCaH = 0.1276; % Ca bound to the high affinity troponin sit
IntegF = 0;
SL = 2.3; % This is changed in SBSolve.m

init = [Nxb XBpreR XBpostR x_XBpreR x_XBpostR TropCaL TropCaH IntegF SL];