DeducciónIntegraciónavanzado

Deducción: Simpson 1/3

De tres nodos equiespaciados a los pesos 1, 4, 1 de Simpson integrando el polinomio cuadrático de Lagrange.

Pesos de la parábola interpolante

  1. Tomamos x0=ax_0=a, x1=a+b2x_1=\frac{a+b}{2} y x2=bx_2=b. Escribimos h=b−a2h=\frac{b-a}{2} y usamos t=x−aht=\frac{x-a}{h}, de modo que tt recorre [0,2][0,2].

    dx=h dtdx=h\,dt

  2. Las bases cuadráticas en la variable t son:

    ℓ0(t)=(t−1)(t−2)2,ℓ1(t)=−t(t−2),ℓ2(t)=t(t−1)2\ell_0(t)=\frac{(t-1)(t-2)}{2},\quad \ell_1(t)=-t(t-2),\quad \ell_2(t)=\frac{t(t-1)}{2}

  3. Integramos cada base en [0,2] y multiplicamos por h:

    h∫02ℓ0(t)dt=h3,h∫02ℓ1(t)dt=4h3,h∫02ℓ2(t)dt=h3h\int_0^2\ell_0(t)dt=\frac{h}{3},\quad h\int_0^2\ell_1(t)dt=\frac{4h}{3},\quad h\int_0^2\ell_2(t)dt=\frac{h}{3}

  4. Al juntar los términos aparece Simpson 1/3:

    ∫abf(x) dx≈h3[f(a)+4f ⁣(a+b2)+f(b)]\int_a^b f(x)\,dx\approx\frac{h}{3}\left[f(a)+4f\!\left(\frac{a+b}{2}\right)+f(b)\right]