Universidad de las Fuerzas Armadas

Métodos Numéricos

Actividad de Aprendizaje #3

Nombre: Mateo Perugachi NRC: 4141

Fecha: 07/08/2022

Ejercicio 5.1. Calcular usando la fórmula del punto medio:

𝑏
𝑎+𝑏
∫ 𝑓(𝑥) ≅ (𝑏 − 𝑎)𝑓 ( )
𝑎 2

La integral:

1
∫ 𝑥𝑑𝑥
0

Calcular la integral exacta y dar el error. Dibujar el área exacta y aproximada y explicar el
error.

Planteamiento:

Se utilizará la fórmula del punto medio dado en el índice del ejercicio, a partir de los límites
dados en la integral de la función y con los resultados obtenidos, se hará una comparación
con el valor real de la integral para analizar la exactitud de este método.

Procedimiento:

El valor de la integral real:

1
1
∫ 𝑥𝑑𝑥 =
0 2

Utilizando la fórmula del punto medio:

1
0+1
∫ 𝑥𝑑𝑥 ≅ (1 − 0)𝑓 ( )
0 2
 1
1
∫ 𝑥𝑑𝑥 ≅ (1) ( )
0 2

1
1
∫ 𝑥𝑑𝑥 ≅
0 2

Se halla el error:

𝐸𝑟𝑟𝑜𝑟 = 𝑉𝑎𝑙𝑜𝑟 𝑅𝑒𝑎𝑙 − 𝑉𝑎𝑙𝑜𝑟 𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜

1 1
𝐸𝑟𝑟𝑜𝑟 = − =0
2 2

La gráfica es:

Conclusión:

A partir de los resultados obtenidos, se puede observar que el error entre el valor real y el
obtenido a partir de la fórmula del punto medio es nula, lo que nos indica que este método es
ideal para realizar aproximaciones de integrales de primer grado.

Ejercicio 5.2. Usando la regla del trapecio compuesta con tres subintervalos, calcular:

1
∫ 𝑥 2 𝑑𝑥
0

Calcular la integral exacta y el error. Hacer un dibujo del área exacta y la aproximada.

Planteamiento:
Se obtendrá el valor de h, necesario para realizar la fórmula del trapecio compuesta, mismo
que se obtiene a partir de los intervalos dados en la integral. A partir del resultado, se hará
una comparación con el valor real tanto numérico y gráficamente.

Procedimiento:

Obtenemos el valor de h:

𝑏−𝑎
ℎ=
𝑛

1−0
ℎ=
3

1
ℎ=
3

Se obtienen loso valores de x:

𝑥0 = 𝑎

𝑥0 = 0

𝑥1 = 𝑎 + ℎ

1
𝑥1 =
3

𝑥2 = 𝑎 + 2 ∙ ℎ

2
𝑥2 =
3

Aplicamos la fórmula del trapecio compuesto:

𝑏 𝑛−1
𝑏−𝑎
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ ∙ [𝑓(𝑎) + 2 ∙ ∑ 𝑓 (𝑥𝑖 ) + 𝑓 (𝑏)]
𝑎 2𝑛
𝑖=1

𝑛=3

1
1−0
∫ 𝑥 2 𝑑𝑥 ≅ ∙ {𝑓 (0) + 2[𝑓 (𝑥1 ) + 𝑓 (𝑥2 )] + 𝑓(1)}
0 6
 1
1 1 4
∫ 𝑥 2 𝑑𝑥 ≅ ∙ [0 + 2 ( + ) + 1]
0 6 9 9

1
1 5
∫ 𝑥 2 𝑑𝑥 ≅ ∙ [2 ∙ ( ) + 1]
0 6 9

1
∫ 𝑥 2 𝑑𝑥 ≅ 0.3519
0

El valor real de la integral es:

1
∫ 𝑥 2 𝑑𝑥 = 0.3333
0

Obtenemos el valor del error:

|𝐸| = 𝑉𝑎𝑙𝑜𝑟 𝑟𝑒𝑎𝑙 − 𝑉𝑎𝑙𝑜𝑟 𝑎𝑝𝑟𝑜𝑥𝑖𝑚𝑎𝑑𝑜

|𝐸| = 0.3333 − 0.3519

|𝐸| = 0.0185

Calculamos el error del trapecio compuesto:

(𝑏 − 𝑎 ) ℎ 2
|𝐸1𝑇 | = ∙ 𝑓′′(𝑐)
12

1
(1 − 0) ∙
|𝐸1𝑇 | = 9 ∙2
12

|𝐸1𝑇 | = 0.0185

Obtenemos la gráfica del área en el programa Octave:

Conclusión:

