Bonjour Monsieur,
Depuis tout à l'heure j'essaye de faire les exercices sur le système de Lorentz et proie prédateur. Je pense avoir trouver le bon X' = f(X) et les points d'équilibre pour les deux systèmes. Pourtant sur Scilab je n'arrive pas à modéliser ces systèmes, je n'arrive pas à calculer X avec la fonction ode(). Dans les deux cas, Scilab m'envoie cette erreur :
Indice invalide.
ode: An error occurred in 'lsoda' subroutine.
J'ai du avoir mal choisi mes valeurs pour les conditions initiales ou les valeurs de t.
Par exemple pour le problème de Lorentz j'ai écrit :
function dXdt=f(t, X) SG = 10; p = 28; b = 8/3; x = X(1); y = X(2); z = X(3); dXdt = [SG(y-x) px - y - xz xy -bz]; endfunction t = linspace(0,10,100) X = ode([0;0;0],0,t,f);
Pouvez vous m'aider s'il vous plait ! :)