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

Author:
Kenneth Tran <k.tran@auckland.ac.nz>
Date:
2017-08-31 15:20:50+12:00
Desc:
Adding image.
Permanent Source URI:
https://models.cellml.org/workspace/4a2/rawfile/b09963ece3950a7e2c1256070cab0b5e13411d90/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