Commit 46a5ca2e by Chiara Antonini

Upload new file

1 parent 30a1d7b9
function dx=modelPeng2016(t,x,p,u,xT)
lambda_A=p(1);
r_p1=p(2);
r_a1=p(3);
r_m=p(4);
k_CX=p(5);
r_p2=p(6);
r_a2=p(7);
alpha_XD=p(8);
pigreco_D=p(9);
alpha_DC=p(10);
k_RC=p(11);
mi_C=p(12);
alpha_DR=p(13);
pigreco_R=p(14);
alpha_IR=p(15);
alpha_XR=p(16);
mi_R=p(17);
alpha_CI=p(18);
mi_I=p(19);
alpha_DCDR=p(20);
mi_D=p(21);
pigreco_C=p(22);
alpha_DRR=p(23);
A=x(1);
X1=x(2);
X2=x(3);
D_m=x(4);
C2=x(5);
R2=x(6);
I2=x(7);
D_C=x(8);
D_R=x(9);
C1=x(10);
R1=x(11);
I1=x(12);
dA = lambda_A*(1-A);
dX1 = r_p1 * A * X1 - r_a1*(1-A)*X1 - r_m*(1-A)*X1 - k_CX*C2*X1;
dX2 = r_p2 * X2 - r_a2 * X2 + r_m*(1-A)*X1 - k_CX * C2 * X2;
dD_m = alpha_XD * (r_a1*(1-A)*X1 + k_CX*C2*X1 + r_a2*X2 + k_CX*C2*X2) - pigreco_D*D_m; %qui ci andrebbe un meno davanti
dC2 = alpha_DC*D_m + 0.5*pigreco_C*C1 - k_RC*R2*C2 - mi_C*C2;
dR2 = alpha_DR*D_m + 0.5 * pigreco_R*R1 + alpha_IR*I2 + alpha_XR*(X1+X2) - mi_R*R2;
dI2 = alpha_CI * C2 - mi_I*I2;
dD_C = 0.5*pigreco_D*D_m - alpha_DCDR*D_C;
dD_R = alpha_DCDR*D_C - mi_D*D_R;
dC1 = alpha_DC*D_C - mi_C*C1 - k_RC*R1*C1 -pigreco_C*C1;
dR1 = alpha_DRR*D_R + alpha_IR*I1 - mi_R*R1 - pigreco_R*R1;
dI1 = alpha_CI * C1 - mi_I*I1;
dx=[dA;dX1;dX2;dD_m;dC2;dR2;dI2;dD_C;dD_R;dC1;dR1;dI1];
end
\ No newline at end of file \ No newline at end of file
Styling with Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!