Hello,

Quelqu'un pourrait-il me renseigner sur la fonction ode45 ?
La réponse que me donne ode45 pour le système d'équations différentielles suivant est différente selon l'ordre des équations dy(...) = ... :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
function dy = dABPQ(t,y,p)
dy = zeros(5,1);
 
if t>p.tdos;
   F=0;
else
   F=p.F;
end
 
dy(1) = F; 
dy(2) = -p.k1*y(2)*y(3) -p.k2*y(2)*y(3) -y(2)/y(1)*dy(1)                  ;
dy(3) = -p.k1*y(2)*y(3) -p.k2*y(2)*y(3) -y(3)/y(1)*dy(1) +p.Bini/y(1)*F   ;
dy(4) =  p.k1*y(2)*y(3)                 -y(4)/y(1)*dy(1)                  ;
dy(5) =                  p.k2*y(2)*y(3) -y(5)/y(1)*dy(1)                  ;
 
 
end
la réponse si dy(1) = F est déplacé à la fin, et donc devient dy(5), est différente

Pour le contexte il s'agit de deux réactions parallèles et dy(1), y(1) sont respectivement la variation de volume et le volume.