Location: Tran 2017 - Cross-bridge model of shortening heat @ eb49874a619c / passiveForces.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/passiveForces.m

% Function to calculate the passive force
% Passive force function is fitted to passive data presented in Fig 2D

function PF = passiveForces(SL, passive)

    if passive

        x = SL/2.3;
        SHAM = [-1212.77 4497.99 -5565.60 2298.68]/110;
        PF = SHAM(1) + SHAM(2).*x + SHAM(3).*x.^2 + SHAM(4).*x.^3;

    else
        
        PF = 0;
end