A partir de los resultados obtenidos, se pudo analizar que es necesario tener en cuenta los
intervalos y el número de iteraciones que nos pida el problema, ya que esto es la base de la
resolución del método del trapecio compuesto. De igual forma, se puede observar que el error
es un poco alto, y esto se debe a que el número es bajo, entre mayor sea el número de
iteraciones mayor precisión en el resultado existirá.

Ejercicio 5.3. Usando la regla de Simpson calcular:

1
∫ 𝑥 3 𝑑𝑥
0

Calcular la integral exacta y calcular el error. ¿Qué curva utiliza para aproximar la función
la regla de Simpson?¿Cómo explicarías el error?

Planteamiento:

Se utilizará la fórmula de Simpson para obtener el valor de la integral dada en el problema,

a partir del resultado obtenido se realizará una comparación con el valor real y conocer el
error que existe entre los dos valores. De igual manera, se conocerá la función utilizada por
este método para obtener el resultado.

Procedimiento:

Aplicamos la fórmula de Simpson:

 𝑏
𝑏−𝑎 𝑎+𝑏
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ [𝑓 (𝑎) + 4𝑓 ( ) + 𝑓 (𝑏 ) ]
𝑎 6 2

1
3
1−0 3 0+1 3
∫ 𝑥 𝑑𝑥 ≅ [1 + 4 ( ) + 03 ]
0 6 2

1
1
∫ 𝑥 3 𝑑𝑥 ≅
0 4

El valor real es:

1
∫ 𝑥 3 𝑑𝑥 = 0.25
0

Calculamos el error:

|𝐸| = 𝑉𝑎𝑙𝑜𝑟 𝑟𝑒𝑎𝑙 − 𝑉𝑎𝑙𝑜𝑟 𝑎𝑝𝑟𝑜𝑥𝑖𝑚𝑎𝑑𝑜

|𝐸| = 0.25 − 0.25

|𝐸 | = 0

Conclusión:

Si lo ponemos dentro del plano gráfico, podemos mencionar que al igual que mi ejercicio 5.1
los errores se logran compensar, ahora bien, si queremos una forma más rigurosa hallamos
el respectivo error de la fórmula de Simpson:

𝑆
(𝑏 − 𝑎 ) 5 4
𝐸 =− 𝑓 (𝑐 )
2880

Al tener que 𝑓(𝑥 ) = 𝑥 3 , la cuarta derivada de dicha función es cero, por lo tanto, el error es
nulo. A partir de esto se puede decir que la fórmula de Simpson presenta una precisión de al
menos tres.

Ejercicio 5.4. Dada la integral:

∞
2
∫ 𝑒 −𝑥 𝑑𝑥 = √𝜋 ≈ 1.77245
−∞
Aproximar la integral de forma que el resultado sea exacto hasta las décimas usando la
fórmula de los trapecios compuesta ¿Cuántos subintervalos son necesarios? (Indicación: se
puede tomar un menor que (−∞, ∞), por ejemplo [-2,2]).

Procedimiento:

Utilizamos la aproximación:

∞ 2
2 2
∫ 𝑒 −𝑥 𝑑𝑥 ≈ ∫ 𝑒 −𝑥 𝑑𝑥
−∞ −2

Ahora bien, realizamos usando de manera repetida la fórmula simple:

Con n=1

2
2 2 − (−2) −(−2)2 2
∫ 𝑒 −𝑥 𝑑𝑥 ≈ [𝑒 + 𝑒 −2 ] = 0.073
−2 2

En este caso no funciona.

Con n=2

2 0 2
−𝑥 2 −𝑥 2 2
∫ 𝑒 𝑑𝑥 = ∫ 𝑒 𝑑𝑥 + ∫ 𝑒 −𝑥 𝑑𝑥
−2 −2 0

0 2
2 2 0 − (−2) −(−2)2 2 2 − 0 −22 2
∫ 𝑒 −𝑥 𝑑𝑥 + ∫ 𝑒 −𝑥 𝑑𝑥 ≈ [𝑒 + 𝑒 −0 ] + [𝑒 + 𝑒 −0 ] = 2.036
−2 0 2 2

De igual manera con este no funciona.

Con n=3

2 2
2 − 2
3 3
−𝑥 2 −𝑥 2 −𝑥 2 2
∫ 𝑒 𝑑𝑥 = ∫ 𝑒 𝑑𝑥 + ∫ 𝑒 𝑑𝑥 + ∫ 𝑒 −𝑥 𝑑𝑥
2 2
−2 −2 −
3 3
 2 2
− 2
3 3
−𝑥 2 −𝑥 2 2
∫ 𝑒 𝑑𝑥 + ∫ 𝑒 𝑑𝑥 + ∫ 𝑒 −𝑥 𝑑𝑥
2 2
−2 −
3 3

