Source Latex
\documentclass[12pt]{article}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage[french]{babel}
\usepackage{amsmath}
\usepackage[utf8]{inputenc}
\usepackage{calc}
\usepackage{enumerate}
\usepackage{array}
\usepackage{pst-all}
\usepackage{hyperref}
\hypersetup{
pdfauthor={Yoann Morel},
pdfsubject={Construction de la fonction exponentielle: méthode d'Euler},
pdftitle={Méthode d'Euler pour construire la fonction exponentielle},
pdfkeywords={Méthode d'Euler}
}
\hypersetup{
colorlinks = true,
linkcolor = red,
anchorcolor = red,
citecolor = blue,
filecolor = red,
urlcolor = red
}
\voffset=-1cm
% Raccourcis diverses:
\newcommand{\dsp}{\displaystyle}
\newcommand{\la}{\left\{}\newcommand{\ra}{\right\}}
\newcommand{\lp}{\left(}\newcommand{\rp}{\right)}
\newcommand{\lb}{\left[}\newcommand{\rb}{\right]}
\def\N{{\rm I\kern-.1567em N}}
\def\R{{\rm I\kern-.1567em R}}
\title{Construction de la fonction exponentielle -- Méthode d'Euler}
\author{Y. Morel}
\date{}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}
\maketitle
\vspace*{2em}
La fonction exponentielle est l'unique fonction dérivable sur $\R$ (cf. cours)
qui vérifie $f'=f$ et $f(0)=1$.
\medskip
\`A partir de ces deux seules informations, on peut construire
approximativement la fonction et sa courbe,
c'est-à-dire qu'on peut calculer approximativement la valeur $f(x)$ pour tout nombre réel $x$.
\section{Méthode d'Euler: approximation des valeurs de $f$}
L'information principale sur $f$ est une information sur sa dérivée $f'$.
Or la dérivée de $f$ en $x$
est la limite du taux de variation:
\[f'(x)=\lim_{h\to0}\dfrac{f(x+h)-f(x)}h\]
Ainsi, pour un nombre $h$ "assez petit" on peut écrire l'approximation
\[f'(x)\simeq\dfrac{f(x+h)-f(x)}h\]
Plus $h$ est petit, meilleure est cette approximation.\\
Maintenant, comme $f'=f$, on a donc
\[f(x)\simeq\dfrac{f(x+h)-f(x)}h\]
soit, en isolant $f(x+h)$,
\[f(x+h)\simeq hf(x)+f(x)\]
ou encore
\[f(x+h)\simeq(1+h)f(x)\]
En particulier, en partant de $x=0$, on obtient
\[f(h)=(1+h)f(0)\]
et comme $f(0)=1$,
\[f(h)=(1+h)\]
On a donc
\begin{itemize}
\item pour $x=0$, $f(h)=(1+h)$
\item pour $x=h$, $f(h+h)=f(2h)\simeq(1+h)f(h)=(1+h)^2$,
\item pour $x=2h$,$f(2h+h)=f(3h)\simeq(1+h)f(2h)=(1+h)^3$,
\item \dots \dots
\item pour tout entier $n$, $f(nh)\simeq(1+h)^n$
\end{itemize}
On balaye ainsi, par pas de $h$, les valeurs prises par la fonction $f$.
\bigskip
Pour calculer une valeur particulière $f(x)$, il suffit d'adapter ce pas:
si $x$ est un nombre réel quelconque, et $n$ un entier naturel,
il suffit de poser $h=\dfrac{x}n$ pour avoir $x=nh$,
et alors, avec le résultat précédent,
\[f(x)=f(nh)\simeq(1+h)^n\]
soit finalement
\[f(x)\simeq\lp1+\dfrac{x}n\rp^n\]
\section{Calcule de $\exp(1)$}
On peut maintenant donner une valeur approchée de $f(1)$.
On arrondissant les valeurs à $10^{-3}$ près, on obtient:
\begin{itemize}
\item Pour $n=10$, $f(1)=\lp1+\dfrac1{10}\rp^{10}\simeq 2,594$
\item Pour $n=100$, $f(1)=\lp1+\dfrac1{100}\rp^{100}\simeq 2,705$
\item Pour $n=1000$, $f(1)=\lp1+\dfrac1{1000}\rp^{1000}\simeq 2,717$
\item \dots
\end{itemize}
Valeurs qu'on compare bien s\^ur à la valeur de
$\exp(1)$, connue par toute calculatrice ou ordinateur**:
\[\exp(1)\simeq 2,718\]
\section{Courbe représentative de l'exponentielle}
On peut de m\^eme calculer toute valeur souhaitée $f(x)$ de la m\^eme fa\c con,
et tracer ainsi l'allure de la courbe représentative de la fonction exponentielle, obtenue comme courbe limite de la courbe de la fonction $f$ lorsque $n\to+\infty$.
\medskip
On trace donc la courbe de la "vraie"*** fonction exponentielle
et celles des fonctions définies par
$f(x)=\lp1+\dfrac{x}n\rp^n$ pour différentes valeurs de l'entier $n$:
\[\psset{xunit=1cm,yunit=.15cm,arrowsize=8pt}\begin{pspicture*}(-1.5,-5)(4,60)
\psline{->}(-1.4,0)(3.9,0)
\psline{->}(0,-.5)(0,56)
\multido{\i=-1+1}{5}{\psline(\i,-1)(\i,1)\rput(\i,-2.4){$\i$}}
\multido{\i=10+10}{5}{\psline(-.1,\i)(.1,\i)\rput[r](-.2,\i){$\i$}}
\psplot[linewidth=1.8pt,linecolor=red]{-1}{8}{2.71828 x exp}
\psplot[linewidth=1.8pt,linecolor=blue]{-1}{6}{1 x 10 div add 10 exp}
\psplot[linewidth=1.8pt,linecolor=cyan]{-1}{6}{1 x 50 div add 50 exp}
\psplot[linewidth=1.8pt,linecolor=magenta]{-1}{7.6}{1 x 5 div add 5 exp}
%\pspolgon(
\psline[linewidth=2.5pt,linecolor=red](.5,45)(1,45)\rput[l](1.2,44){$\exp$}
\psline[linewidth=2.5pt,linecolor=cyan](.5,41)(1,41)\rput[l](1.2,41){$n=50$}
\psline[linewidth=2.5pt,linecolor=blue](.5,37)(1,37)\rput[l](1.2,37){$n=10$}
\psline[linewidth=2.5pt,linecolor=magenta](.5,33)(1,33)\rput[l](1.2,33){$n=5$}
\end{pspicture*}\]
\vspace{4em}
**: Pour utiliser la fonction exponentielle en python (et d'autres comme le logarithme, les fonctions trigonométriques, \dots ), il faut la charger depuis la bibliothèque des fonctions mathéamtiques:
\textit{from math import exp}
\textit{exp(1)}
***: La fonction exponentielle n'a pas d'expression algébrique "simple" (comme la fonction carré, un polyn\^ome \dots). \\
La "vraie courbe" de l'exponentielle est en fait tracée en utilisant des calculs approchés, comme ceux menés ici. \\
Les méthodes de calcul des valeurs de l'exponentielle par une calculatrice ou un ordinateur sont juste plus efficaces: c'est-à-dire la précision meilleure qu'avec la méthode d'Euler et pour un nombre plus faible de calculs. \\
Trouver ce genre de méthodes plus efficaces est le champ mathématique qui s'appelle l'anlayse numérique.
\end{document}
Télécharger le fichier source 
