Cuadratura gaussiana

En análisis numérico un método de cuadratura es una aproximación de una integral definida de una función. Una cuadratura de Gauss n, es una cuadratura construida para obtener el resultado exacto al integrar polinomios de grado 2n-1 o menos. Para esto selecciona los puntos de evaluación x_i y los pesos w_i de forma conveniente. La regla suele expresarse para una integral en el intervalo [−1, 1], y viene dada por la siguiente expresión:

En el caso particular en que ${displaystyle f(x)}$ es un polinomio de grado 2n-1 o menos, la cuadratura de Gauss da el valor exacto de la integral. En el caso general, tal cuadratura dará buenas aproximaciones si ${displaystyle f(x)}$ puede ser bien aproximada por un polinomio de grado 2n-1 o menos, en el intervalo [−1, 1].

También conocido como método de Gauss-Legendre, los coeficientes están dados por

donde ${displaystyle P_{n}}$ es el polinomios de Legendre de grado n en el intervalo [−1, 1], y los x_i son las raíces de dicho polinomio. La siguiente tabla muestra los valores de los x_i y los pesos asociados w_i, para distintos valores de n.

Aplicando un cambio de variable lineal, se puede convertir una integral en el intervalo [a, b], en una integral en [−1, 1]:

Luego se puede aplicar la Cuadratura de Gauss para aproximar la integral en [−1, 1]:

Aproxime la integral ${displaystyle f(x)=x^{3}+2x^{2}}$ de 1 a 5 cuando n = 2 mediante el método de cuadratura de Gauss y después comparelo con el resultado exacto.

${displaystyle n=2}$

${displaystyle 2n-1=2(2)-1=3}$

Con ${displaystyle n=2}$ podemos resolver la integral con exactitud para todos los polinomios de grado igual o menor a 3 para f(x)