2 2 2
− 3 − (−2) 2 2
3 − (− 3) −(−2)2 2 2
−(−2)2 −(− ) −( )
≈ [𝑒 + 𝑒 3 ] + [𝑒 3 +𝑒 3 ]
2 2
2
2−3 2 2 2
+ [𝑒 −(2) + 𝑒 −(3) ] = 1.734
2

Con 𝑛 = 3, el resultado de la integral es el mismo que el dado por el problema.

Conclusión:

A partir de lo realizado en el ejercicio, se pudo analizar que si se escoge un intervalo

específico es posible obtener el resultado de una integral que tenga intervalos mayores. De
igual manera, se pudo concluir que entre mayor número de particiones se realice, mayor
exactitud existirá en el resultado respecto al valor real.

Ejercicio 5.5. Supongamos que la fórmula de la integración numérica

3
∫ 𝑓(𝑥 )𝑑𝑥 ≅ 𝐴0 𝑓(−1) + 𝐴1 𝑓(0) + 𝐴2 𝑓 (1)
−2

es de tipo interpolatorio. Determinar los coeficientes (pesos) de la fórmula utilizando las

condiciones del grado de precisión. ¿Cuál es el grado de precisión de dicha fórmula?
Razónese la respuesta. Calcular la integral si 𝑓 (𝑥 ) = cosh (𝑥). ¿Con qué curva crees que
hemos interpolado la función para aproximar su integral?

Planteamiento:

En este caso, se impondrá que la función presenta una precisión de 0 y 1, a partir de esto se
construirá un sistema de ecuaciones del cual se obtendrá los valores de 𝐴0 , 𝐴1 𝑦 𝐴2 , y con
ello se calculará el grado de precisión exacto de la formula dada en el problema.

Procedimiento:

Para que tenga precisión 0, además, la fórmula ha de ser exacta para f(x)=1 cumplir:
 2
∫ 1𝑑𝑥 = 𝐴0 (1) + 𝐴1 (1) + 𝐴2 (1)
−2

Para que tenga precisión 1, además, la fórmula ha de ser exacta para f(x)=x cumplir:

2
∫ 𝑥𝑑𝑥 = 𝐴0 (−1) + 𝐴1 (0) + 𝐴2 (1)
−2

Para que tenga precisión 2, además, la fórmula ha de ser exacta para 𝑓 (𝑥 ) = 𝑥 2 cumplir:

2
∫ 𝑥 2 𝑑𝑥 = 𝐴0 (−1)2 + 𝐴1 (0)2 + 𝐴2 (1)2
−2

Ahora tenemos un sistema lineal de 3 ecuaciones con 3 incógnitas:

𝐴0 + 𝐴1 + 𝐴2 = 4

−𝐴0 + 𝐴2 = 0

16
𝐴0 + 𝐴2 =
3
8 4
En nuestro caso la solución es: 𝐴0 = 𝐴2 = 3 𝑦 𝐴1 = − 3. Y por lo tanto la fórmula queda:

2
8 4 8
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ 𝑓(−1) − 𝑓 (0) + 𝑓(1)
−2 3 3 3

Por su construcción, sabemos que tiene precisión de al menos 2. Veamos cual es la

precisión para ver si tiene una precisión de al menos 3 y probamos con la 𝑓 (𝑥 ) = x3

2
8 4 8
∫ 𝑥 3 𝑑𝑥 = (−1)3 − (0)3 + (1)3
−2 3 3 3

Operamos:

8 8
0= − +0+
3 3

Y si tiene precisión al menos de 3

Vamos a ver la precisión 4 con la función de 𝑓 (𝑥 ) = 𝑥 4

2
8 4 8
∫ 𝑥 4 𝑑𝑥 = (−1)4 − (0)4 + 8(1)4
−2 3 3 3

Operamos y tenemos que:

64 8 8
≠ −0+
5 3 3

Por lo tanto, el grado de precisión de la fórmula es 3.

2
8 4 8
∫ cosh(𝑥 ) 𝑑𝑥 ≈ cosh(−1) − cosh(0) + cosh(1) = 6.89
−2 3 3 3

El resultado exacto es mayor, como se puede ver en el grafico realizado:

2
∫ cosh(𝑥 ) 𝑑𝑥 = sinh(𝑥 )2−2 = 7.25
−2

Y el error es E= 0.35

Conclusión:

Se pudo observar que la curva que interpola el problema dado, al usar tres puntos, es una
parábola. De igual manera, se pudo analizar que es necesario imponer grados de precisión a
una formula interpoladora para conocer los valores de sus coeficientes, y con los valores
obtenidos se puede obtener el grado de precisión exacto de la misma.

Ejercicio 5.6. Calcular los nodos y los pesos de la fórmula de Newton-Cotes abierta con dos
nodos para el intervalo [-1,1]. Determinar su grado de precisión. Calcular la integral si
𝑓 (𝑥 ) = cosh (𝑥). ¿Con qué curva crees que hemos interpolado la función para aproximar
su integral?

