Location: Tran 2017 - Cross-bridge model of shortening heat @ 1ff97ce0ff62 / passiveForces.m

Author:
Kenneth Tran <k.tran@auckland.ac.nz>
Date:
2017-08-31 15:06:57+12:00
Desc:
Uploading again due to openCor bug.
Permanent Source URI:
https://models.cellml.org/workspace/4a2/rawfile/1ff97ce0ff62973dcc4434f5d0d3acfc5285a1a8/passiveForces.m
Compare... Download Source 

% 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