Calculatrice Méthode d Euler ODE (Explicite, Implicite, Heun)

Calculatrice méthodes d Euler pour ODE y' = f(t,y): Euler explicite, implicite, Heun et point median. Convergence ordre 1-2, exemples logistique et Newton.

Paramètres de calcul
Ce calculateur vous a-t-il été utile ?

Méthodes d Euler pour équations differentielles ordinaires

Les méthodes d Euler sont la classe la plus simple d algorithmes pour resoudre numeriquement les équations differentielles ordinaires (ODE) de forme y' = f(t,y), y(t_0) = y_0. Developpees par Leonhard Euler en 1768, elles sont le point de depart de l analyse numérique des ODE.

Euler explicite (Forward Euler)

Formule: y_{n+1} = y_n + h f(t_n, y_n) ou h = (T - t_0)/N est le pas. Convergence globale ordre 1: erreur = O(h). Rapide mais imprecis. Instable pour problèmes stiff.

Euler implicite (Backward Euler)

Formule: y_{n+1} = y_n + h f(t_{n+1}, y_{n+1}). Necessite resoudre équation implicite a chaque pas. Coût double mais A-stable: pas de limite sur h. Idéal pour ODE stiff.

Méthode de Heun

Predicteur-correcteur: predicteur Euler explicite, correcteur trapezoidal. Ordre 2: erreur = O(h^2). Coût: 2 evaluations par pas. Compromis optimal précision/coût.

Point median (Modified Euler)

k1 = f(t_n, y_n), y_{n+1} = y_n + h f(t_n + h/2, y_n + h/2 k1). Ordre 2 comme Heun. Stable pour ODE oscillantes. Variante Runge-Kutta deux etages.

Convergence et ordre

Erreur troncature locale Euler explicite: -h^2/2 y''(xi). Sommant N = (T-t0)/h pas: erreur globale O(h). Heun et point median: O(h^2). RK4 classique: O(h^4). Compromis: ordre supérieur = plus d evaluations par pas, h plus grand possible.

Stabilité

Pour y' = lambda y, Euler explicite stable si |1 + h lambda| < 1. Implicite A-stable toujours. Pour ODE stiff implicite essentiel.

Problèmes stiff

ODE stiff ont échelles temporelles tres differentes. Explicite necessite h minuscule. Implicite n a pas cette limite. Exemple: équation Robertson en cinétique chimique.

ODE remarquables

  • Décroissance exponentielle: y' = -k y.
  • Croissance logistique: y' = r y (1 - y/K).
  • Newton cooling: T' = -k (T - T_amb).
  • Riccati: y' = y^2 - t.
  • Lotka-Volterra: proie-predateur.
  • Oscillateur harmonique: y'' + y = 0.

Quand utiliser chaque méthode

  • Euler explicite: didactique, non stiff, h petit acceptable.
  • Euler implicite: stiff (chimie, fluides, equilibres).
  • Heun: code educatif, précision moyenne.
  • RK4: standard pratique pour non stiff.
  • Adaptatifs (RK45): vitesses variables.

Histoire

Euler 1768 introduit méthode. Cauchy 1840 prouve rigoureusement convergence. Heun 1900 et Kutta 1901 etendent aux ordres superieurs. Dahlquist 1956 introduit A-stabilité. Logiciels modernes (LSODE, scipy.integrate.solve_ivp) incluent dizaines de méthodes specialisees.

Erreur locale et globale : d'ou vient l'ordre

La distinction entre erreur locale et globale explique pourquoi l'Euler explicite est dit methode d'ordre un. L'erreur locale de troncature est l'erreur commise en un seul pas, en supposant exacte la valeur de depart : en developpant la solution en serie de Taylor on voit qu'Euler neglige le terme (h^2/2) y'', donc son erreur locale est proportionnelle a h^2. L'erreur globale, accumulee sur tout l'intervalle divise en environ 1/h pas, est en revanche proportionnelle a h : diviser le pas par deux divise par deux l'erreur finale.

