Location: BG_LCC_diadic @ 3b58efb1f452 / matlab_parameter_fitting / dx_f.m

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2022-06-03 11:31:35+12:00
Desc:
Change name of V_di
Permanent Source URI:
https://staging.physiomeproject.org/workspace/891/rawfile/3b58efb1f452d11b3fb8d105e9b1cb48b4e39cda/matlab_parameter_fitting/dx_f.m

function dx = dx_f(x,V)

alpha1_0 = 5;
beta1_0 = exp(3);

z_f1 = 10;
z_r1 = 100;

alpha2_0 = 5;
beta2_0 = exp(3);

z_f2 = -10;
z_r2 = -100;

alpha1 = alpha1_0*exp(z_f1*V);
beta1 = beta1_0*exp(z_r1*V);

alpha2 = alpha2_0*exp(z_f2*V);
beta2 = beta2_0*exp(z_r2*V);

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

dx = [v1; -(v1+v2); v2];
end