Location: BG_LCC_diadic @ 811cbd9176a1 / matlab_parameter_fitting / dx_f.m

Author:

Shelley Fong <sfon036@UoA.auckland.ac.nz>

Date:

2022-07-07 11:10:27+12:00

Desc:

Update readme

Permanent Source URI:

https://staging.physiomeproject.org/workspace/891/rawfile/811cbd9176a1dcddd93f0a389a5d99a75ce0523d/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