0% encontró este documento útil (0 votos)
32 vistas14 páginas

Integracion - Numerica

Este documento describe tres métodos de integración numérica: el método del trapecio, el método del trapecio de aplicación múltiple y el método de Simpson. Explica cómo calcular la integral aproximada y el error de cada método para diferentes funciones.

Cargado por

201912592d
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
32 vistas14 páginas

Integracion - Numerica

Este documento describe tres métodos de integración numérica: el método del trapecio, el método del trapecio de aplicación múltiple y el método de Simpson. Explica cómo calcular la integral aproximada y el error de cada método para diferentes funciones.

Cargado por

201912592d
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

INTEGRACION NUMERICA

Las formulas de Newton Cotes son los tipos de integración


mas comunes. Se basan en la estrategia de reemplazar
una función complicada o datos tabulados por un
polinomio de aproximación que es fácil de integrar:
𝑏 𝑏
𝐼 = ∫ 𝑓(𝑥) 𝑑𝑥 ≡ ∫ 𝑓𝑛 (𝑥) 𝑑𝑥
𝑎 𝑎
1 Método del Trapecio o la Regla del Trapecio
El método del trapecio consiste en emplear la función 𝑓𝑛 (𝑥)
de primer grado, es decir 𝑓1 (𝑥) es de primer grado
𝑏 𝑏
𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 ≡ ∫ 𝑓1 (𝑥) 𝑑𝑥
𝑎 𝑎
Sabemos que una recta se puede representar como:
𝑓(𝑏) − 𝑓(𝑎)
𝑓1 (𝑥) = 𝑓(𝑎) + (𝑥 − 𝑎)
𝑏−𝑎
𝑓(𝑏) − 𝑓(𝑎) 𝑎 𝑓(𝑏) − 𝑎 𝑓(𝑎)
𝑓1 (𝑥) = 𝑥 + 𝑓(𝑎) −
𝑏−𝑎 𝑏−𝑎
𝑓(𝑏) − 𝑓(𝑎) 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑎) − 𝑎 𝑓(𝑏) + 𝑎 𝑓(𝑎)
𝑓1 (𝑥) = 𝑥+
𝑏−𝑎 𝑎−𝑏
𝑓(𝑏) − 𝑓(𝑎) 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑏)
𝑓1 (𝑥) = 𝑥+
𝑏−𝑎 𝑏−𝑎
𝑏
𝑓(𝑏) − 𝑓(𝑎) 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑏)
𝐼=∫ [ 𝑥+ ] 𝑑𝑥
𝑎 𝑏 − 𝑎 𝑏−𝑎
𝑏
𝑓(𝑏) − 𝑓(𝑎) 𝑥 2 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑏)
𝐼= + 𝑥|
𝑏−𝑎 2 𝑏−𝑎 𝑎
𝑓(𝑏) − 𝑓(𝑎) (𝑏2 − 𝑎2 ) 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑏)
𝐼= + (𝑏 − 𝑎)
𝑏−𝑎 2 𝑏−𝑎

pág. 1
𝑎+𝑏
𝐼 = [𝑓(𝑏) − 𝑓(𝑎)] + 𝑏 𝑓(𝑎) − 𝑎 𝑓(𝑏)
2
𝑎𝑓(𝑏) + 𝑏𝑓(𝑏) − 𝑎𝑓(𝑎) − 𝑏𝑓(𝑎) + 2𝑏𝑓(𝑎) − 2𝑎𝑓(𝑏)
𝐼=
2
𝑏𝑓(𝑏) − 𝑎𝑓(𝑏) + 𝑏𝑓(𝑎) − 𝑎𝑓(𝑎)
𝐼=
2
(𝑏 − 𝑎)𝑓(𝑏) + (𝑏 − 𝑎)𝑓(𝑎)
𝐼=
2
𝑓(𝑎) + 𝑓(𝑏)
𝐼 = (𝑏 − 𝑎)
2
Error del método del trapecio es:
Para determinar el error del método del trapecio se utiliza
la expresión:
1
𝐸𝑎 = − 𝑓´´(𝜉)(𝑏 − 𝑎)3
12
Donde 𝑓´´(𝜉) se calcula con la expresión
𝑏
∫ 𝑓(𝑥) 𝑑𝑥
𝑓´´(𝜉) = 𝑀𝑒𝑑𝑖𝑎(𝑓(𝑥)) = 𝑎
𝑏−𝑎

