Bonjour,
J'ai écrit un programme sous Scilab afin de tracer la solution approchée d'une suite d'Euler vn d'une équation différentielle.
Mon problème n'est pas tant théorique mais plutôt lié à la pratique de Scilab !
Voici mon programme :
clf()
clear()
N=50;
h=40/N;
n=1:h:41
m=80;
g=9.81;
k1=g*80/(200000/3600);
v(1)=0;
for i=1:(length(n)-1)
v(i+1)=h*g+v(i)*(1-(h*k1)/m)
end
plot2d(n', v(n+1))
Dans le programme ci-dessus, j'obtient bien une courbe lorsque N = 50, mais dès que je fait varier N en dessous de 40 j'ai l'erreur "Indice invalide." sur mon instruction plot2d.
Je m'arrache les cheveux dessus car je dois tracer la même courbe pour N=10 !
Merci d'avance pour vos conseils et votre aide
-----