s2Eva_IT2018_T2 Deducir Simpson 1/3

Ejercicio: 2Eva_IT2018_T2 Deducir Simpson 1/3

Para el problema, se usan los puntos: [a,f(a)], [b,f(b)] y [c,f(c)]
por donde pasa la curva f(x) aproximada a un polinomio de grado 2, $f(x) \approx p(x)$

\int_a^b f(x) dx \approx \int_a^b p(x) dx

p(x) = L_a f(a) + L_c f(c) + L_b f(b)

\int_a^b p(x) dx = \int_a^b \Big[ L_a f(a) +L_c f(c) + L_b f(b) \Big] dx =

= \int_a^b L_a f(a) dx + \int_a^b L_c f(c) dx + \int_a^b L_b f(b) dx

\int_a^b p(x) dx = I_1 + I_2 + I_3

Como referencia se usa la gráfica para relacionar a, b, c y h:

Primer Integral

Para el primer integral $I_1= \int_a^b L_a f(a) dx$ se tiene que:

L_a = \frac{(x-b)(x-c)}{(a-b)(a-c)} = \frac{(x-b)(x-c)}{(-2h)(-h)}

L_a = \frac{(x-b)(x-c)}{2h^2}

se convierte en:

I_1 = \int_a^b \frac{(x-b)(x-c)}{2h^2} f(a) dx = \frac{f(a)}{2h^2} \int_a^b (x-b)(x-c)dx

Para dejar la parte del integral en función de h, a y b, teniendo que c está en la mitad de [a,b], es decir c = (a+b)/2 , se usa:

u = x-a

por lo que $\frac{du}{dx}=1$ y $du = dx$

x-c = (u+a) - \frac{a+b}{2}

= u+ \frac{a-b}{2} = u - \frac{b-a}{2}

x-c = u-h

x-b = (u+a)-b = u-2\Big(\frac{b-a}{2}\Big) = u-2h

Se actualiza el integral de x entre [a,b] usando u = x-a, que se convierte el rango para u en [0, b-a] que es lo mismo que [0,2h]

\int_a^b (x-b)(x-c)dx = \int_0^{2h} (u-2h)(u-h)du =

= \int_0^{2h} \Big( u^2 - 2hu - uh + 2h^2 \Big) du

= \int_0^{2h} \Big( u^2 - 3hu + 2h^2 \Big) du

= \frac{u^3}{3}- 3h\frac{u^2}{2}+ 2h^2u \Big|_0^{2h}

= \frac{(2h)^3}{3}- 3h\frac{(2h)^2}{2} + 2h^2(2h) -(0-0+0)

= \frac{8h^3}{3}- 6h^3 + 4h^3 =\frac{8h^3}{3}- 2h^3

= \frac{2h^3}{3}

resultado que se usa en I₁

I_1= \frac{f(a)}{2h^2}\frac{2h^3}{3} =\frac{h}{3} f(a)

que es el primer término de la fórmula general de Simpson 1/3

Segundo Integral

Para el Segundo integral $I_2= \int_a^b L_c f(c) dx$ se tiene que:

L_c = \frac{(x-a)(x-b)}{(c-a)(c-b)} = \frac{(x-a)(x-b)}{(h)(-h)}

L_c = \frac{(x-a)(x-b)}{-h^2}

se convierte en:

I_2 = \frac{f(c)}{-h^2} \int_a^b (x-a)(x-b) dx

= \frac{f(c)}{-h^2} \int_0^{2h} (u)(u-2h) du

siendo:

\int_0^{2h}(u^2-2hu)du=\Big(\frac{u^3}{3}-2h\frac{u^2}{2}\Big)\Big|_0^{2h}

=\frac{(2h)^3}{3}-h(2h)^2-(0-0)

=\frac{8h^3}{3}-4h^3 = -\frac{4h^3}{3}

usando en I₂

I_2 = \frac{f(c)}{-h^2}\Big(-\frac{4h^3}{3}) = \frac{h}{3}4f(c)

Tarea: Continuar las operaciones para y tercer integral para llegar a la fórmula general de Simpson 1/3:

I = \frac{h}{3} \Big( f(a)+4f(c) + f(b) \Big)