Ejemplo
Hallar la integral de la función 𝑓(𝑥) = 𝑠𝑒𝑛(𝑥 2 ) desde 0 hasta
1.6 por el método del trapecio y su error verdadero y
aproximado si su verdadero valor es 0.845269, con seis
decimales.
Solución
𝑎=0 𝑓(𝑎) = 0
𝑏 = 1.6 𝑓(𝑏) = 0.549355
𝑓(𝑎) + 𝑓(𝑏)
𝐼 = (𝑏 − 𝑎)
2

pág. 2
0 + 0.549355
𝐼 = (1.6) = 0.439484
2
Error verdadero
𝐸𝑡 = 0.845269 − 0.439484 = 0.405785
0.845269 − 0.439484
𝐸𝑡 = | | × 100 = 48.01%
0.845269
Error aproximado es:
1
𝐸𝑎 = − 𝑓´´(𝜉)(𝑏 − 𝑎)3
12
𝑏
∫ 𝑓´´(𝑥)
𝑓´´(𝜉) = 𝑓´´(𝑥) = 𝑎
𝑏−𝑎
𝑓´(𝑥) = 2𝑥 cos(𝑥 2 ) 𝑓´´(𝑥) = 2cos(𝑥 2 ) − 4𝑥 2 𝑠𝑒𝑛 (𝑥 2 )
1.6
∫ [2cos(𝑥 2 ) − 4𝑥 2 𝑠𝑒𝑛 (𝑥 2 )] 𝑑𝑥 = −2.676876
0
−2.676876
𝑓´´(𝜉) = = −1.671172
1.6
1
𝐸𝑎 = − (−1.671172)(1.6)3 = 0.5704227
12

2 Método del Trapecio de Aplicación Múltiple


Este método consiste en dividir el intervalo dado en
intervalos mas pequeños igualmente espaciados
𝑏−𝑎
ℎ=
𝑛
luego los puntos que determina cada intervalo igualmente
espaciados se determina como:
𝑥𝑖 = 𝑎 + ℎ𝑖, 𝑖 = 1, . . . , 𝑛 − 1
La integral esta dado por:
𝑥1 𝑥2 𝑏
𝐼 = ∫ 𝑓1 (𝑥) 𝑑𝑥 + ∫ 𝑓1 (𝑥) 𝑑𝑥 + ⋯ + ∫ 𝑓1 (𝑥) 𝑑𝑥
𝑎 𝑥1 𝑥𝑛−1

pág. 3
𝑓(𝑎) + 𝑓(𝑥1 ) 𝑓(𝑥1 ) + 𝑓(𝑥2 )
𝐼 = (𝑥1 − 𝑎) + (𝑥2 − 𝑥1 ) +⋯
2 2
𝑓(𝑥𝑛−1 ) + 𝑓(𝑏)
+ (𝑏 − 𝑥𝑛−1 )
2
𝑓(𝑎) + 𝑓(𝑥1 ) 𝑓(𝑥1 ) + 𝑓(𝑥2 ) 𝑓(𝑥𝑛−1 ) + 𝑓(𝑏)
𝐼=ℎ +ℎ + ⋯+ ℎ
2 2 2
𝑛−1

𝐼 = [𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏)]
2
𝑖=1
𝑛−1
𝑛ℎ
𝐼= [𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏)]
2𝑛
𝑖=1
𝑓(𝑎) + 2 ∑𝑛−1
𝑖=1 𝑓(𝑥𝑖 ) + 𝑓(𝑏)
𝐼 = (𝑏 − 𝑎)
2𝑛
El error de este método se calcula por:
(𝑏 − 𝑎)3
𝐸𝑎 = − 𝑓´´(𝜉)
12𝑛2

Ejemplo
Hallar la integral de la función 𝑓(𝑥) = 𝑥 3 − 6𝑥 2 + 11𝑥 − 6
desde 1 hasta 2 por el método del trapecio múltiple con
𝑛 = 4 y su error aproximado con seis decimales.
Solución
2−1
ℎ= = 0.25 → 𝑥1 = 1.25 𝑥2 = 1.5 𝑥3 = 1.75
4
𝑎 = 1 𝑓(𝑎) = 0
𝑥1 = 1.25 𝑓(𝑥1 ) = 0.328125
𝑥2 = 1.5 𝑓(𝑥2 ) = 0.375000
𝑥3 = 1.75 𝑓(𝑥3 ) = 0.234375
𝑏=2 𝑓(𝑏) = 0

