Question sur la fonction ode pour le système de Lorentz et proie prédateur

Question sur la fonction ode pour le système de Lorentz et proie prédateur

par Stourbe Alban,
Nombre de réponses : 3

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 ! :)
En réponse à Stourbe Alban

Re: Question sur la fonction ode pour le système de Lorentz et proie prédateur

par Mottelet Stephane,

Bonjour,

Je pense qu'il manque un "*" ici : SG(y-x)

Donc plutot ecrire SG*(y-x).

Cordialement

En réponse à Mottelet Stephane

Re: Question sur la fonction ode pour le système de Lorentz et proie prédateur

par Mottelet Stephane,
Et il manque aussi toutes les autres multiplications !
px -> p*x, etc.
S.