En matemáticas, la fórmula de Euler-Maclaurin relaciona a integrales con series. Esta fórmula puede ser usada para aproximar integrales por sumas finitas o, de forma inversa, para evaluar series (finitas o infinitas) resolviendo integrales. La fórmula fue descubierta independientemente por Leonhard Euler y Colin Maclaurin en 1735. Euler usó esta fórmula para calcular valores de series infinitas con convergencia lenta y Maclaurin la utilizó para calcular integrales.
La fórmula
Si z es un número correlacional y
es una función suave (suficientemente derivable) definida
, entonces, la integral
![{\displaystyle I=\int _{0}^{n}f(x)\,dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/899c923440fe4d1a416850c6e0b0395101102f86)
puede ser aproximada por la siguiente suma:
![{\displaystyle S={\frac {f\left(0\right)}{2}}+f\left(1\right)+\cdots +f\left(n-1\right)+{\frac {f\left(n\right)}{2}}={\frac {f\left(0\right)+f\left(n\right)}{2}}+\sum _{k=1}^{n-1}f\left(k\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d87a8b5938e1cf7dccbb19a54ffb317851e348b7)
(ver regla del trapecio). La fórmula de Euler-Maclaurin nos da una expresión para la diferencia entre la suma y la integral en función de derivadas de
en los extremos del intervalo de integración (0 y n). Para cualquier entero positivo p, tenemos que se cumple:
![{\displaystyle S-I=\sum _{k=1}^{p}{\frac {B_{k+1}}{(k+1)!}}\left(f^{(k)}(n)-f^{(k)}(0)\right)+R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7301fb9b3b1202a76b422ab8efbfc04583dfd8de)
donde
son los números de Bernoulli y R es una estimación del error normalmente pequeña.
Realizando un cambio de variable en la integral, se puede modificar esta fórmula para funciones
definidas en otros intervalos de la recta real.
El término de error
El término de error R es:
![{\displaystyle R=(-1)^{p}\int _{0}^{n}f^{(p+1)}(x){B_{p+1}(x-\lfloor x\rfloor ) \over (p+1)!}\,dx,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/38991aa6fcb86b30250c43328eba416f4dbd03f6)
donde
son los polinomios de Bernoulli periódicos. El término de error se puede acotar por:
![{\displaystyle \left|R\right|\leq {\frac {2}{(2\pi )^{2p}}}\int _{0}^{n}\left|f^{(p+1)}(x)\right|\,dx.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8dc86615fdab8b9508a253e598613e87ddf7b908)
Usos
Sumas de polinomios
Si
es un polinomio y p es suficientemente grande, entonces el término de error R se anula, por lo que se pueden resolver series de polinomios de forma exacta. Por ejemplo, si
, escogiendo p = 2 se obtiene:
![{\displaystyle \sum _{i=0}^{n}i^{3}=\left({\frac {n(n+1)}{2}}\right)^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/999bbaa6239c87a130eb92b9c9c74d4d1771c40c)
(ver fórmula de Faulhaber).
Integración numérica
La fórmula de Euler-Maclaurin se usa también para el análisis de errores en integraciones numéricas, de hecho, los métodos de extrapolación se basan en esta fórmula.
Expansión asintótica de series
Cuando se quiere calcular la expansión asintótica de series, la forma más cómoda de la fórmula de Euler-Maclaurin es:
![{\displaystyle \sum _{n=a}^{b}f(n)\sim \int _{a}^{b}f(x)\,dx+{\frac {f(a)+f(b)}{2}}+\sum _{k=1}^{\infty }\,{\frac {B_{2k}}{(2k)!}}\left(f^{(2k-1)}(b)-f^{(2k-1)}(a)\right)\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd074803fd60d43f4dc9f4309423dcb8c9497c0b)
donde
y
son enteros. Puede ocurrir que esta fórmula siga siendo válida incluso tomando el límite
o
, o ambos. En muchos casos, la integral de la derecha es resoluble mediante funciones elementales de forma cerrada incluso cuando la serie de la izquierda no puede ser resuelta. Entonces, todos los términos de la serie asintótica pueden ser expresados mediante funciones elementales, por ejemplo:
![{\displaystyle \sum _{k=0}^{\infty }{\frac {1}{(z+k)^{2}}}\sim \underbrace {\int _{0}^{\infty }{\frac {1}{(z+k)^{2}}}\,dk} _{=1/z}+{\frac {1}{2z^{2}}}+\sum _{t=1}^{\infty }{\frac {B_{2t}}{z^{2t+1}}}\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e1a3210a7ef7dae4b3c335dc7016d1dc56bb2b43)
Donde la serie de la izquierda es igual a la suma de
y
, donde la serie de la derecha es la función poligamma de primer orden.
Restando
a los dos lados de la expresión, obtenemos una serie asintótica de
. De hecho, esta serie es el punto inicial de una de las posibles derivaciones de la fórmula de Stirling del factorial.
Demostración
Demostración por inducción matemática
Se seguirá la demostración que aparece en Apostol.[1]
Los polinomios de Bernoulli
con
se pueden definir recursivamente como sigue:
![{\displaystyle B_{0}(x)=1,\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/566c70fbefc86d0d98ab0f629c507cc951fc3fef)
![{\displaystyle B_{n}'(x)=nB_{n-1}(x){\mbox{ y }}\int _{0}^{1}B_{n}(x)dx=0{\mbox{ para }}n\geq 1.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/842f7c0aff08b547bac73238cd7adbd247975030)
Los primeros 4 son:
![{\displaystyle B_{1}(x)=x-1/2,\quad B_{2}(x)=x^{2}-x+1/6,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77f5f52190845b5a9b84e217606073ab41dd3465)
![{\displaystyle B_{3}(x)=x^{3}-{\frac {3}{2}}x^{2}+{\frac {1}{2}}x,\quad B_{4}(x)=x^{4}-2x^{3}+x^{2}-{\frac {1}{30}},\dots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/dac4a6c4446f2cbb758dc7dfa17df3cb2b85cf38)
Los valores
son los números de Bernoulli. Para
se cumple
.
Las funciones periódicas de Bernoulli
se definen como:
![{\displaystyle P_{n}(x)=B_{n}(x-\lfloor x\rfloor ){\mbox{ for }}0<x<1,\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8aa657ca7414d73e788a46aecd0e9ced9d22c12c)
Es decir, son iguales a los polinomios de Bernoulli en el intervalo (0,1), pero son funciones periódicas de periodo 1 en el resto del eje real.
Sea la integral
![{\displaystyle \int _{k}^{k+1}f(x)\,dx=\int u\,dv,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f2b836d26ec64810c87cfcc5940bef8131405ef)
donde
![{\displaystyle {\begin{aligned}u&{}=f(x),\\du&{}=f'(x)\,dx,\\dv&{}=P_{0}(x)\,dx\quad ({\mbox{con }}P_{0}(x)=1),\\v&{}=P_{1}(x).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ef03f79dd233a57f9e5b77b99a80f47e2e747f5e)
Integrando por partes obtenemos
![{\displaystyle {\begin{aligned}uv-\int v\,du&{}={\Big [}f(x)P_{1}(x){\Big ]}_{k}^{k+1}-\int _{k}^{k+1}f'(x)P_{1}(x)\,dx\\\\&{}={f(k)+f(k+1) \over 2}-\int _{k}^{k+1}f'(x)P_{1}(x)\,dx.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/07d550bba53a05f62adabbb43e748eff2014fc44)
Sumando desde
hasta
se obtiene:
![{\displaystyle \int _{1}^{n}f(x)\,dx={f(1) \over 2}+f(2)+\cdots +f(n-1)+{f(n) \over 2}-\int _{1}^{n}f'(x)P_{1}(x)\,dx.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a85ab387a8ab1671cb280781eb6e52db8d8253e)
Sumando
a ambos lados de la igualdad y reagrupando términos se obtiene:
![{\displaystyle \sum _{k=1}^{n}f(k)=\int _{1}^{n}f(x)\,dx+{f(1)+f(n) \over 2}+\int _{1}^{n}f'(x)P_{1}(x)\,dx.\qquad (1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/288cfa9f1a994c2dfc8873452ff9742ac0530b1d)
Por tanto, los dos últimos términos nos dan el error cuando la integral se toma como aproximación de la serie.
Consideremos ahora a la siguiente integral:
![{\displaystyle \int _{k}^{k+1}f'(x)P_{1}(x)\,dx=\int u\,dv,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/db67032ebdab4f28877a464d9dab248e6f707dd1)
donde
![{\displaystyle {\begin{aligned}u&{}=f'(x),\\du&{}=f''(x)\,dx,\\dv&{}=P_{1}(x)\,dx,\\v&{}=P_{2}(x)/2.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e285321f5c572ef45d7cf040332a0563f914ef3)
Integrando otra vez por partes se obtiene
![{\displaystyle {\begin{aligned}uv-\int v\,du&{}=\left[{f'(x)P_{2}(x) \over 2}\right]_{k}^{k+1}-{1 \over 2}\int _{k}^{k+1}f''(x)P_{2}(x)\,dx\\\\&{}={f'(k+1)-f'(k) \over 12}-{1 \over 2}\int _{k}^{k+1}f''(x)P_{2}(x)\,dx.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7561af0438d0769990766bcda8f0f16b220a323c)
Sumando desde
hasta
y reemplazando la última integral en (1) por el resultado que se acaba de obtener, tenemos:
![{\displaystyle \sum _{k=1}^{n}f(k)=\int _{1}^{n}f(x)\,dx+{f(1)+f(n) \over 2}+{f'(n)-f'(1) \over 12}-{1 \over 2}\int _{1}^{n}f''(x)P_{2}(x)\,dx.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11799fcfd9425a07287ef77c65098a7ce7b29bd1)
Obviamente, este procedimiento puede ser iterado. De esta manera se obtiene una demostración de la fórmula de Euler-Maclaurin por inducción, en la que los pasos de la inducción constan de una integración por partes y en el uso de las propiedades de las funciones periódicas de Bernoulli.
Para acotar el tamaño del error cuando la suma se aproxima por la integral, se tiene en cuenta que, en el intervalo
, los polinomios de Bernoulli alcanzan sus valores máximos absolutos en los puntos finales del intervalo (véase D.H. Lehmer en la referencias) y que
La fórmula de Euler-Maclaurin puede ser obtenida como una aplicación de algunas ideas de espacios de Hilbert y análisis funcional. Sea
un polinomio de Bernoulli. Un conjunto de funciones duales a los polinomios de Bernoulli está dado por
![{\displaystyle {\tilde {B}}_{n}(x)={\frac {(-1)^{n+1}}{n!}}\left[\delta ^{(n-1)}(1-x)-\delta ^{(n-1)}(x)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6133525d77d432d6e8f7c1776a96e59cc9c39ce9)
donde δ es la función delta de Dirac. Esta fórmula no es más que una notación formal de la idea de tomar derivadas en un punto, entonces se tiene
![{\displaystyle \int _{0}^{1}{\tilde {B}}_{n}(x)f(x)\,dx={\frac {1}{n!}}\left[f^{(n-1)}(1)-f^{(n-1)}(0)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfd85df2f66c7c1563057deb8602cdd6d295d1de)
para
y una función diferenciable cualquiera
en el intervalo unidad, para el caso en el que
se define
. Los polinomios de Bernoulli, como sus duales, forman un conjunto ortogonal de estados en el intervalo unidad, así se tiene:
![{\displaystyle \int _{0}^{1}{\tilde {B}}_{m}(x)B_{n}(x)\,dx=\delta _{mn}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0065bf1ba70c084d7e4333dcf2e8182d70816333)
y
![{\displaystyle \sum _{n=0}^{\infty }B_{n}(x){\tilde {B}}_{n}(y)=\delta (x-y).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e644b080dc8c8f64b96e920371b2b5d468474b8b)
La fórmula de Euler-MacLaurin se obtiene multiplicando la última igualdad por la función a sumar
e integrando el resultado sobre el intervalo unidad:
![{\displaystyle f(x)=\int _{0}^{1}\sum _{n=0}^{\infty }B_{n}(x){\tilde {B}}_{n}(y)f(y)\,dy}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95d5fe29713cd169c139c5ce7a52743b8f030c15)
![{\displaystyle =\int _{0}^{1}f(y)\,dy+\sum _{n=1}^{N}B_{n}(x){\frac {1}{n!}}\left[f^{(n-1)}(1)-f^{(n-1)}(0)\right]-{\frac {1}{(N+1)!}}\int _{0}^{1}B_{N+1}(x-y)f^{(N)}(y)\,dy.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/15b204247b009cba32ac1f0dfda2bf7521bff290)
Tomando
y reagrupando términos se obtiene la fórmula buscada junto con el término de error. Nótese que los números de Bernoulli se definen como
y que estos se anulan para n impares mayores que 1. Nótese también que en esta derivación se asume que la función
es suficientemente diferenciable, en particular,
ha de ser una función analítica
La fórmula de Euler-MacLaurin puede verse como una representación de funciones en el intervalo unidad por el producto directo de los polinomios de Bernoulli y sus duales. Sin embargo, esta representación no es completa en el conjunto de funciones cuadrado integrables. La expansión en término de polinomios de Bernoulli tiene una núcleo no trivial. En particular,
pertenece al núcleo, pues la integral de
se anula en el intervalo unidad, así como la diferencia de sus derivadas en los extremos del intervalo.
Referencias
Bibliografía
- Pierre Gaspard, "r-adic one-dimensional maps and the Euler summation formula", Journal of Physics A, 25 (letter) L483-L485 (1992). (Describe las autofunciones del operador de transferencia para el mapa de Bernoulli)
- Xavier Gourdon and Pascal Sebah, Introduction on Bernoulli's numbers, (2002)
- D.H. Lehmer, "On the Maxima and Minima of Bernoulli Polynomials", American Mathematical Monthly, volume 47, pages 533–538 (1940)
Datos: Q282023