pág. 4
𝑓(𝑎) + 2(𝑓(𝑥1 ) + 𝑓(𝑥2 ) + 𝑓(𝑥3 )) + 𝑓(𝑏)
𝐼 = (𝑏 − 𝑎)
2𝑛
0 + 2(0.328125 + 0.375000 + 0.234375) + 0
𝐼=1
8
𝐼 = 0.234375
El error aproximado es:
(𝑏 − 𝑎)3
𝐸𝑎 = − 𝑓´´(𝜉)
12𝑛2
𝑏
∫ 𝑓´´(𝑥)
𝑓´´(𝜉) = 𝑎
𝑏−𝑎
𝑓´(𝑥) = 3𝑥 2 − 12𝑥 + 11 𝑓´´(𝑥) = 6𝑥 − 12
2
∫ (6𝑥 − 12) 𝑑𝑥 = −3
1
−3
𝑓´´(𝜉) = = −3
1
13
𝐸𝑎 = − (−3) = 0.015625
12(16)

3 Método de Simpson o Regla de Simpson 1/3


La regla de Simpson 1/3 resulta cuando un polinomio de
interpelación de segundo grado.
𝑏 𝑏
𝐼 = ∫ 𝑓(𝑥) 𝑑𝑥 ≡ ∫ 𝑓2 (𝑥) 𝑑𝑥
𝑎 𝑎
Para determinar 𝑓2 (𝑥) utilizamos los polinomios de
interpolación de Lagrange, la cual es:
𝑛

𝑓𝑛 (𝑥) = ∑ 𝐿𝑖 (𝑥)𝑓(𝑥𝑖 )
𝑖=0
Donde:

pág. 5
𝑛
𝑥 − 𝑥𝑗
𝐿𝑖 (𝑥) = ∏
𝑥𝑖 − 𝑥𝑗
𝑗=0
𝑗≠𝑖

Para 𝑛 = 2 se tiene
(𝑥 − 𝑥1 )(𝑥 − 𝑏) (𝑥 − 𝑎)(𝑥 − 𝑏)
𝑓2 (𝑥) = 𝑓(𝑎) + 𝑓(𝑥1 )
(𝑎 − 𝑥1 )(𝑎 − 𝑏) (𝑥1 − 𝑎)(𝑥1 − 𝑏)
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
+ 𝑓(𝑏)
(𝑏 − 𝑎)(𝑏 − 𝑥1 )
𝑏 (𝑥 − 𝑥1 )(𝑥 − 𝑏) (𝑥 − 𝑎)(𝑥 − 𝑏)
𝐼=∫ [ 𝑓(𝑎) + 𝑓(𝑥1 )
𝑎 (𝑎 − 𝑥1 )(𝑎 − 𝑏) (𝑥1 − 𝑎)(𝑥1 − 𝑏)
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
+ 𝑓(𝑏)] 𝑑𝑥
(𝑏 − 𝑎)(𝑏 − 𝑥1 )
Después de integrar y simplificar se obtiene:

𝐼 ≡ [𝑓(𝑎) + 4𝑓(𝑥1 ) + 𝑓(𝑏)]
3
𝑏−𝑎 𝑓(𝑎) + 4𝑓(𝑥1 ) + 𝑓(𝑏)
ℎ= → 𝐼 ≡ (𝑏 − 𝑎)
2 6
el error de este método está dado por:
(𝑏 − 𝑎)5 (4)
𝐸𝑎 = − 𝑓 (𝜉)
2880

4 Método de Simpson 1/3 de aplicación Múltiple


Así como en la regla del trapecio, la regla de Simpson se
mejora al dividir el intervalo de integración en varios
segmentos de un mismo tamaño
𝑏−𝑎
ℎ=
𝑛
la integral se puede expresar como:

pág. 6
𝑥1 𝑥2 𝑏
𝐼 = ∫ 𝑓2 (𝑥) 𝑑𝑥 + ∫ 𝑓2 (𝑥) 𝑑𝑥 + ⋯ + ∫ 𝑓2 (𝑥) 𝑑𝑥
𝑎 𝑥1 𝑥𝑛−1
ℎ ℎ
𝐼≡ [𝑓(𝑎) + 4𝑓(𝑥1 ) + 𝑓(𝑥2 )] + [𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + 𝑓(𝑥4 )] + ⋯
3 3