Planteamiento:

El procedimiento es el mismo que en el anterior problema, ya que es necesario imponer

grados de precisión al método dado, en este caso el de Newton-Cotes, y de esta manera,
obtener los valores de los nodos y los coeficientes de este. Obtenido dichos resultados, es
posible calcular el grado de precisión exacto.

Procedimiento:

1
1 1
∫ 𝑓(𝑥 )𝑑𝑥 ≈ 𝐴0 𝑓 (− ) + 𝐴1 𝑓 ( )
−1 3 3

Razonando planteamos el sistema:

1
𝑓 (𝑥 ) = 1 ∫ 1𝑑𝑥 = 𝐴0 (1) + 𝐴1 (1)
−1

1
1 1
𝑓 (𝑥 ) = 𝑥 ∫ 𝑥𝑑𝑥 = 𝐴0 (− ) + 𝐴1 ( )
−1 3 3

Y operando, tenemos:

𝐴0 + 𝐴1 = 2

1 1
− 𝐴0 + 𝐴1 = 0
3 3

Y resulta 𝐴0 = 𝐴1 = 1. Y para el intervalo [−1,1] la formula abierta de Newton-Cotes con

dos nodos es:

1
1 1
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ 𝑓 (− ) + 𝑓 ( )
−1 3 3

Para que tena precisión 2, además, la fórmula ha de ser exacta para 𝑓 (𝑥 ) = 𝑥 2 y cumplir:

1
1 2 1 2
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ (− ) + ( )
−1 3 3
Pero no es así, puesto que:

2 2
≠
3 9

Por lo tanto, el grado de precisión de la fórmula es 1.

2
1 1
∫ cosh(𝑥 ) 𝑑𝑥 ≈ cosh (− ) + cosh ( ) = 2.11
−2 3 3

El resultado exacto es:

1
∫ cosh(𝑥 ) 𝑑𝑥 = sinh(𝑥 )1−1 = 2.35
−1

Y nuestro error es:

𝐸 = 0.24

Conclusión:

A partir de los resultados obtenidos, se pudo analizar que el grado de precisión del método
de Newton-Cotes con nodos en el intervalo dado es de grado 1, sin embargo, para obtener
los valores de los coeficientes se supuso que su precisión era mayor. De igual manera, se
pudo observar que el error entre el valor real y calculado a partir de este método es
relativamente alto.

Ejercicio 5.7. 1. Obtener 𝑥0 para la fórmula de cuadratura

1
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ 2𝑓 (𝑥0 )
−1

tenga grado de precisión al menos 1. ¿Cuál es su grado de precisión?

2. Obtener 𝑥0 𝑦 𝑥1 para que la fórmula de cuadratura

1
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ 𝑓 (𝑥0 ) + 𝑓(𝑥1 )
−1

tenga grado de precisión al menos dos. ¿Cuál es su grado de precisión? ¿Es de Newton-
Cotes?
3. Usar la fórmula obtenida anteriormente para calcular un valor aproximado de

2
∫ (𝑥 2 + 1)𝑑𝑥
−2

Realizando previamente un cambio de variable adecuado. ¿Qué error se comete al aplicar

dicha fórmula en este ejercicio? Justificar la respuesta sin hacer la integral exacta.

Planteamiento:

Para poder resolver los 3 items dado en este problema, es necesario realizar un sistema de
ecuaciones para poder obtener los valores indicados para cada uno, imponiendo los grados
de precisión indicados. Y a partir de lo obtenido, se obtendrá el grado de exactitud real de la
fórmula obtenida.

Procedimiento:

1
1. Para que tena precisión al menos de uno, la fórmula de ∫−1 𝑓 (𝑥 )𝑑𝑥 ≈ 2𝑓 (𝑥0 ) ha de
ser exacta para 𝑓 (𝑥 ) = 1 𝑦 𝑓 (𝑥 ) = 𝑥:

1
𝑓 (𝑥 ) = 1 ∫ 1𝑑𝑥 = 2(1)
−1

1
𝑓 (𝑥 ) = 𝑥 ∫ 𝑥𝑑𝑥 = 2𝑥0
−1

Y tenemos las ecuaciones:

2=2

0 = 2𝑥0

Por lo tanto, podemos decir que 𝑥0 = 0 y la fórmula queda:

𝟏
∫ 𝒇(𝒙)𝒅𝒙 ≈ 𝟐𝒇(𝟎)
−𝟏

Que es la fórmula del punto medio para el intervalo de [-1,1]. Veamos si tiene precisión 2.
Tendría que ser exacta para 𝑓(𝑥 ) = 𝑥 2 . Pero:
 1
