Bonjour !
Je suis étudiante en 2e année de licence de bio et je participe à un projet de recherche de microbiologie ( en gros ). Nous voulons modéliser la dynamique de 4 populations : B, I, P et R.
On a le système suivant :
dB/dt=αB(t)-βB(t)P(t)-ωB(t)P(t)
dI/dt=γI(t)- μ(auto)I(t)-h(t) μ(i)I(t)
dP/dt=μ(auto)λ(auto)I(t)+h(t) μ(i) λ(i)I(t)+ β λ(I)B(t)P(t)- ωB(t)P(t)
dR/dt=φR(t)+ ωB(t)P(t)
avec pas mal de constantes.
J'ai tenté d'écrire ce code scilab en donnant certaines valeurs à mes constantes:
function ydot=f(t,y)
ydot[1]=50*B-45*B*P-5*B*P
ydot[2]=50*I-40*I
ydot[3]=10*40*I+45*10*B*P-5*B*P
ydot[4]=50*R+5*B*P+40*R
endfunction
y0[1]=10;y0[2]=10;y0[3]=0;y0[4]=0;t0=0.1;t=0.1:0.1:%pi;
y=ode(y0,t0,t,f)
plot(t,y)
J'obtiens que t et y0 ne sont pas définis, ainsi qu'une erreur 276 sur "ydot[1]=50*B-45*B*P-5*B*P" !
Quelqu'un comprend pourquoi ça ne marche pas ?
Merci de votre aide !
-----