+ [𝑓(𝑥𝑛−2 ) + 4𝑓(𝑥𝑛−1 ) + 𝑓(𝑏)]
3
Efectuando y reduciendo se tiene:
𝑛−1 𝑛−2

𝐼 ≡ [𝑓(𝑎) + 4 ∑ 𝑓(𝑥𝑖 ) + 2 ∑ 𝑓(𝑥𝑗 ) + 𝑓(𝑏)]
3
1=1,3,5 𝑗=2,4,
𝑛−1 𝑛−2
(𝑏 − 𝑎)
𝐼≡ [𝑓(𝑎) + 4 ∑ 𝑓(𝑥𝑖 ) + 2 ∑ 𝑓(𝑥𝑗 ) + 𝑓(𝑏)]
3𝑛
1=1,3,5 𝑗=2,4,

El error de este método está dado por:


(𝑏 − 𝑎)5 (4)
𝐸𝑎 = − 𝑓 (𝜉)
180𝑛4
Donde 𝑓 (4) (𝜉) es el promedio de la cuarta derivada en el
intervalo

5 Método de Simpson 3/8


De manera similar a la obtención de la regla del trapecio
y Simpson 1/3 es posible ajustar un polinomio de
Lagrange de tercer grado como:
𝑏 𝑏
𝐼 ≡ ∫ 𝑓(𝑥) 𝑑𝑥 ≡ ∫ 𝑓3 (𝑥) 𝑑𝑥
𝑎 𝑎
Para determinar 𝑓3 (𝑥) utilizamos los polinomios de
interpolación de Lagrange, la cual es:
𝑛

𝑓𝑛 (𝑥) = ∑ 𝐿𝑖 (𝑥)𝑓(𝑥𝑖 )
𝑖=0

pág. 7
Donde:
𝑛
𝑥 − 𝑥𝑗
𝐿𝑖 (𝑥) = ∏
𝑥𝑖 − 𝑥𝑗
𝑗=0
𝑗≠𝑖