∫ 𝑥 2 𝑑𝑥 ≠ 2(0)2
−1

Por lo tanto, la fórmula tiene precisión 1.

2. Tendría que ser exacta para 𝑓(𝑥 ) = 1, 𝑓(𝑥 ) = 𝑥 𝑦 𝑓 (𝑥 ) = 𝑥 2 :

1
∫ 1𝑑𝑥 = 1 + 1
−1

1
∫ 𝑥𝑑𝑥 = 𝑥0 + 𝑥1
−1

1
∫ 𝑥 2 𝑑𝑥 = 𝑥02 + 𝑥12
−1

Es decir:

2=1+1

0 = 𝑥0 + 𝑥1

2
= 𝑥02 + 𝑥21
3
1 1
Resolviendo el sistema obtenemos que 𝑥0 = − 𝑦 𝑥1 = , con lo que la fórmula de
√3 √3

cuadratura es:

1
1 1
∫ 𝑓(𝑥 )𝑑𝑥 ≈ 𝑓 (− )+ 𝑓( )
−1 √3 √3

Veamos que si tiene precisión 3:

1 3 3
3
1 1
∫ 𝑥 𝑑𝑥 = (− ) +( )
−1 √3 √3

1
≠0
2

3. Vamos a hacer un cambio de variable lineal que pase del intervalo [−2,2] 𝑎 [−1,1].
Será de la forma 𝑥 = 𝑚𝑡 + 𝑛:
 𝑆𝑖 𝑥 = −2 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑡 = −1 − 2 = −𝑚 + 𝑛

𝑆𝑖 𝑥 = 2 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑡 = 1 2=𝑚+𝑛

Con lo que m=2 y n=0 y el cambio de variable es 𝑥 = 2𝑡 𝑦 𝑑𝑥 = 2𝑑𝑡:

2
𝟐 𝟏 2
2
1 1
∫ (𝒙 + 𝟏)𝒅𝒙 = ∫ ((𝟐𝒕 )2 + 1)2𝑑𝑡 ≈ 2 [(2 (− )) + 1] + [(2 ∗ ) + 1] = 9.33
−𝟐 −𝟏 √3 √3

Conclusión:

A partir de lo obtenido, en el literal 2, la fórmula calculada presenta un grado de precisión 2,

por lo que se puede decir que tiene mayor precisión que la regla de los trapecios, la cual tiene
una precisión de 1 con 2 nodos. Por otro lado, no se puede decir que es la fórmula de Newton-
Cotes, ya que, contando los extremos del intervalo, no tiene los nodos igualmente separados.
En el literal 3 la solución obtenida es exacta, y el error es cero, debido a que el integrando es
un polinomio de grado dos y el grado de precisión es dos.

Ejercicio 5.8. 1. Aproximar mediante las reglas del trapecio y de Simpson el valor de la
integral

3
𝐼 = ∫ (𝑥 3 + 1)𝑑𝑥
0

2. Comparar los valores aproximados con el valor exacto. ¿Se podría haber predicho alguno
de los errores?

3. Utilizando la regla de los trapecios compuesta para aproximar I, ¿qué número de

subintervalos será suficiente para que el error, en valor absoluto, sea menor que 10−6?

Planteamiento:

Se hará uso de la fórmula de los trapecios para el literal 1, para el literal 2 se usará la fórmula
de Simpson y la aplicamos de igual manera y a partir de esto se dirá que el valor exacto es el
dado por la fórmula de Simpson, y por último el error de la fórmula de los trapecioes
compuesta la aplicamos de manera inmediata.

Procedimiento:
 1. La fórmula de los trapecios es:

𝑏
𝑏−𝑎
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ [𝑓 (𝑎) + 𝑓 (𝑏)]
𝑎 2

La aplicamos:

3
3−0 3
∫ (𝑥 3 + 1)𝑑𝑥 ≈ [(0 + 1) + (33 + 1))] = 43.50
0 2

2. La fórmula de Simpson es:

𝑏
𝑏−𝑎 𝑎+𝑏
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ [𝑓(𝑎) + 4𝑓 ( ) + 𝑓 (𝑏 ) ]
𝑎 6 2

La aplicamos

𝑏
3−0 3
∫ 𝑓 (𝑥 )𝑑𝑥 ≈ [(0 + 1) + 4(1.53 + 1)) + (33 + 1))] = 23.25
𝑎 6

3. El error de la fórmula de los trapecios compuesta es:

(𝑏 − 𝑎)ℎ2 ′′
𝐸Tℎ = − 𝑓 (𝑐), 𝑐 ∈ (𝑎, 𝑏)
12

