Calculateur de la Méthode d'Euler

Calculateur de la Méthode d'Euler

Le Calculateur de la Méthode d'Euler résout numériquement tout problème de valeur initiale de premier ordre de la forme \( y' = f(x, y), \; y(x_0) = y_0 \) en utilisant la méthode d'Euler classique (explicite). Il renvoie un tableau d'itération complet, trace le polygone d'Euler sur un champ de vecteurs en direct, compare la solution à trois tailles de pas différentes pour que vous puissiez visualiser la convergence de la méthode, et — si vous fournissez la solution exacte sous forme fermée — produit une analyse d'erreur par étape.

Qu'est-ce que la méthode d'Euler ?

La méthode d'Euler est l'algorithme le plus simple pour approximer la solution d'un problème de valeur initiale. À partir d'un point connu \( (x_0, y_0) \) sur la courbe de solution, elle progresse de manière répétée par un petit pas de taille h le long de la pente locale \( f(x, y) \) :

Géométriquement, chaque étape est un court segment rectiligne dont la pente est égale à la valeur de l'équation différentielle au point actuel. La ligne brisée résultante — le polygone d'Euler — est une approximation de la solution réelle (généralement courbe).

Quelle est sa précision ?

La méthode d'Euler est une méthode du premier ordre. L'erreur de troncature locale à chaque étape est \( O(h^2) \) et l'erreur globale après intégration sur un intervalle fixe est \( O(h) \). En pratique :

• Réduire de moitié la taille du pas divise approximativement par deux l'erreur globale.
• L'erreur croît linéairement avec la longueur de l'intervalle d'intégration.
• L'erreur est maximale là où la solution présente une forte courbure.

La comparaison intégrée de la taille du pas (h, h/2, h/4) vous permet de voir directement cette convergence linéaire : activez l'option et vérifiez que les trois valeurs finales s'approchent d'une limite commune, chaque valeur étant environ deux fois plus proche de la limite que la précédente.

Lire le graphique

La visualisation superpose quatre types d'informations sur un seul plan de coordonnées :

• Champ de vecteurs gris — de courts segments de ligne dont l'inclinaison est égale à \( f(x, y) \) en ce point. Considérez-le comme "le flux dicté par l'EDO". Toute courbe de solution doit être tangente au champ en chaque point.
• Polygone d'Euler indigo — la solution numérique par étapes. Chaque segment part du point précédent de la grille et pointe selon \( f(x_n, y_n) \) sur une distance de h.
• Courbe exacte en pointillés verts — présente uniquement lorsque vous fournissez la solution sous forme fermée. Les pointillés oranges verticaux sont les erreurs locales signées \( y_n - y_{\text{exact}}(x_n) \).
• Courbes de comparaison orange et verte — le même problème ré-exécuté à h/2 et h/4, affiché lorsque la comparaison de taille de pas est activée.

Comment utiliser ce calculateur

1. Entrez le membre de droite de l'EDO dans le champ marqué y' =. Utilisez x et y comme variables. Les opérateurs supportés sont + − × ÷ ^, et les fonctions supportées incluent sin, cos, tan, asin, acos, atan, exp, ln, log, log10, log2, sqrt, abs.
2. Définissez les conditions initiales : la valeur de départ x₀, le y₀ initial à ce point, la taille du pas h (positif pour intégrer vers l'avant, négatif pour intégrer vers l'arrière), et le nombre d'étapes n.
3. (Optionnel) Fournissez la solution exacte y(x) si vous la connaissez. Le calculateur calculera \( |y_n - y(x_n)| \) à chaque étape et indiquera les erreurs maximale et finale.
4. Basculez les options de visualisation : le champ de vecteurs est activé par défaut ; la comparaison de taille de pas superpose deux courbes supplémentaires à h/2 et h/4.
5. Cliquez sur Exécuter. La section des résultats affiche les statistiques résumées, le graphique, un panneau de comparaison de convergence et le tableau d'itération complet. Le survol d'une ligne met en évidence le point correspondant sur le graphique (et vice versa).

Exemple concret

