UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INDUSTRIAL Y SISTEMAS(FIIS)

ESCUELA PROFESIONAL:
INGENIERÍA DE SISTEMAS

CURSO: MÉTODOS NÚMERICOS

TEMA: “INTEGRACIÓN NUMÉRICA”

Prof. Dr.: GUILLERMO ANTONIO MAS AZAHUANCHE

2021

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
Integración de SIMPSON

18
19
20
21
22
23
24
EJERCICIOS METODO TRAPECIO
EJERICIO 1
Estime, mediante el método de los trapecios, el
valor de la siguiente integral definida, con una
partición de siete (7) elementos, y estime el error
cometido.

Se asume una partición regular, para ello determinamos la longitud de los

subintervalos,

y luego escribimos la partición

correspondiente,

Empleando la expresión de la fórmula de los trapecios para particiones

regulares, se obtiene una estimación de la integral definida

25
Para estimar el error cometido, empleamos la fórmula

correspondiente , con , , es decir que es

el máximo en valor absoluto de la segunda derivada de la función .
Determinamos entonces el valor de .
Para obtenemos:

; ; ;
Con la tercera derivada de la función, determinamos los números críticos
de la segunda derivada en el intervalo de integración , y calculamos
los extremos relativos que en ese intervalo se pueda alcanzar.
Números críticos de la segunda derivada:

Se evalúan estos números críticos, y los extremos del intervalo de

integración en la segunda derivada de la función, para obtener el máximo
en valor absoluto de :

Por tanto, .

Al sustituir en la fórmula de error, se obtiene .

26
EJERICIO 2
−3 𝑑𝑥
Calcule ∫0 usando la regla del trapecio con n=6.
𝑥 2 +16

Solucion:
𝑏−𝑎 3−0 1 𝑏−𝑎 1
𝐴𝑥 = = = luego =
𝑛 6 2 𝑛 4

Sustituyendo la formula
−3 𝑑𝑥 1
∫0 = 4 (𝑓(𝑥0)+ 2𝑓(𝑥1)+ 2𝑓(𝑥2)+ 2𝑓(𝑥4)+ 2𝑓(𝑥5)+ 𝑓(𝑥6))
𝑥 2 +16
1
= 4 (𝑓(0)+ 2𝑓(1/2)+ 2𝑓(1)+ 2𝑓(3/2)+ 2𝑓(2)+ 2(5/2)+2 𝑓(3))
1
=4 (0.0625+ 0.1230+ 0.1176 + 0.1096+ 0.1 + 0.09+0.04)
1
=4(0.6427) = 0.1607

EJERICIO 3
𝑥2
2 1 ( 2 )𝑑𝑥
Calcule ∫0 2𝜋 𝑒− usando la regla del trapecio con n=4.
√

Solucion:
𝑏−𝑎 2−0 1 𝑏−𝑎 1
𝐴𝑥 = = = 2 luego =4
𝑛 4 2𝑛

Sustituyendo la formula
𝑥2
2 1 ( 2 )𝑑𝑥 1
∫0 2𝜋 𝑒− = (𝑓(0)+ 2𝑓(1/2)+ 2𝑓(1)+ 2𝑓(3/2)+ 2𝑓(2))
√ 4
1 1 9
1
= 4 ((𝑒 0)+ 2(𝑒 −8 )+ 2(𝑒 −2 )+ 2(𝑒 −8 )+ 2𝑓(𝑒 −2)) = 1.1907
𝑥2
−2 1 ( 2 )𝑑𝑥
𝐿𝑢𝑒𝑔𝑜 ∫0 2𝜋 𝑒− = 0.47501.
√

EJERSICIO 4
Utilizando método trapecio
1
∫0 𝑥 4 𝑑𝑥 usando la regla del trapecio con n=9.
27
 i xi F(xi)
