ATLAS Ouvrir dans l'atlas →

ExerciceTerminale, chapitre 8 — Équations différentielles · Pour aller plus loin

Exercice 35

DémonstrationCalcul

La méthode d'Euler. On considère l'équation y′=yy'=y avec la condition y(0)=1y(0)=1, dont la solution est la fonction exponentielle. On cherche à en calculer une valeur approchée en un réel x>0x>0 sans utiliser l'exponentielle. Pour cela, on fixe un entier n⩾1n\geqslant 1, on pose h=xnh=\dfrac{x}{n}, et l'on construit les nombres y0y_{0}, y1y_{1}, …, yny_{n} par y0=1etyk+1=yk+h ykpour 0⩽k⩽n−1y_{0}=1 \qquad\text{et}\qquad y_{k+1}=y_{k}+h\,y_{k} \quad\text{pour } 0\leqslant k\leqslant n-1

1. Expliquer pourquoi yky_{k} est une valeur approchée de ekh\mathrm{e}^{kh}. On pourra utiliser l'approximation d'une fonction par sa tangente, f(t+h)≈f(t)+h f′(t)f(t+h)\approx f(t)+h\,f'(t).

2. Montrer que yn=(1+xn)ny_{n}=\left(1+\dfrac{x}{n}\right)^{n}.

3. En déduire, à l'aide de l'exercice 32 du chapitre 5, que yny_{n} tend vers ex\mathrm{e}^{x} quand nn tend vers +∞+\infty.

4. Pour x=1x=1, calculer yny_{n} pour n=10n=10 et n=1 000n=1\,000, et comparer à e\mathrm{e}. Que peut-on dire de l'erreur commise ?

5. Reprendre la question 2 pour l'équation y′=ayy'=ay, y(0)=1y(0)=1, où aa est un réel donné.

Source : d'après le poly « Mathématiques : du lycée aux CPGE scientifiques » (LLG / H4), § 5.2, note 21 (méthode d'Euler)

Corrigé

1. Pour hh petit, la courbe d'une fonction dérivable est proche de sa tangente : f(t+h)≈f(t)+h f′(t)f(t+h)\approx f(t)+h\,f'(t). Si ff est solution de y′=yy'=y, cela devient f(t+h)≈f(t)+h f(t)f(t+h)\approx f(t)+h\,f(t). Partant de f(0)=1=y0f(0)=1=y_{0}, la relation yk+1=yk+hyky_{k+1}=y_{k}+hy_{k} reproduit cette approximation de proche en proche : y1≈f(h)y_{1}\approx f(h), y2≈f(2h)y_{2}\approx f(2h), et ainsi de suite, yky_{k} étant une valeur approchée de f(kh)=ekhf(kh)=\mathrm{e}^{kh}.

2. La relation s'écrit yk+1=(1+h) yky_{k+1}=(1+h)\,y_{k} : la suite (yk)(y_{k}) est géométrique de raison 1+h1+h et de premier terme 11. Donc yn=(1+h)ny_{n}=(1+h)^{n}, et comme h=xnh=\frac{x}{n} : yn=(1+xn)n\boxed{y_{n}=\left(1+\frac{x}{n}\right)^{n}}

3. D'après l'exercice 32 du chapitre 5, (1+xn)n\left(1+\frac{x}{n}\right)^{n} tend vers ex\mathrm{e}^{x}. lim⁡n→+∞yn=ex\boxed{\lim_{n\to+\infty}y_{n}=\mathrm{e}^{x}} La méthode d'Euler converge donc vers la vraie solution lorsque le pas hh tend vers 00.

4. Pour x=1x=1 :

nynerreur e−yn102,593 740,124 51 0002,716 920,001 36\begin{array}{|c|c|c|} \hline n & y_{n} & \text{erreur }\mathrm{e}-y_{n} \\ \hline 10 & 2{,}593\,74 & 0{,}124\,5 \\ \hline 1\,000 & 2{,}716\,92 & 0{,}001\,36 \\ \hline \end{array}

L'approximation est toujours inférieure à e\mathrm{e} : la courbe de l'exponentielle, convexe, est au-dessus de ses tangentes (chapitre 3), si bien que chaque pas sous-estime la croissance. Quand nn est multiplié par 100100, l'erreur est divisée par environ 100100 : elle est à peu près proportionnelle au pas h=1nh=\frac{1}{n}, et vaut environ e2n\frac{\mathrm{e}}{2n}. Pour gagner une décimale, il faut donc dix fois plus de calculs, ce qui rappelle la lenteur observée à l'exercice 32 du chapitre 5.

5. Pour y′=ayy'=ay, la relation devient yk+1=yk+h a yk=(1+ah) yky_{k+1}=y_{k}+h\,a\,y_{k}=(1+ah)\,y_{k}, d'où yn=(1+axn)n\boxed{y_{n}=\left(1+\frac{ax}{n}\right)^{n}} qui tend vers eax\mathrm{e}^{ax}, toujours d'après l'exercice 32 du chapitre 5.

Remarque. La méthode d'Euler s'applique à des équations que l'on ne sait pas résoudre explicitement, comme l'équation logistique de l'exercice 4 : c'est ainsi que travaillent les logiciels de simulation, avec des variantes plus précises.

Dans l'édifice

Thèmes

équation différentielle, équation y' = ay, condition initiale, méthode d'euler

Notions

équation différentielle, exponentielle

Ouvrir dans l'atlas interactif →