Método de Romberg
# Método de Romberg

En análisis numérico, el Método de Romberg genera una matriz triangular cuyos elementos son estimaciones numéricas de la integral definida siguiente:

    $\int _{a}^{b}f(x)\,dx$

usando la extrapolación de Richardson de forma reiterada en la regla del trapecio. El método de Romberg evalúa el integrando en puntos equiespaciados del intervalo de integración estudiado. Para que este método funcione, el integrando debe ser suficientemente derivable en el intervalo, aunque se obtienen resultados bastante buenos incluso para integrandos poco derivables. Aunque es posible evaluar el integrando en puntos no equiespaciados, en ese caso otros métodos como la cuadratura gaussiana o la cuadratura de Clenshaw–Curtis son más adecuados.

## Método

El método se define de forma recursiva así:

    $R(0,0)={\frac {1}{2}}(b-a)(f(a)+f(b))$
    $R(n,0)={\frac {1}{2}}R(n-1,0)+h_{n}\sum _{k=1}^{2^{n-1}}f(a+(2k-1)h_{n})$
    $R(n,m)=R(n,m-1)+{\frac {1}{4^{m}-1}}(R(n,m-1)-R(n-1,m-1))$

o

    $R(n,m)={\frac {1}{4^{m}-1}}(4^{m}R(n,m-1)-R(n-1,m-1))$

donde

    $n\geq 1$
    $m\geq 1$
    $h_{n}={\frac {b-a}{2^{n}}}.$

La cota superior asintótica del error de R(n,m) es:

    $O\left(h_{n}^{2^{m+1}}\right).$

La extrapolación a orden cero $R(n,0)$ es equivalente a la Regla del trapecio con $n+2$ puntos. a orden uno $R(n,1)$ es equivalente a la Regla de Simpson con $n+2$ puntos.

Cuando la evaluación del integrando es numéricamente costosa, es preferible reemplazar la interpolación polinómica de Richardson por la interpolación racional propuesta por Bulirsch & Stoer.

## Implementación en Python del método de Romberg

La siguiente es una implementación del método de Romberg en Python:

## Ejemplo

Como ejemplo, se le integra la función gaussiana en el intervalo $[O,1]$, esto es la función error evaluada en 1, cuyo valor es ${\rm {erf}}(1)\doteq 0.842700792949715$. Se calculan los elementos de la matriz triangular fila a fila, terminando los cálculos cuando la diferencia entre las dos últimas filas es menor que 1E-8.

El resultado en la esquina inferior derecha de la matriz triangular es el resultado correcto con la precisión deseada. Nótese que este resultado se deriva de aproximaciones mucho peores obtenidas con la regla del trapecio mostradas aquí en la primera columna de la matriz triangular.

