Base de travail : le TD3 de 2015-16 (Y. Demay)
1. Méthode solve existante dans Sagemath ?
Lire le livre Sagemath, le tutoriel, faire une recherche sur internet "sagemath equations différentielles" ou "sagemath solve differential equations"
Exemple avec l'équation différentielle et la condition initiale .
2. Implémentatation de la méthode d'Euler
Ecrire une fonction prenant comme argument une fonction , les nombres , , , , rendant une liste des coordonnées où et vérifie la relation de récurrence
Dessiner la réunion des segments reliant deux points consécutifs
Essayer avec , , , différentes valeurs de ,
Une équation différentielle d'ordre 2 est une équation de la forme
On peut la convertir en une équation diff. vectorielle d'ordre 1 : en posant L'équation différentielle vectorielle d'ordre 1 admet une solution approchée par la méthode d'Euler exposée plus haut.1. Ecrire une fonction prenant comme argument une fonction , les nombres , rendant une liste de triplets correspondant par la méthode d'Euler à une approximation de la solution de l'équation vérifiant les conditions initiales et .
Ecrire une fonction dessinant la réunion des segments reliant deux points consécutifs.
2. Tester les deux fonctions avec l'équation
Qu'obtient t-on si on varie ?
La fonction sagemath desolve donne t-elle une solution formelle de cette équation ?
Comment se compare le graphe de la solution formelle avec le graphe de la solution approchée de l'équation pour puis ?