Como 𝑓(𝑥 ) = 𝑥 3 tenemos que 𝑓 ′′(𝑐 ) = 6𝑐 que para 𝑐 ∈ (0,3), como en ese intervalo la
derivada segunda es creciente, |𝑓 ′′(𝑐 )| < 6(3) = 18. Teniendo además en cuenta que ℎ =
𝑏−𝑎
donde n es el número de subintervalos y usando la fórmula del error:
𝑛

3−0 2
(3 − 0) (
|𝐸ℎT | = 𝑛 ) 18 < 10−6
12

Es decir:

3(3)2 18
2
< 10−6
12𝑛

Y despejando n:
 3(3)2 18(10)6
√ ≤𝑛
12

Conclusión:

En el literal 2, el valor exacto es el obtenido a partir del método de Simpson, ya que esta
presenta un grado de precisión 3 y se está realizando la integral de un polinomio de grado 3.

En el literal 3, se dice que 6363.96 ≤ 𝑛. Es decir, para garantizar, con este criterio, que el
error es menor que 10−6 , necesitamos al menos 6364 subintervalos. Esta es una condición
suficiente, no necesaria.

Ejercicio 5.9. Sea la integral

1
2
𝐼 = ∫ 𝑒 −𝑥 𝑑𝑥
0

1. Obtener el valor aproximado de la integral mediante la regla del trapecio compuesto

con dos subintervalos. Acotar el error en valor absoluto usando la fórmula del error.
2. Determinar el número n de subintervalos de modo que la regla del trapecio
compuesta aproxima el valor de I con un error menor que 10−6.

Planteamiento:

Para el literal 1, se utilizará la regla compuesta del trapecio para obtener el valor de la integral,
a partir de este resultado se comparará con el valor real y así obtener el error que existe entre
ambas y de igual forma se calculará el error dentro de la regla en sí. Para el literal 2, se
utilizará la formula del error, la cual será igualada al error dado por el problema y con esto
se obtendrá el valor de n.

Procedimiento:

1. Usamos la regla compuesta del trapecio:

b n−1
(b − a)
∫ f(x)dx = [f(a) + 2 ∑ f(xi ) + f(b)]
a 2n
i=1

2 1
Con n=2, f(x) = e−x , a = 0, b = 1 entonces x1 = 2 y:
 1
2 1 − 2 −02 2 2
∫ e−x dx ≈ [e + 2e−0.5 + e−1 ] = 0.7314
0 2(2)

Calculamos el valor exacto:

1
2
∫ e−x dx = 0.7468
0

Estimamos el error usando la formula del error de los trapecios compuesta:

(b − a)h2 ′′
EhT = − f (c), cϵ(a, b)
12

Calculamos h con n=2

b−a
h=
n

1
h=
2
2 2
Sabemos f(x) = e−x tenemos que f ′′ (c) = e−c (4c2 − 2).

Calculamos el error con la formula

12
(1 − 0)
|EhT | ≤ 2 2 = 0.042 < 0.1
12

Entonces sabemos que al menos una cifra decimal es exacta.

2. Usando la fórmula del error:

1−0 2
(1 − 0) (
n )
|EhT | = 2 < 10−6
12

(1)(1)2
2 ≤ 10−6
12 ∗ n2
 2(10)6
n≥√
12

n ≥ 408.2

Conclusión:

Para el literal 1, el valor exacto es 0.7468, el valor aproximado 0.7314 y efectivamente, una
cifra decimal es exacta. Y no hay dos cifras decimales exactas, por lo tanto, la estimación del
error fue muy buena.

Para el literal 2, para garantizar, con este criterio, que el error es menor que 10−6 es necesario
al menos 409 subintervalos. Esta es una condición suficiente, no necesaria.

Ejercicio 5.10. Es conocido que

2
1
𝑙𝑜𝑔2 = ∫ 𝑑𝑥
1 𝑥

Y damos como dato que si f tiene derivada cuarta continua en [a,b] entonces el error de a la
fórmula de Simpson compuesta viene dado por:

(𝑏 − a)h4 (iv)
𝐸ℎS =− 𝑓 (𝑐 ), 𝑐 ∈ (a, b)
180

1. Aproximar el valor de log 2 mediante la fórmula de Simpson compuesta que utiliza

dos subintervalos (es decir, cinco nodos) y acotar el error en valor absoluto usando
la fórmula del error. Comparar esta cota con el error verdadero.
2. Determinar el número total de nodos que serán suficientes para que la fórmula de
Simpson compuesta proporcione un valor aproximado de 𝑙𝑜𝑔2 con un error menor
que 10−4 .

Planteamiento:

Para el literal 1, es necesario conocer los valores de x a partir del valor de h, mismo que se
calcula con el intervalo dado en la integral. Con estos datos obtenidos, se calcula la integral
con la fórmula de Simpson compuesta y se compara el resultado con el valor real para conocer
el error que existe.

