Formule D'Euler-Maclaurin

Formule D'Euler-Maclaurin

Formule d'Euler-Maclaurin

En mathématiques, la formule d'Euler-Maclaurin est une relation entre sommes discrètes et intégrales. Elle fut découverte indépendamment, aux alentours de 1735, par le mathématicien suisse Leonhard Euler et l'écossais Colin Maclaurin. Elle peut être utilisée pour approcher des intégrales par un procédé discret, par exemple dans la méthode des trapèzes ou celle de Romberg, ou à l'inverse pour transformer une somme discrète (finie ou non) et lui appliquer les techniques du calcul infinitésimal.

Énoncé

Soient deux entiers relatifs p et q. Pour une fonction 2k fois continûment dérivable sur le segment [p,q], la formule s'énonce ainsi :

\frac{f\left( p\right) +f\left( q\right) }{2}+\sum_{j=p+1}^{q-1}f\left(
j\right) =\int_p^q f(x)\,dx
+\sum_{j=1}^k\frac{b_{2j}}{(2j)!}\left(f^{(2j-1)}(q)-f^{(2j-1)}(p)\right)+R

avec :

 R = - \int_p^q f^{(2k)}(x) {B_{2k}(x-\lfloor x \rfloor) \over (2k)!}\,dx,

La notation Bi désigne le i-ème polynôme de Bernoulli, et B_i(x-\lfloor x \rfloor) en est une version périodisée. Les nombres bi désignent les nombres de Bernoulli : b1 = −1/2, b2 = 1/6, b3 = 0, b4 = −1/30, b5 = 0, b6 = 1/42, b7 = 0, b8 = −1/30.

Bien sûr, un simple changement de variable permet d'obtenir une formule analogue pour une fonction définie sur un segment à bornes non entières.

Démonstration

On se contentera de faire la démonstration sur l'intervalle [n,n + 1] avec n \in \mathbb{Z}  ; la formule précédente s'en déduit par sommation.

Soit g une fonction continûment dérivable sur [n,n + 1] . En utilisant la propriété des polynômes de Bernoulli :  \forall k \in \mathbb{N} B_{k+1}' = \left(k+1\right) B_{k} , on trouve en faisant une intégration par parties :

\int_n^{n+1} g \left( t \right) B_k \left( t-n \right) dt = \left[ \frac{g \left( t \right) B_{k+1} \left( t-n \right)}{k+1} \right]_n^{n+1} - \frac{1}{k+1} \int_n^{n+1} g' \left( t \right) B_{k+1} \left( t-n \right) dt

Or, sachant que pour  k \ge 2 , on a B_k \left( 1 \right) = B_k \left( 0 \right) = b_k , on en déduit :

\int_n^{n+1} g \left( t \right) B_k \left( t-n \right) dt = \frac{b_{k+1}}{k+1} \left( g \left( n+1 \right) - g \left( n \right) \right) - \frac{1}{k+1} \int_n^{n+1} g'\left( t\right) B_{k+1} \left( t-n \right) dt

Par récurrence sur k de 0 à 2p, en prenant g = f(2p), on obtient :

\int_n^{n+1} f \left( t \right)  dt = \frac{f\left( n\right) +f\left( n+1\right) }{2}+\sum_{k=2}^{2p} \frac{\left( -1 \right)^{k-1} b_k}{k!} \left( f^{(k-1)}\left(n+1\right) - f^{(k-1)}\left(n\right) \right) + \frac{1}{(2p)!} \int_n^{n+1} f^{(2p)} \left( t\right) B_{2p} \left( t-n \right) dt

Enfin, avec la propriété : \forall k \ge 1 , b_{2k+1} = 0 , on en déduit :

\int_n^{n+1} f \left( t \right)  dt = \frac{f\left( n\right) +f\left( n+1\right) }{2}+\sum_{k=2}^{\lfloor \frac{p}{2} \rfloor} \frac{b_{2k}}{(2k)!} \left( f^{(2k-1)}\left(n+1\right) - f^{(2k-1)}\left(n\right) \right) + \frac{1}{(2p)!} \int_n^{n+1} f^{(2p)} \left( t\right) B_{2p} \left( t-n \right) dt

Références

  • Analyse numérique et équations différentielles, J.-P. Demailly, Presses universitaires de Grenoble
  • Portail des mathématiques Portail des mathématiques
Ce document provient de « Formule d%27Euler-Maclaurin ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Formule D'Euler-Maclaurin de Wikipédia en français (auteurs)

Игры ⚽ Нужно сделать НИР?

Regardez d'autres dictionnaires:

  • Formule d'euler-maclaurin — En mathématiques, la formule d Euler Maclaurin est une relation entre sommes discrètes et intégrales. Elle fut découverte indépendamment, aux alentours de 1735, par le mathématicien suisse Leonhard Euler et l écossais Colin Maclaurin. Elle peut… …   Wikipédia en Français

  • Formule d'Euler-Maclaurin — En mathématiques, la formule d Euler Maclaurin (appelée parfois formule sommatoire d Euler) est une relation entre sommes discrètes et intégrales. Elle fut découverte indépendamment, aux alentours de 1735, par le mathématicien suisse Leonhard… …   Wikipédia en Français

  • Formule De Stirling — Pour les articles homonymes, voir Stirling. La formule de Stirling, du nom du mathématicien James Stirling, donne un équivalent de la factorielle au voisinage de l infini réel (quand n tend vers l infini) : que l on trouve souvent …   Wikipédia en Français

  • Formule de stirling — Pour les articles homonymes, voir Stirling. La formule de Stirling, du nom du mathématicien James Stirling, donne un équivalent de la factorielle au voisinage de l infini réel (quand n tend vers l infini) : que l on trouve souvent …   Wikipédia en Français

  • Euler — Leonhard Euler « Euler » redirige ici. Pour les autres significations, voir Euler (homonymie). Leonhard Euler …   Wikipédia en Français

  • Formule de Stirling — Pour les articles homonymes, voir Stirling. La formule de Stirling, du nom du mathématicien écossais James Stirling, donne un équivalent de la factorielle au voisinage de l infini réel (quand n tend vers l infini) : que l on trouve souvent… …   Wikipédia en Français

  • EULER (L.) — Avec Joseph Louis Lagrange, son émule plus jeune, Leonhard Euler est l’un des deux géants mathématiques qui ont dominé la science du XVIIIe siècle. Ses travaux, d’une abondance inégalée, couvrent tout le champ des mathématiques, de la mécanique… …   Encyclopédie Universelle

  • Leonard Euler — Leonhard Euler « Euler » redirige ici. Pour les autres significations, voir Euler (homonymie). Leonhard Euler …   Wikipédia en Français

  • Leonhard Euler — « Euler » redirige ici. Pour les autres significations, voir Euler (homonymie). Leonhard Euler Portrait par Johann Georg Brucker Naissance …   Wikipédia en Français

  • Leonhard Paul Euler — Leonhard Euler « Euler » redirige ici. Pour les autres significations, voir Euler (homonymie). Leonhard Euler …   Wikipédia en Français

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”