Méthode de Runge Kutta sous scilab

[invite55056469]

Voilà je viens de débuter les équations diff sous scilab et je commence déjà par zigzager. Auriez vous l'amabilité de jeter un coup d'oeil à ce programme et me donner un coup de pousse les topos?
Merci d'avance...

function yprime=f(x,y);
yprime=3 .*cos(2*x) - y;
endfunction;
plot2d(x,y);
function [x,y]=RK1(f,x0,y0,h,n);
h= 1 ./n;
x(1)=0 //x0 + h;
y(1)=1//y0 + (h .*f(x0,y0));
for i=1:n
y(i+1)=y(i-1)+ (h .*f(x(i-1),y(i-1)));
x(i+1)=x0+ (i+1) .*h;
end
endfunction

n=10;
h=0.1;
x0=0;
y0=1;
[xm,ym] = RK1(f,x0,y0,h,n)
-----

[invite55056469]

Petite correction les potos.

clear
// méthode d’Euler ou Runge Kutta 1
// f: la fonction, (x0,y0): le point initial
// h: le pas et n: le nombre de points
// le résultat est un tableau contenant deux listes x et y
function yprime=f(x,y);
yprime=3 .*cos(2*x) - y;
endfunction;

function [x,y]=RK1(f,x0,y0,h,n);
h= 1 ./n;
x(1)=0 //x0 + h;
y(1)=1//y0 + (h .*f(x0,y0));
for i=1:n
y(i+1)=y(i-1)+ (h .*f(x(i-1),y(i-1)));
x(i+1)=x0+ (i+1) .*h;
end
endfunction

n=10;
h=0.1;
x0=0;
y0=1;
[x,y] = RK1(f,x0,y0,h,n)

Nota:svp faites vite. ça urge.