Para el literal 2, se iguala la ecuación del error de la fórmula de Simpson al error dado por el
literal y de esta manera se despejará n para así conocer su valor.

Procedimiento:

1. Dividimos [1, 2] en 2 subintervalos, y luego cada uno de ellos en otros 2 otra vez. Es
decir, en 4. Como la fórmula es de Newton-Cotes los nodos tienen que estar
igualmente separados y h = 1/4 = 0, 25:

𝑎 = 𝑥0 = 1, 𝑥1 = 1, 25𝑥2 = 1, 5𝑥3 = 1, 75 𝑦𝑥4 = 𝑏 = 2

Los nodos de los extremos, a y b suman con peso, o factor 1. A los interiores, si son pares,
en este caso x2, se les multiplica por 2 y si son impares, es decir x1 y x3 se suman y se
multiplican por 4. Es decir, el valor aproximado:

3 𝑚−1 𝑚
ℎ
∫ 𝑒 𝑥 𝑑𝑥 ≈ [𝑓(𝑎) + 2 · ∑ 𝑓 (𝑥2𝑖 ) + 4 · ∑ 𝑓(𝑥2𝑖 − 1) + 𝑓(𝑏)] =
0 3
𝑖=1 𝑖=1

1 2
0.25
= [𝑓(𝑎) + 2 · ∑ 𝑓 (𝑥2𝑖 ) + 4 · ∑ 𝑓 (𝑥2𝑖 − 1) + 𝑓(𝑏)] =
3
𝑖=1 𝑖=1

0.25 1 1 1 1 1 0.25 1 1 1 1 1
= [ + 2 + 4( + ) + ] = ·[ +2 + 4( + )+ ]
3 𝑎 𝑥2 𝑥1 𝑥3 𝑏 3 1 1.5 1.25 1.75 2
= 0.6933

El error de la fórmula de Simpson compuesta es:

(𝑏 − 𝑎)ℎ4 (𝑖𝑣)
𝐸ℎ𝑆 = − 𝑓 (𝑐 ), 𝑐 𝜖 (𝑎, 𝑏)
180

Como 𝑓 (𝑥) = 1/𝑥 tenemos que 𝑓 4) (𝑐) = 24/𝑐 5 que para c ∈ (1, 2), como en ese intervalo
la derivada cuarta es decreciente, |𝑓 4) (𝑐) | ≤ 24/1 = 24. Usando la fórmula del error:

0.254
|𝐸𝑛𝑆 | ≤ (2 − 1)5 · 24 = 0.0005 < 0.001
180
Y el error es menor que una milésima y tiene tres cifras decimales exactas.

2. El error de la fórmula de Simpson compuesta es:

(𝑏 − 𝑎) ℎ4 (𝑖𝑣)
𝐸ℎ𝑆 =− 𝑓 (𝑐 ), 𝑐 𝜖 (𝑎, 𝑏)
180
𝑏−𝑎
Como ℎ = 𝑛

(𝑏 − 𝑎)5 (𝑖𝑣)
𝐸ℎ𝑆 =− 𝑓 (𝑐), 𝑐 𝜖 (𝑎, 𝑏)
180𝑛4

1
|𝐸𝑛𝑆 | ≤ 24 < 0.0001
180𝑛4

Despejando n:

1
240000 4
( ) <𝑛
180

Conclusión:

Para el literal 1, el valor exacto es 0.6931, el valor aproximado 0.6933 y efectivamente, tres
cifras decimales son exactas. Y no hay cuatro cifras decimales exactas, por lo que, en este
caso, la estimación del error fue muy buena.

Para el literal 2, 6, 04 < n y si tomamos n = 8 (y como n = 2m) m = 4 subintervalos,

garantizamos que el error sea menor que 0, 0001. N = 7 no vale porque tiene que ser par.

Ejercicio 5.11. Si se utiliza la regla de los trapecios compuesta para aproximar:

2
𝐼 = ∫ ln2 (𝑥 ) 𝑑𝑥
1

¿Qué número de subintervalos será suficiente elegir para que el error fuese menor que
10−3 ?

Planteamiento:
Se utilizará la formula del error para el método de trapecios compuestos, y esta será iguala
con el error dado en el problema, como los demás datos ya están dados por el problema, es
posible despejar n y así conseguir su valor.

Procedimiento:

Utilizamos la fórmula de trapecios compuestos:

(𝑏 − 𝑎 )ℎ 2
𝐸ℎ𝑇 = − 𝑓"(𝑐)
12

Como 𝑓(𝑥 ) = ln2 (𝑥 ):

2 ln 𝑐 − 2
𝑓"(𝑐) =
𝑐2

Como 𝑓 ′′(1) = 2 y 𝑓 ′′(2) = 0.15, el máximo en valor absoluto es 2 y |𝑓 ′ ′(𝑐) ≤ 2.

