Bonjour à tous
Je suis débutant sur Scilab et je dois résoudre par la méthode d’Euler un système proie-prédateur vectoriellement.
Je dois tracer le portrait de phase en utilisant la méthode Euler explicite.
J'y arrive en utilisant la méthode ode mais pas en écrivant moi-même la methdoe d'euler.

Voici ce que j'ai déjà écrit :

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
18
functionXp=F(X)
    a=log(2)/90
    b=log(2)/900
    c=log(10)/30
    d=log(10)/3000
    Xp(1)=a*X(1)-b*X(1)*X(2)
    Xp(2)=-c*X(2)+d*X(1)*X(2)
endfunction
 
function [x,y]=orbite(N,x0,y0)
    dt=1/N
    t=linspace(0,1,N+1)
    for i=1:N
        x(i)=x(i-1)+dt*x(i-1)
        y(i)=y(i-1)+dt*y(i-1)
    [X]=F(X)
    end
endfunction
Je vous remercie