what's wrong my code about ode model

4 visualizaciones (últimos 30 días)
기훈 강
기훈 강 el 9 de Jun. de 2023
Comentada: Walter Roberson el 9 de Jun. de 2023
I had to training course 'solve ordinary differential equation in MATLAB online course'
and I wrote that one but this code doesn't run
how can I modify that??
The first error has line 3 ' C=Y(1)'
function dY=ESdynamics(t,Y,r_C,Cprime,r_max,k,r_A,delta_A,r_I,delta_I,r_E,Eprime,r_S,Sprime,beta,gamma)
C=Y(1);
A=Y(2);
I=Y(3);
E=Y(4);
S=Y(5);
f(C)=min(r_C*(1-C/Cprime),r_max);
dC=f(C)*C-k*C*E;
dA=r_A*C-delta_A*A;
dI=r_I*C*E-delta_I*I;
dE=-r_E*(E-Eprime)+beta*A*I*E*S-gamma*E*S;
dS=-r_S*(S-Sprime)-beta*A*I*E*S+gamma*E*S;
end
  3 comentarios
기훈 강
기훈 강 el 9 de Jun. de 2023
thanks
Walter Roberson
Walter Roberson el 9 de Jun. de 2023
Note that ode45 and related functions will produce incorrect answers when min() or max() are used.

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Ordinary Differential Equations en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by