Besoin d'aide pour tracé sous Maple (équa-diff "cosmologique")

[invitebd8dbca5]

Bonjour.
Je poste dans ce forum suite à un topic que j'ai créé dans le forum astrophysique. Le sujet dans lequel je voulais trouvé une équation a ensuite basculé sur un problème de tracé sous Maple et je pense que trouverai plus d'aide ici.
Le topic en question est :
http://forums.futura-sciences.com/po...51.html#828051

Le plus important des messages est sans doute le n°13 les messages suivant :

n°13

On a les différentes équations (qui sont toutes égales c'est juste une différence de notation )

------------------------------------------------

(Eq10) :

------------------------------------------------

(Eq11) :

------------------------------------------------

(Eq9) :

------------------------------------------------

(Eq12) :

------------------------------------------------

Avec

Le but est de tracer a(x) sous Maple.

Pour s'aider j'ai trouvé ce document (voir partie 5) :
http://www.isteem.univ-montp2.fr/GRA...rie/index.html

Ils obtiennent ces graphes (cliquer sur les ovales):
http://www.isteem.univ-montp2.fr/GRA...)/tableau.html

Maintenant un exemple de ce que j'obtiens en prenant :






donc



Normalement je devrais obtenir d'après le site :
http://www.isteem.univ-montp2.fr/GRA...ES/1_0.51.html

et j'obtiens :
http://img124.imageshack.us/img124/6...omaple1qa6.jpg

En plus avec certaines valeurs j'obtiens des moitiés de graphes (genre avec 0.7 et -1) et pas les graphes complets

(Q12) Je pense qu'en fait le problème vient plus de la rédaction sous maple qu'autre chose. Donc si quelqu'un avait la bonté de m'expliquer comment faire ce serait vraiment sympa . (ou mieux encore si quelqu'un avait le temps de rédiger un truc sous Maple qui marche mais bon...) Merci.

n°15

Oui je sais que ça a globalement la même tête, mais le problème c'est que je n'obtiens pas le graphe qu'il obtient (sur le site). Je m'explique : la plupart du temps, sous maple, avec les exos classiques (et donc adaptés), on obtient de belles courbes comme sur le témoin. Là le problème c'est que je n'ai pas une seule courbe (ça trace un espèce d'ensemble de courbes (même si une seule apparait en jaune)), et pour certaines valeurs, je n'obtiens pas de courbe complète.
J'ai essayer avec la valeur absolue, mais pour certaines valeurs des omegas, ça change complètement l'équation, et je n'obtiens plus de courbes ressemblant au témoin.
Pour par exemple 0,7 et -1, je devrais obtenir :
http://www.isteem.univ-montp2.fr/GRA...S/0.7_-11.html
et je n'obtiens que la première partie, c'est à dire jusqu'à ce que a(x) atteigne le maximum. Et pareil je n'ai pas une seule courbe mais un espèce d'ensemble de courbes.
Voilà le problème .

n°16 :

Bon juste une petite précision :
Dans maple, il y a deux façons de procéder pour résoudre une équa-diff.
Soit on utilise DEplot (ce que j'ai fait) qui trace directement le résultat.
Soit on utilise Dsolve qui calcule la solution de l'équa-diff, puis on trace le graphique qui correspond au résultat.
Le problème c'est que la deuxième solution qui marche en général, ne fonctionne pas chez moi : Maple moulinne mais n'a pas l'air de trouver la solution.
J'ai essayer vaguement sous Mathematica, mais apparemment ça ne fonctionne pas non plus.

Voilà merci beaucoup à ceux qui m'aideront .
-----

[invitebd8dbca5]

un petit up

[pephy]

bonjour
je peux essayer de voir ce qui cloche mais il me faudrait le fichier Maple, parce que bonjour les équations....

[pephy]

Re,
pas moyen de faire passer les fichiers Maple en pièce jointe.J'ai pu quand même récupérer l'Eq11 et faire quelques essais;chez moi çà ne marche ni avec DEplot , ni avec dsolve.
Dans ce cas je sors mon arme fatale la résolution numérique.
http://img166.imageshack.us/img166/2...ysique1nz8.jpg
Peut être que c'est bon comme çà? Je n'ai pas testé d'autres valeurs

[A voir en vidéo sur Futura]

[invitebd8dbca5]

Merci de t'être penché sur le problème. En faite, j'ai découvert et essayé la résolution numérique sous Maple hier soir. Pour les graphes où l'univers est en expansion tout le temps ça marche. Mais avec par exemple (Mm=0.7,Mr=0,Lambda=-1) l'univers s'étend puis se recontracte (euh oui je ne l'avais pas dit mais a(x) c'est le facteur d'échelle de l'univers et x est proportionnel au temps), on obtient juste la phase d'expansion.
En dessous et au dessus de certaines valeurs de x, Maple n'a pas l'air de s'en sortir (il me sort une "singularity").
Faut trouver une arme encore plus fatale que cela...