0 0 0
1 0.1111 1.52354833*10−4
2 0.2222 2.43767733*10−3
3 0.3333 0.0123407414
4 0.4444 0.0390028372
5 0.5555 0.0952217706
6 0.6666 0.1974518437
7 0.7777 0.3658039543
8 0.8888 0.6240453964
9 1 1
0.1111
I= 2 (0+2(1.52354833*10−4 + 2.43767733*10−3 +
0.0123407414 + 0.0390028372 + 0.0952217706 +
0.1974518437 + 0.3658039543+0.6240453964)
= 0.204030.
EJERSICIOS 5
3 1
∫0 𝑑𝑥 usando la regla del trapecio con n=9.
1+𝑥
𝑏−𝑎 3−0
h= = = 0.333
𝑛 9

i xi F(xi)
0 0 1
1 0.333 0.7501
2 0.666 0.6002
3 0.999 0.5002
4 1.332 0.4288
5 1.665 0.3752
6 1.998 0.33
7 2.331 0.3002
8 2.664 0.2729
9 2.99=3 0.250

28
 0.333
I= (1+0.25+2(0.6002 + 0.5002+ 0.4288 + 0.375 +
2
0.0952217706 + 0.1974518437 + 0.3658039543)
= 0.204030.
EJERSICIOS 6
2 𝑋 3 𝑑𝑥
∫1 1+𝑥2

Aplicamos la siguiente formula

𝑏 ℎ
I=∫𝑎 𝑓 (𝑥)𝑑𝑥 = 2 [𝑓 (𝑎) + 𝑓(𝑏 )]
𝑏−𝑎
h= 𝑛
2−1
h= =1
1

f(a)=f(1)=0.5
f(b)=f(2)=3.313708
reemplazando la formula
I=1/2(0.5+3.313708)
I=1.906854
METODO SIMPSON 1/3
Ejercicio 1

1 𝑥𝑑𝑥
Calcular el valor aproximado de la integral ∫0 (x + 1)(x + 2)
,
1
utilizando la regla de Simpson de 3 con n = 8.

Solución:
Ordenamos los valores de la función en un cuadro

29
De manera que
ℎ
I≈ [f(x0) + 4 (f(x1) + f(x3) + f(x5) + f(x7)) + 2 (f(x2) + f(x4) + f(x6)) +
3
f(x8)]
≈ 0.125 3 [4(0.05228 + 0.11483 + 0.14652 + 0.162319) + 2(0.0888 +
0.1333 + 0.15584) + 0.1666]
≈ 0.117773
Ejercicio 2.
1
Aproximar la siguiente integral, aplicando la regla de Simpson de 3 y
subdividiendo en 5 intervalos.
1
2
∫ 𝑒 𝑥 𝑑𝑥
0

Solución:
En este caso, tenemos que n=5, y la partición que se genera es: 𝑃 =
{0,0.2,0.4,0.6,0.8,1} Además, los puntos medios de cada
subintervalo son: 𝑃𝑀 = {0.1,0.3,0.5,0.7}
Por lo tanto, sustituimos los datos en la fórmula para obtener:

Ejercicio 3.
1
Aproximar la siguiente integral, utilizando la regla de Simpson de 3 y

subdividiendo en 4 intervalos.
4
𝑒𝑥
∫ 𝑑𝑥
2 𝑥
30
Solución:
En este caso, tenemos que n=4, y la partición que se genera es: 𝑃 =
{2,2.5,3,3.5,4} Además, los puntos medios de cada
subintervalo son: 𝑃𝑀 = {2.25,2.75,3.25,3.75}
Sustituyendo todos estos datos en la fórmula obtenemos la siguiente
aproximación:

Ejercicio 4.
1
Estime mediante el regla de Simpson de 3
el valor de la
siguiente integral definida, con una partición de siete (7) elementos, y
estime el error cometido.
Solución:
Para aplicar el método de Simpson, la partición debe ser regular y el
número de subintervalos debe ser par, como la partición debe tener
siete elementos se cumple está condición y
determinamos, entonces, la longitud de los
subintervalos:
Construyendo la partición correspondiente,

Empleando la fórmula del método de Simpson, se obtiene una

estimación de la integral definida

31
Ejercicio 5.
Teniendo en cuenta que no es conocida una primitiva de la función
2 1 2
𝑓 (𝑥) = 𝑒 𝑥 , calcúlese el valor de la integral definida ∫0 𝑒 𝑥 𝑑𝑥 Con un
error menor a 0.003.
Solución:
2
La función con la que debemos trabajar es 𝑒 𝑥 . Aplicando la regla de
Simpson de cometemos un error dado por
ℎ4
𝐸 ≤ |(𝑏 − 𝑎) 180 𝑀4 | , 𝑀4 ≥ f (4) (x), ∀x ∈ [0, 1]

Calcularemos las derivadas correspondientes:

2
𝑓′(𝑥) = 2𝑥𝑒 𝑥 (𝑏 − 𝑎)5
𝐸 (𝑁) ≤ ( 228)
2 180𝑁4
𝑓′′(𝑥) = 2(1 + 2𝑥 2 )𝑒 𝑥 E(1) ≈ 1.2666
𝑓′′′(𝑥) = 4(3𝑥 + 2𝑥 3 )𝑒 𝑥
2
19 E(2) ≈ 0.0791
= E(3) ≈ 0.0156
𝑓 ′ ′′′(𝑥) = 4(4𝑥 4 + 12𝑥 2 + 3)𝑒 𝑥
2 15𝑁4
E(4) ≈ 0.0049
{ E(5) ≈ 0.0020

Se puede observar que 𝑓 ′ ′′′(𝑥)

es creciente en [0, 1] de modo
que el máximo valor de dicha
funci´on coincide con el valor en
′′′(1)
x = 1, esto es 𝑓 ′ = 4𝑒(4 +
12 + 3) < 4𝑥3𝑥19 = 228
, por lo que consideraremos
que 𝑀4 ≤ 228 Por ello,

32
 De
modo
que
para
que el
número de subintervalos sea par hemos de tomar N = 6
Finalmente:
I ≈ h 3 [f(x0) + 4 (f(x1) + f(x3) + f(x5)) + 2 (f(x2) + f(x4)) + f(x8)]
≈ 1/6 3 [1 + 4(1.02817 + 1.28403 + 2.0026) + 2(1.11752 +
1.55962) + 2.71828]
≈ 1.4628
Ejercicio 6.
Una cuerda vibra adoptando la forma, 𝑦 = 𝑠𝑒𝑛(𝑥) entre las
abscisas 𝑥 = 0 𝑦 𝑥 = 4 en un instante 𝑡0. Calcúlese
aproximadamente la longitud de la cuerda, utilizando la regla de
Simpson con n = 8.
Solución:
Dado que tenemos que calcular la longitud de la función f(x) = sen
x, entre x = 0 y x = 4, aplicaremos la fórmula

𝑏 2 4
𝐿 = ∫𝑎 √1 + (𝑓 ′ (𝑥)) 𝑑𝑥 = ∫0 √1 + 𝑐𝑜𝑠 2 𝑥𝑑𝑥 , que nos proporciona la
integral que debemos estimar numéricamente mediante la regla de
Simpson con n = 8 como propone el enunciado.

𝑔(𝑥) = √1 + 𝑐𝑜𝑠 2 𝑥

I ≈ h 3 [g(x0) + 4 (g(x1) + g(x3) + g(x5) + g(x7)) + 2 (g(x2) + g(x4) +

g(x6)) + g(x8)]
0.5
≈ [1.41421 + 4(1.33047 + 1.0025 + 1.28134 + 1.37002) +
3
2(1.13663 + 1.08313 + 1.40715) + 1.19468]
≈ 4.96667

33
EJERCICIOS METODO SIMPSON 3/8
EJERCICIO 1
Hallar una aproximación a la integral por medio del método de
Simpson 3/8
𝑥 2 𝑑𝑥
∫ 1 [1,2]
1 + 𝑥2
Se procede a hallar el punto medio a partir de los intervalos dados
a=1 b=2 n=3
2−1
ℎ= = 0.33, 𝑥0 = 1
3

𝑥1 = 1.33 𝑒𝑠 𝑒𝑙 𝑝𝑢𝑛𝑡𝑜 𝑚𝑒𝑑𝑖𝑜

𝑥2 = 1.66
Reemplazamos la formula inicial
𝑥0 = 1
12 𝑑𝑥
1 = 0.5
∫
1+ 12

Luego para hallar el nuevo valor de x utilizaremos la formula

𝑓 = (𝑥0 + ℎ)
𝑥1 = 1 + 0.33 = 1.33
1.332 𝑑𝑥
1 = 0.48
∫
1+ 1.332

𝑥2 = 1.33 + 0.33 = 1.66

1.662 𝑑𝑥
∫ 1 = 0.44
1+ 1.662
Tenemos que por el método de Simpson 3/8 el resultado de la
integral es 0.44

34
EJERCICIO 2
2 𝑥 3 𝑑𝑥
Obtenga ∫1 1
1+𝑥 2
Se utilizan las siguientes formulas:

Donde
a=1, b=2, n=3
Se debe notar que para el método Simpson 3/8 simple el valor de
n es igual a 3
Luego

Lo que se puede resumir en la siguiente tabla

Reemplazando en la formula
3(0.333333)
𝐼= [0.5 + 3(1.10092 + 2.020793) + 3.313708]
8
𝐼 = 1.647045

35
EJERCICIO 3
Integre la siguiente función entre los limites a=-1 y b=1, utilizando
𝑥2
1 1
6 intervalos ∫−1 2𝜋 𝑒 − 2 𝑑𝑥
√

36
 𝐼 = 0.682851
EJERCICIO 4
Evalue la integral de los siguientes datos tabulados:

𝐼 = 2.220833

37
EJERCICIO 5
Aproximar la siguiente integral, usando la regla de Simpson de 3/8
4
∫ 𝑒 𝑥 𝑙𝑛𝑥𝑑𝑥
1
Solucion:
En este caso tenemos los siguientes datos:

EJERCICIO 6
4
Aproximar la siguiente integral: ∫1 𝑒 𝑥 𝑙𝑛𝑥𝑑𝑥
Aplicando la regla de simpsion de 3/8, y subdiviendo en 3
intervalos.
Solucion:
Identificamos n=3 y la partición correspondiente:
P={1,2,3,4}
Al considerar los puntos que dividen en tres partes iguales a cada
subintervalo, tenemos los siguientes datos :

38
EJERCICIOS ADICIONALES

39
40
41
42
43