Location: BG_LCC @ 7f2940423be8 / matlab_parameter_fitting / df_BG.m

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2022-07-06 16:27:15+12:00
Desc:
Updating to LRd 34.4 pL cell parameters
Permanent Source URI:
https://models.physiomeproject.org/workspace/6d7/rawfile/7f2940423be8efc20786e70c4edc7ead28f346f5/matlab_parameter_fitting/df_BG.m
Compare... Download Source 

function dxdt = df_BG(x,t,t_vec,V_vec,params)
R = 8.314;
T = 310;
F = 96485;

V = interp1(t_vec,V_vec,t);

alpha1 = params(1)*exp(params(2)*F*V/R/T);
beta1 = params(3)*exp(params(4)*F*V/R/T);
alpha2 = params(5)*exp(params(6)*F*V/R/T);
beta2 = params(7)*exp(params(8)*F*V/R/T);
k3f = beta1*alpha2/alpha1/beta2 * 1e5;
k3r = 1e5;

v1 = alpha1*x(2) - beta1*x(1);
v2 = alpha2*x(2) - beta2*x(3);
v3 = k3f*x(1) - k3r*x(3);

dxdt = [v1-v3; -v1-v2; v2+v3];

end