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

% Function to generate Cai transient

function Cai = getCai(t)

    tau1 = 20; %(ms)
    tau2 = 33; %(ms)
    Ca_amp = 1.5; %(uM)
    Ca_dias = 0.09; %(uM)
    t_start = 50; %(ms)
    
    beta = (tau1/tau2)^(-1/(tau1/tau2 -1)) - (tau1/tau2)^(-1/(1 - tau2/tau1));
    
    if (t>t_start)
        Cai = ((Ca_amp - Ca_dias)/beta)*(exp(-(t-t_start)/tau1) - exp(-(t-t_start)/tau2)) + Ca_dias;
    else
        Cai = Ca_dias;
    end
    
end