équation différentielle scilab

Fermé
mansy1010 Messages postés 1 Date d'inscription dimanche 5 décembre 2010 Statut Membre Dernière intervention 5 décembre 2010 - Modifié par mansy1010 le 5/12/2010 à 17:01
Bonjour,

Je suis une débutante en scilab et j'ai un soucis avec la méthode de numerov, j'ai réussi à écrire le programme mais j'ai un bug! Pouvez vous m'aider svp car j'ai un contrôle dessus cette semaine :(
Voici ce que j'ai fait

function z=f(x,y)
z=4*y
endfunction

format("v",15);
n=10.0;
x0=0;
y0=0.2;
x1=0.1;
y1=0.5;
yb0=y0;
h=x1-x0;

x=linspace(x0,x0+n*h,n+1);
yc=linspace(y0,y0+n*h,n+1);
ye=-0.63505643.*exp(-2*x1)+0.83505643.*exp(2*x1);
---> plot2d(x,ye,1); LE SOUCIS EST LA, on me dit que les 2 sont de dimensions incompatibles!!!
for i=3:n+1
x2=x1+h;
y2p=2*y1-y0+(h^2)/11*(10*f(x1,y1)+f(x0,y0));
y2=2*y1-y0+(h^2)/12*(f(x2,y2p)+10*f(x1,y1)+f(x0,y0));
yc(i)=y2;
x1=x2;
y1=y2;
end;
plot2d(x,yc,5);

C'est peut être bête comme truc mais je n'y connais rien donc ça n'aide pas!

Merci d'avance