Usamos la fórmula del error:

2−1 2
(2 − 1) (
|𝐸ℎ𝑇 | ≤ 𝑛 ) 2 < 0.001
12

2
|𝐸ℎ𝑇 | ≤ < 0.001
12𝑛2

Despejamos n:

2000
√ <𝑛
12

Como sabemos que n=13

12.9 < 13

Conclusión:
Se dice que, 12.9 < 𝑛 y se necesita 𝑛 = 13 subintervalos para garantizar un error menor que
una milésima. (De hecho, el valor exacto es 0.1883 y el valor aproximado con 𝑛 = 13 es
0.1887 y vemos que aproximamos de forma exacta tres cifras decimales).

Ejercicio 5.12. 1. Encontrar las constantes 𝐴0 , 𝐴1 , x1 de modo que la fórmula de cuadratura

simple:

1
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ A0 𝑓 (0) + A1 𝑓(𝑥1 )
0

Tenga el grado de precisión más alto posible. ¿Cuál es el grado de precisión?¿Es una
fórmula de Newton-Cotes?

2.- Realizar un cambio de variable en la integral

2
2
∫ (1 + e−x )dx
0

Para poder utilizar la fórmula obtenida anteriormente y calcular de esa manera una
aproximación del valor de dicha integral.

Planteamiento:

Para el literal 1, se impondrá que el grado de precisión de la fórmula es de 3, a partir de esto,

se construirá un sistema de ecuaciones para poder obtener los valores de las variables
impuestas por el problema, se armará la fórmula final y se calculará su grado de precisión
exacto.

Para el literal 2, se utilizará un cambio de variable, de manera que 𝑥 = mt + n, y con esto se

podrá calcular la integral a partir de la fórmula obtenida en el literal anterior.

Procedimiento:

Imponemos tres condiciones; para que la fórmula sea exacta f(x)=1; f(x)=x y f(x)=𝑥 2 .

1
∫ 1 𝑑𝑥 = 𝐴0 (1) + 𝐴1 (1)
0
 1
∫ 𝑥 𝑑𝑥 = 𝐴0 (0) + 𝐴1 (𝑥0 )
0

1
∫ 𝑥 2 𝑑𝑥 = 𝐴0 (02 ) + 𝐴1 (𝑥02 )
0

Es decir:

1 = 𝐴0 + 𝐴1

1
= 𝐴1 𝑥0
2

1
= 𝐴1 (𝑥02 )
3

Los valores de A son:

1 3
𝐴0 = , 𝐴1 =
4 4

2
𝑥0 =
3

Reemplazamos en la fórmula:

1
1 3 2
∫ 𝑓 (𝑥 )𝑑𝑥 = 𝑓 (0) + 𝑓 ( )
0 4 4 3

No tiene nodos igualmente separados, por ello investigamos f(x)=𝑥 3 :

1
1 3 2 3
∫ 𝑥 3 𝑑𝑥 = (0)3 + ( )
0 4 4 3

1 3 2 3
≠ ( )
4 4 3

1 2
≠
4 9
El grado de precisión de la fórmula es dos, teniendo en cuenta que la fórmula de los trapecios,
con dos nodos tienen precisión uno.

Literal 2:

Para encontrar la integral utilizamos la fórmula obtenida:

𝑥 = 𝑚𝑡 + 𝑛

𝑆𝑖 𝑥 = 0, 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑡 = 0; 0 = 𝑛

𝑆𝑖 𝑥 = 2, 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑡 = 1; 2 = 𝑚

Como sabemos m=2 y n=0, el cambio de la variable es x=2t y dx=2dt:

2 1
2 2
∫ (1 + 𝑒 −𝑥 )𝑑𝑥 = ∫ (1 + 𝑒 −(2𝑡) )2𝑑𝑡
0 0

2
2 2
−𝑥 2 )𝑑𝑥
1 −(2(0))
2 3 −(2( ))
3
∫ (1 +𝑒 = 2[ (1 + 𝑒 + (1 + 𝑒 ]
0 4 4

2
2
∫ (1 + 𝑒 −𝑥 )𝑑𝑥 = 2.7535
0

4
La solución exacta es 2.8821. Los nodos para intervalos [0,2] son 𝑡0 = 0 𝑦 𝑡1 = 3.

Conclusión:

Para el literal 1, el grado de precisión de la fórmula es dos, que no está mal, teniendo en
cuenta que la fórmula de los trapecios, con dos nodos, tiene precisión uno.

Para el literal 2, la solución obtenida muestra un error no muy alto con respecto al valor real,
por lo que se puede decir que este método es útil para realizar la integral de una ecuación, sin
embargo, es necesario aumentar el número de subintervalos para obtener mayor precisión en
los resultados.