Apprendre Maple

[laure12345]

Bonjour! J'ai un petit probleme avec mon programme maple. Je m'explique: J'ai résolu une equation non linéaire avec le type numeric bien sur, sauf que maintenant je veux integrer mes solutions, j'ai donc créer une procedure avec la méthode des trapèzes (on m'a aidé) et maintenant j'aimerai en tracer la courbe (en dimension 3d), et je pense qu'il faut le faire avec la fonction pointplot3d, mais je ne sais pas l'utiliser avec une procédure, seulement avec des points et leurs coordonnées (et l'aide ne m'aide pas plus que ça...). Merci d'avance

[zozo]

Bonjour Laure,
Si tu crées dans une procédure P ta liste de points L, tu insères pointplot3d(L) en fin de procédure, et donc tu auras le graphique comme résultat de la procédure P.
Donc le programme a la structure suivante, à compléter par tes données:

with(plots);

P:=proc(....)

L:=......
pointplot3d(L);

end proc;

@+

[laure12345]

Merci beaucoup!! Je suis confrontée a un nouveau probleme, rien a voir avec ma procédure. J'ai résolu mon systeme d'equa diff avec l'aide de la fonction numéro, et maintenant j'aimerai déterminer précisement t, pour lequel x(t)=0, j'écris donc betement s:=solve(x(t)=0,t), et forcement ça ne marche pas, avez une idée ?

[laure12345]

*numérique.

[zozo]

Bonjour,
Tout dépend de l'équation concernée: si elle est vraiment compliquée, Maple ne pourra pas la résoudre avec solve et soit affichera aucun résultat ou parfois RootOf(x(t),index=1),....
Dans ce cas, il faut essayer une résolution approchée du style

s:=fsolve(x(t),t);

ou

s:=fsolve(x(t),t=a..b); si l'on sait que la solution est entre a et b.

A plus tard.

[laure12345]

Je n'y arrive pas, si j'écris seulement s:=fsolve(x(t),t=a..b) maple ne peut pas comprendre que je veux déterminer t pour x(t)=0 non?
J'ai donc écris:
G:=dsolve({eqx,eqy,eqz,ei},{x(t),y(t),z(t)},numeric,output=listprocedure);
(puis je trace ma courbe) (donc ça OK)
puis
S:=fsolve(z(t)=0,t=0.5..0.6);
Mais maple me recopie ce que je viens d'écrire... (t=0.55 environ)

[zozo]

Essayez

G:=dsolve({eqx,eqy,eqz,ei},{x(t),y(t),z(t)},numeric,output=listprocedure);

puis

assign(%);

Ceci permettra de mettre dans la variablre z(t) la solution trouvée.

Enfin faire

S:=fsolve(z(t)=0,t=0.5..0.6);

ZOZO