Ejercicio: trapecio y punto medio

Comparación de trapecio compuesto, punto medio simple y punto medio compuesto, con errores absolutos y relativos.

Comparación en una integral corta

EjemploIntegral de x^3 en [0,1]

Calcula I=01x3dxI=\int_0^1 x^3\,dx. Compara trapecio compuesto con n=2n=2, punto medio simple y punto medio compuesto con n=2n=2.

  1. El valor exacto es:

    I=01x3dx=[x44]01=14I=\int_0^1x^3\,dx=\left[\frac{x^4}{4}\right]_0^1=\frac14
  2. Para trapecio compuesto con n=2n=2, h=1/2h=1/2 y los nodos son 0,1/2,10,1/2,1:

    T2=14[f(0)+2f ⁣(12)+f(1)]=14[0+218+1]=516\begin{aligned}T_2&=\frac14\left[f(0)+2f\!\left(\frac12\right)+f(1)\right]\\&=\frac14\left[0+2\cdot\frac18+1\right]=\frac{5}{16}\end{aligned}
  3. Para punto medio simple, el centro del intervalo completo es m=1/2m=1/2:

    M1=(10)f ⁣(12)=18M_1=(1-0)f\!\left(\frac12\right)=\frac18
  4. Para punto medio compuesto con n=2n=2, los centros son 1/41/4 y 3/43/4:

    M2=12[f ⁣(14)+f ⁣(34)]=12(164+2764)=732\begin{aligned}M_2&=\frac12\left[f\!\left(\frac14\right)+f\!\left(\frac34\right)\right]\\&=\frac12\left(\frac{1}{64}+\frac{27}{64}\right)=\frac{7}{32}\end{aligned}
  5. Los errores absolutos y relativos son:

    meˊtodoIQIQ/IT21/1625%M11/850%M21/3212.5%\begin{array}{c|c|c} \text{método} & |I-Q| & |I-Q|/I \\ \hline T_2 & 1/16 & 25\% \\ M_1 & 1/8 & 50\% \\ M_2 & 1/32 & 12.5\% \end{array}

En este ejemplo, punto medio compuesto con dos subintervalos es el más preciso de los tres. La mejora viene de usar dos centros en lugar de un solo rectángulo para todo [0,1].

Usar solo puntos interiores

EjemploPunto medio con n=4 y n=8

Calcula I=∫_0^{π/2} sin(x)e^{-x} dx con punto medio compuesto usando n=4 y n=8.

  1. El valor exacto es el mismo del ejercicio anterior:

    I=0.396060211824619I=0.396060211824619
  2. Las aproximaciones de punto medio son:

    M4=0.409710138362011,M8=0.401008515076530M_4=0.409710138362011,\qquad M_8=0.401008515076530
  3. Sus errores absolutos son:

    E4=0.013649926537392,E8=0.004948303251911|E_4|=0.013649926537392,\qquad |E_8|=0.004948303251911

El error baja al aumentar el número de subintervalos, pero aquí sigue siendo mayor que Simpson para el mismo refinamiento.