Les methodes d'ordre deux, comme Heun et le point median, annulent aussi le terme quadratique et ont une erreur locale h^3 et globale h^2 : diviser h par deux reduit l'erreur au quart. Un moyen pratique d'estimer l'erreur sans connaitre la solution exacte est l'extrapolation de Richardson, qui compare les resultats obtenus avec deux pas differents (h et h/2) pour deduire l'ordre effectif et une estimation quantitative de l'erreur.

Vers Runge-Kutta : l'echelle des methodes a un pas

Les methodes d'Euler sont les premiers echelons d'une hierarchie plus large, celle des methodes de Runge-Kutta a un pas. L'Euler explicite est le Runge-Kutta d'ordre un ; le point median et la methode de Heun sont des variantes d'ordre deux qui evaluent la pente en des points intermediaires pour doubler la precision. Le celebre RK4, d'ordre quatre, combine quatre evaluations par pas et offre un excellent compromis entre precision et cout, et est le cheval de bataille de la resolution numerique d'EDO.

L'etape suivante est le controle adaptatif du pas : des methodes comme Runge-Kutta-Fehlberg (RKF45) calculent simultanement une estimation d'ordre quatre et une d'ordre cinq, et de leur difference deduisent l'erreur locale, adaptant automatiquement h pour le maintenir sous un seuil. Ainsi le solveur utilise de grands pas la ou la solution est lisse et de petits pas la ou elle varie rapidement, optimisant le nombre d'evaluations : une idee qui descend directement des simples methodes d'Euler.

Mis à jour en juillet 2026

Questions fréquentes

Qu'est-ce que la methode d'Euler ?

C'est la methode numerique la plus simple pour resoudre une equation differentielle ordinaire y' = f(t,y) avec condition initiale. Elle approche la solution par la formule y_{n+1} = y_n + h f(t_n, y_n), avancant par petits pas h selon la direction de la pente donnee par l'equation.

Quelle difference entre Euler explicite et implicite ?

L'Euler explicite evalue la pente au debut du pas, y_{n+1} = y_n + h f(t_n,y_n), et est simple mais peu stable. L'implicite l'evalue a la fin, y_{n+1} = y_n + h f(t_{n+1},y_{n+1}), exigeant de resoudre une equation a chaque pas mais offrant une stabilite bien superieure pour les problemes raides.

Quel ordre a la methode d'Euler ?

L'Euler explicite est d'ordre un : l'erreur globale est proportionnelle au pas h, donc diviser h par deux divise l'erreur par deux. Les methodes de Heun et du point median sont d'ordre deux (erreur proportionnelle a h^2), bien plus precises a pas egal.

Qu'est-ce qu'un probleme raide (stiff) ?

C'est une equation ou coexistent des echelles de temps tres differentes, avec des composantes qui decroissent rapidement a cote d'autres lentes. Les methodes explicites comme Euler exigent des pas minuscules pour rester stables, tandis que les methodes implicites gerent ces problemes avec des pas bien plus grands.

Qu'est-ce que la methode de Heun ?

C'est une methode d'ordre deux aussi appelee Euler ameliore : elle utilise Euler comme predicteur pour estimer y_{n+1}, puis corrige en moyennant les pentes au debut et a la fin du pas (regle des trapezes). Elle double l'ordre de precision par rapport a Euler au cout d'une evaluation supplementaire.

Comment estimer l'erreur sans la solution exacte ?

Avec l'extrapolation de Richardson : on resout le probleme avec un pas h et un pas h/2 et on compare les resultats. La difference, combinee a l'ordre connu de la methode, fournit une estimation quantitative de l'erreur et est la base du controle adaptatif du pas.

Quel rapport avec Runge-Kutta ?

Les methodes d'Euler sont les cas les plus simples de la famille Runge-Kutta a un pas : Euler est d'ordre un, Heun et le point median d'ordre deux, le celebre RK4 d'ordre quatre. En augmentant les evaluations de pente par pas on obtient une plus grande precision.

Quand utiliser quelle methode ?

Euler explicite pour des problemes simples et pedagogiques ; Heun ou point median quand plus de precision est necessaire a faible cout ; Euler implicite pour les problemes raides ou la stabilite est cruciale. Pour des applications serieuses on passe a RK4 ou a des methodes adaptatives comme RKF45.

Commentaires (0)

Pas encore de commentaire. Soyez le premier !