calcPassiveForceLengthCurve.m

% This function is part of the NMSM Pipeline, see file for full license.
%
% This function calculates the passive force of a muscle fiber based on
% the normalized fiber length and the maximum isometric force.
%
% Inputs:
%   normalizedFiberLength: 3d mat (1, numMuscles, 1) of normalized fiber
%   maxIsometricForce: 3d mat (1, numMuscles, 1) of maximum isometric force
%   pennationAngle: 3d mat (1, numMuscles, 1) of pennation angle
%
% Outputs:
%   passiveForce: 3d mat (1, numMuscles, 1) of passive passiveForce
%
% (3d mat, 3d mat, 3d mat) -> 3d mat
% returns normalized passive force


function passiveForce = calcPassiveForceLengthCurve( ...
    normalizedFiberLength, maxIsometricForce, pennationAngle)
passiveForce = maxIsometricForce  .* cos(pennationAngle) .* ...
    passiveForceLengthCurve(normalizedFiberLength);
end