Location: Tran 2017 - Cross-bridge model of shortening heat @ eb49874a619c / Fevents.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://staging.physiomeproject.org/workspace/4a2/rawfile/eb49874a619c53bd0b587ea6854313751e76edd3/Fevents.m
Compare... Download Source 


function [value,isterminal,direction] = Fevents(t,y,Phase,Params)
% Locate the time when the active force equals the afterload at 
% which point, the integration should stop.

    [dy F_total ATPase F_net] = XBModel(t,y,Phase,Params);

    value = F_net;     % Detect height = 0
    isterminal = 1;   % Stop the integration
    direction = 0;   % Negative direction only
end