Considérons \( y' = x + y, \; y(0) = 1 \) avec h = 0,1 et 10 étapes. La solution exacte est \( y(x) = -x - 1 + 2e^x \). L'application de la formule d'Euler donne :

L'erreur finale est d'environ 0,249. Diviser h par deux (0,05) ramène l'erreur finale à environ 0,13, et diviser à nouveau par deux (0,025) la ramène à environ 0,067 — une convergence linéaire nette, exactement comme le prédit la théorie.

Méthode d'Euler vs autres méthodes numériques

Quand Euler fait défaut

La méthode d'Euler explicite peut mal se comporter dans trois situations :

• Taille de pas trop grande — le polygone oscille ou diverge. La solution est de réduire h ; la comparaison h, h/2, h/4 rend cela instantanément visible.
• EDO raides — des équations avec des modes à décroissance rapide et lente simultanés obligent h à être minuscule pour la stabilité. Passez à une méthode implicite (Euler rétrograde) ou BDF.
• Singularités dans f(x, y) — une division par zéro, une sqrt d'un nombre négatif ou un ln d'un nombre non positif arrêtera l'intégration. Le calculateur signale clairement l'étape incriminée.

Applications courantes

• Physique — Deuxième loi de Newton sous forme de système de premier ordre, désintégration radioactive \( \dot{N} = -\lambda N \), loi du refroidissement de Newton.
• Biologie et épidémiologie — Croissance logistique \( \dot{y} = r\,y(1 - y/K) \), modèles compartimentaux SIR.
• Économie — Intérêts composés continus, modèles de croissance de Solow simples.
• Chimie — Cinétique de réaction de premier ordre \( \dot{c} = -k c \).
• Enseignement — Introduction du concept d'intégration numérique avant de passer à RK4 ou aux solveurs adaptatifs.

Foire Aux Questions

Qu'est-ce que la méthode d'Euler ?

La méthode d'Euler est la procédure numérique la plus simple pour résoudre un problème de valeur initiale y' = f(x, y), y(x0) = y0. À chaque étape, elle fait progresser la solution par y_{n+1} = y_n + h · f(x_n, y_n), suivant effectivement la pente au point actuel sur une courte distance h. Elle est précise au premier ordre, ce qui signifie que l'erreur globale est O(h).

Quelle est la précision de la méthode d'Euler ?

La méthode d'Euler présente une erreur de troncature locale O(h²) et une erreur globale O(h). Diviser par deux la taille du pas réduit approximativement de moitié l'erreur globale. C'est pourquoi la comparaison de convergence à h, h/2 et h/4 dans ce calculateur est si instructive : vous pouvez voir l'erreur diminuer de manière approximativement linéaire avec h.

Quand la méthode d'Euler échoue-t-elle ?

La méthode d'Euler peut devenir instable pour les problèmes raides ou lorsque la taille du pas est trop grande par rapport à la courbure locale de la solution. Vous pourriez voir la solution numérique osciller, exploser vers l'infini ou s'écarter visiblement de la solution réelle. Réduire h aide généralement ; pour les équations raides, les méthodes implicites telles que la méthode d'Euler implicite sont préférées.

Comment choisir la taille du pas ?

Commencez par un h qui donne environ 10 à 50 étapes sur l'intervalle d'intérêt. Si le polygone d'Euler dévie visiblement du champ de vecteurs ou de votre solution exacte, divisez h par deux et relancez. Utilisez la comparaison intégrée h, h/2, h/4 pour vérifier que les trois courbes convergent l'une vers l'autre.

Quelle est la différence entre la méthode d'Euler et Runge-Kutta (RK4) ?

La méthode de Runge-Kutta d'ordre quatre évalue la pente en quatre points par étape et les combine avec des poids (1, 2, 2, 1)/6, donnant une erreur globale O(h⁴) — plusieurs ordres de grandeur plus précise que l'O(h) d'Euler pour le même nombre d'étapes. Euler reste précieux pour l'enseignement du concept d'intégration numérique et pour des applications très simples ou de faible précision.

Puis-je utiliser ceci pour des systèmes d'EDO ?

Ce calculateur gère une seule EDO scalaire de premier ordre y' = f(x, y). Pour les systèmes ou pour les EDO d'ordre supérieur, vous pouvez réécrire l'équation sous forme de système de premier ordre et utiliser un solveur de système dédié, ou convertir une équation de second ordre en deux équations de premier ordre et les résoudre composante par composante.

Puis-je intégrer vers l'arrière dans le temps ?

Oui — entrez une taille de pas h négative. Le calculateur progressera à partir de x₀ dans la direction négative pendant n étapes. Ceci est utile pour reconstruire le passé à partir d'un état présent connu.

Lectures complémentaires

• Méthode d'Euler — Wikipédia
• Méthodes de Runge–Kutta — Wikipédia
• Champ de vecteurs — Wikipédia
• Équation différentielle raide — Wikipédia