MAPLE et methode d'euler

[invite42abb461]

Bonjour, je souhaite utiliser MAPLE pour tracer une solution de l'equadiff suivante :

[1+y'²]y=k²

J'ai donc choisi k =1 pour simplifier et j'ai ecrit la procédure suivante :

> restart:with(plots):

> f:=x-> ( 1 - 1/x)

> x[0] := 1; y[0] := 5;

> for k to 10 do y[k + 1] := y[k] + 0.1 f(x[k]) od;

> l := [[x[n], y[n]] $ (n = 1 .. 10)];

> plot(l, style = line)




Alors, MAPLE me répond :
"Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

Error, empty plot"

Le probleme semble venir de ma fonction f mais je ne sais comment le corriger. Certains d'entre vous peuvent ils me venir en aide ?
Merci d'avance.
-----

[cedbont]

Bonjour,
j'ai exécuté ça (ça ressemble un peu à ce que tu postes) :

> restart:with(plots):

> f:=x-> 1 - 1/x:

> x[0] := 1;

> y[0] := 5;

> for k from 0 to 10 do y[k + 1] := y[k] + 0.1*f(x[k]) od;

> l := [[x[n], y[n]] $ (n = 1 .. 10)];

> plot(l, style = line);

J'ai obtenu ce que j'ai mis en pièce jointe. Manifestement ton problème vient du fait que tu ne définis pas les x[k].

[invite72eff345]

faire du MAPLE le samedi, vous êtes tarés les gars!

[invite42abb461]

Ok merci mais sur ton fichier, le plot ne marche pas...

[A voir en vidéo sur Futura]

[cedbont]

Exact, puisque tu lui demandes de t'afficher les points de coordonnées (x[k],y[k]). Or, il ne connaît pas les x[k].

[invite42abb461]

Comment corriger ca au niveau de la syntaxe ? Désolé mais je suis une quiche en MAPLE.

[cedbont]

Dans ton for tu dois rajouter une ligne pour dire comment tu incrémentes tes x[k], par exemple;
x[k+1] := x[k] + dx
Tu devras bien sûr préciser la valeur du dx, puisque ce que la commande plot attend est un couple de deux nombres.

[invite42abb461]

j'ai essayé avec dx = 0.01 il me dit toujours la meme chose. Cependant il comprend la boucle now, mais ne trace toujours rien.

[cedbont]

Peux-tu mettre ce que tu as écrit et ce que Maple te répond ?

[invite42abb461]

En fait j'ai trouvé, cela venait du fait que je partais de 1 et pas de zero dans la boucle
Sinon qqun sait faire un programme pour tracer la meme chose avec Runge Kutta ?