Para 𝑛 = 3 se tiene
(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑏)
𝑓3 (𝑥) = 𝑓(𝑎)
(𝑎 − 𝑥1 )(𝑎 − 𝑥2 )(𝑎 − 𝑏)
(𝑥 − 𝑎)(𝑥 − 𝑥1 )(𝑥 − 𝑏)
+ 𝑓(𝑥1 )
(𝑥1 − 𝑎)(𝑥1 − 𝑥2 )(𝑥1 − 𝑏)
(𝑥 − 𝑎)(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
+ 𝑓(𝑥2 )
(𝑥2 − 𝑎)(𝑥2 − 𝑥1 )(𝑥1 − 𝑏)
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
+ 𝑓(𝑏)
(𝑏 − 𝑎)(𝑏 − 𝑥1 )(𝑏 − 𝑥2 )
𝑏 (𝑥
− 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑏)
𝐼≡∫ [ 𝑓(𝑎)
𝑎 (𝑎 − 𝑥1 )(𝑎 − 𝑥2 )(𝑎 − 𝑏)
(𝑥 − 𝑎)(𝑥 − 𝑥1 )(𝑥 − 𝑏)
+ 𝑓(𝑥1 )
(𝑥1 − 𝑎)(𝑥1 − 𝑥2 )(𝑥1 − 𝑏)
(𝑥 − 𝑎)(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
+ 𝑓(𝑥2 )
(𝑥2 − 𝑎)(𝑥2 − 𝑥1 )(𝑥1 − 𝑏)
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )
+ 𝑓(𝑏)] 𝑑𝑥
(𝑏 − 𝑎)(𝑏 − 𝑥1 )(𝑏 − 𝑥2 )
Integrando se tiene:
3ℎ
𝐼≡ [𝑓(𝑎) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓(𝑏)]
8
(𝑏 − 𝑎)
𝐼≡ [𝑓(𝑎) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓(𝑏)]
8
3
𝐸𝑎 = − ℎ5 𝑓 (4) (𝜉)
80
𝑏−𝑎 (𝑏 − 𝑎)5 (4)
ℎ= → 𝐸𝑎 = − 𝑓 (𝜉)
3 6480

pág. 8
Ejemplo
Integrar la función 𝑓(𝑥) = 𝑥 5 − 5𝑥 4 − 15𝑥 3 + 85𝑥 2 − 26𝑥 − 120
Por el método de Simpson 3/8 desde 2.1 hasta 2.9, si su
vedadero valor es I=8.899053646
Solución
0.8
Se requiere 4 puntos equidistantes ℎ = = 0.2667
3
𝑎 = 2.1 → 𝑓(2.1) = 4.935510
𝑥1 = 2.3667 → 𝑓(2.3667) = 13.108094
𝑥2 = 2.6334 → 𝑓(2.6334) = 13.244781
𝑏 = 2.9 → 𝑓(2.9) = 2.085990
Por lo tanto la integra será:
3ℎ
𝐼≡ [𝑓(𝑎) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓(𝑏)]
8
3(0.2667)
𝐼≡ [4.935510 + 3(13.108094) + 3(13.244781)
8
+ 2.085990]
𝐼 = 8.609088
El error verdadero es:
8.899053646−8.609088
𝐸𝑡 = | | 100 =3.25 %
8.899053646
Para hallar su error aproximado
3
𝐸𝑎 = − ℎ5 𝑓 (4) (𝜉)
80
Como 𝑓(𝑥) = 𝑥 − 5𝑥 − 15𝑥 3 + 85𝑥 2 − 26𝑥 − 120
5 4

𝑓´(𝑥) = 5𝑥 4 − 20𝑥 3 − 45𝑥 2 + 170𝑥 − 26


𝑓´´(𝑥) = 20𝑥 3 − 60𝑥 2 − 90𝑥 + 170
𝑓´´´(𝑥) = 60𝑥 2 − 120𝑥 − 90
𝑓 (4) (𝑥) = 120𝑥 − 120

pág. 9
1 2.9 1
𝑓 (4) (𝜉) = ∫ (120𝑥 − 120) 𝑑𝑥 = [60𝑥 2 − 120𝑥]2.9
2.1
0.8 2.1 0.8
387.5 − 12.6
𝑓 (4) (𝜉) = = 468.625
0.8
3
𝐸𝑎 = − (0.8)5 (468.625) = −5.758464
80

Ejemplo
Evalúe la integral siguiente:
4
∫ (1 − 𝑒 −2𝑥 )𝑑𝑥
0
a) En forma analítica
b) con la regla de trapecio
c) con la regla de trapecio múltiple, con 𝑛 = 2 𝑦 4
d) con la regla de Simpson 1/3
e) con la regla de Simpson múltiple 1/3, con 𝑛 = 4
Solución
4
4 𝑒 −2𝑥
a) 𝐼 = ∫0 (1 − 𝑒 −2𝑥 )𝑑𝑥 = (𝑥 + |
2 0

𝑒 −8 1
𝐼 =4+ − = 3.500168 𝑣𝑎𝑙𝑜𝑟 𝑣𝑒𝑟𝑑𝑎𝑑𝑒𝑟𝑜
2 2

𝑓(𝑎)+𝑓(𝑏) 𝑓(0)+𝑓(4)
b) 𝐼 = (𝑏 − 𝑎) = 4( )
2 2
0+ (1 − 𝑒 −8 )
𝐼 = 4( ) = 1.999329
2
3.500168 − 1.999329
𝑒𝑟 = | | ∗ 100 = 42.88 %
3.500168

c) Regla del trapecio múltiple n=2 y 4

pág. 10
𝑛−1
𝑏−𝑎
𝐼= (𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏))
2𝑛
𝑖=1
n=2
4
𝐼 = (𝑓(0) + 2𝑓(2) + 𝑓(4))
4
𝐼 = (0 + 2(1 − 𝑒 −4 ) + (1 − 𝑒 −8 )) = 2.728994
3.500168 − 2.728994
𝑒𝑟 = | | ∗ 100 = 22.03 %
3.500168
n=4
4
𝐼 = (𝑓(0) + 2𝑓(1) + 2𝑓(2) + 2𝑓(3) + 𝑓(4))
8
1
𝐼 = (0 + 2(1 − 𝑒 −2 ) + 2(1 − 𝑒 −4 ) + 2(1 − 𝑒 −6 ) + (1 − 𝑒 −8 ))
2
1
𝐼 = (7 − 2𝑒 −2 − 2𝑒 −4 − 2𝑒 −6 − 𝑒 −8 ) = 3.343702
2
3.500168 − 3.343702
𝑒𝑟 = | | ∗ 100 = 4.47 %
3.500168

𝑏−𝑎
d) 𝐼 = (𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 𝑓(𝑥2 ))
6
2
𝐼 = (𝑓(0) + 4𝑓(2) + 𝑓(4))
3
2
𝐼 = (0 + 4(1 − 𝑒 −4 ) + (1 − 𝑒 −8 ))
3
2
𝐼 = (5 − 4𝑒 −4 − 𝑒 −8 ) = 3.289268
3
3.500168 − 3.289268
𝑒𝑟 = | | ∗ 100 = 6.03 %
3.500168

e) Regla de Simpson Múltiple n=4

pág. 11
𝑛−1 𝑛−2
𝑏−𝑎
𝐼= (𝑓(𝑥0 ) + 4 ∑ 𝑓(𝑥𝑖 ) + 2 ∑ 𝑓(𝑥𝑗 ) + 𝑓(𝑥𝑛 ))
3𝑛
𝑖=1,3,5 𝑗=2,4,6

1
𝐼 = (𝑓(0) + 4𝑓(1) + 4𝑓(3) + 2𝑓(2) + 𝑓(4))
3
1
𝐼 = (0 + 4(1 − 𝑒 −2 ) + 4(1 − 𝑒 −6 ) + 2(1 − 𝑒 −4 ) + (1 − 𝑒 −8 ))
3
1
𝐼 = (11 − 4𝑒 −2 − 4𝑒 −6 − 2𝑒 −4 − 𝑒 8 ) = 3.470592
3
3.500168 − 3.470592
𝑒𝑟 = | | ∗ 100 = 0.89 %
3.500168
Ejemplo
Calcula el volumen del solido de revolución generado al
girar alrededor del eje OX las siguientes
a) 𝑦 = 1 − 𝑥 2 , 𝑥 = −1, 𝑥 = 1.
1 1
b) 𝑦 = − 2 𝑥 2 + 2𝑥, 𝑦 = 2 𝑥.
Solución
a) 𝑦 = 1 − 𝑥 2 , 𝑥 = −1, 𝑧 = 1

Por el metodo de Simpson:


𝑛−1 𝑛
𝑏
ℎ 2ℎ 4ℎ
∫ 𝑓(𝑥)𝑑𝑥 ≈ (𝑓(𝑎) + 𝑓(𝑏)) + ∑ 𝑓(𝑥2𝑖 ) + ∑ 𝑓(𝑥2𝑖−1 )
𝑎 3 3 3
𝑖=1 𝑖=1

pág. 12
1
0.1 0.2
∫ (1 − 𝑥 2 )2 𝑑𝑥 ≈ [𝑓(−1) + 𝑓(1)] + [𝑓(−0.8) + 𝑓(−0.6) +
−1 3 3
0.4
⋯ 𝑓(0.6) + 𝑓(0.8)] + [𝑓(−0.9) + 𝑓(−0.7) + ⋯ + 𝑓(0.7) +
3
𝑓(0.9)] = 1.06669333
𝑉 = 1.06669333
1 1
b) 𝑦 = − 2 𝑥2 + 2𝑥, 𝑦 = 2𝑥

Hallando los puntos de intersección


1 1
= − 𝑥 2 + 2𝑥 = 𝑥 → −𝑥 2 + 4𝑥 = 𝑥 → 𝑥 2 − 3𝑥 = 0
2 2
𝑥(𝑥 − 3) = 0 → 𝑥 = 0 𝑦 𝑥 = 3
Luego
𝑏 − 𝑎 = 𝑛ℎ → 3 = 0.1 𝑛 → 𝑛 = 30
Por el metodo de Simpson:
𝑛−1 𝑛
𝑏
ℎ 2ℎ 4ℎ
𝑉 = ∫ 𝑓(𝑥)𝑑𝑥 ≈ (𝑓(𝑎) + 𝑓(𝑏)) + ∑ 𝑓(𝑥2𝑖 ) + ∑ 𝑓(𝑥2𝑖−1 )
𝑎 3 3 3
𝑖=1 𝑖=1
3
𝑥4 15 2 0.1
𝑉=∫ ( − 2𝑥 3 + 𝑥 ) 𝑑𝑥 ≈ [𝑓(3) + 𝑓(0)] +
0 4 4 3
0.2
+ [𝑓(0.2) + 𝑓(0.4) + 𝑓(0.6) + ⋯ + 𝑓(2.4) + 𝑓(2.6) + 𝑓(2.8)] +
3

pág. 13
0.4
+ [𝑓(0.1) + 𝑓(0.3) + ⋯ + 𝑓(2.5) + 𝑓(2.7) + 𝑓(2.9)]
3
𝑉 = 5.40001 𝜋

pág. 14

También podría gustarte