Tecnológico Nacional de México

Instituto Tecnológico De Pochutla

Trabajo: Ejercicios: ecuaciones de integración de Newton Cotes para funciones

continuas y discretas.

Materia: Métodos Numéricos.

Equipo:

Alcántara Narciso Marco Antonio……181160072

Bastida Ruiz Abimael…………………...181160084

Rojas Cárdenas José Alfredo…………181160064

Rosales Calderón José Enrique………181160069

Silva Hernández Jesús Manuel………181160071

Carrera: Ingeniería Civil

Semestre: 4°

Grupo: “A”

Docente: Ing. Ebert Matus Hernández.

Fecha de entrega: 20 de abril del 2020

1.- Utilizar la regla del trapecio para aproximar la integral:

𝑏
∫ 𝑒 𝑥 𝑑𝑥
𝑎

Donde con los siguientes datos resolver:

𝑎=0

𝑏=1
2
𝑓(𝑥) = 𝑒 𝑥

Por lo tanto, tenemos que:

1
2 𝑓(0) + 𝑓(1) 1+𝑒
∫ 𝑒 𝑥 𝑑𝑥 ≈ (1 − 0) [ ]= = 1.85914
0 2 2

2.- Usar la regla del trapecio para aproximar la integral:

𝑏
𝑓(𝑥)
∫ 𝑑𝑥
𝑎 𝑥

Con los datos:

𝑎=2

𝑏=4

𝑒𝑥
𝑓(𝑥) =
𝑥

Por lo tanto, tenemos que al sustituir los valores obtenemos lo siguiente:

4
𝑒𝑥 𝑓(2) + 𝑓(4) 𝑒2 𝑒4
∫ 𝑑𝑥 = (4 − 2) [ ]= + = 17.3441
2 𝑥 2 2 4
3.- consideremos la integral por el método de newton-Cotes (aproximación por
Trapecio)

Formula:

𝑏
𝐼 = ∫ 𝑓(𝑥) 𝑑𝑥
𝑎

ℎ
𝐹𝑇(𝑓) = [𝑓(𝑥0) + 𝑓(𝑥1)]
2

𝑥0 = 𝑎, 𝑥1 = 𝑎 + ℎ

ℎ =𝑏−𝑎

Datos:

𝑎=1

𝑏=2

1
𝑓(𝑥) =
𝑥

Sustitución:

2
1
𝐼=∫ 𝑑𝑥
1 𝑥

Solución:

2−1 1 1 3 3
𝐹𝑡(𝑓) = (1 + ) = ∗ = = 0.75
2 2 2 2 4

4.- consideremos la integral por newton-cotes (aproximación por Simpson).

Formulas:

𝑏
𝐼 = ∫ 𝑓(𝑥) 𝑑𝑥
𝑎

ℎ
𝐹𝑠(𝑓) = [𝑓(𝑥0) + 4𝑓(𝑥1) + 𝑓(𝑥2)]
3
 𝑥0 = 𝑎, 𝑥1 = 𝑎 + ℎ, 𝑥2 = 𝑎 + 2ℎ

𝑏−𝑎
ℎ=
2

Datos:

1
𝑎 = 1, 𝑏 = 2, 𝑓(𝑥) =
𝑥

Sustitución:

2
1
𝐼=∫ 𝑑𝑥
1 𝑥

Para la solución tenemos lo siguiente:

2−1
ℎ= = 0.5
2

𝑥0 = 1, 𝑥1 = 1.5, 𝑥2 = 2

0.5 1 1
𝐹𝑠(𝑓) = (1 + 4 ( ) + ) = 0.69444
3 1.5 2

2 𝑥3
5.- considere la siguiente integral ∫1 1 𝑑𝑥 por newton-cotes (método del
1+𝑥 2

trapecio)

Formulas:

𝑏
ℎ
𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 = [𝑓(𝑎) + 𝑓(𝑏)]
𝑎 2

𝑏−𝑎
ℎ=
𝑛

Dónde:

𝑎 = 1, 𝑏 = 2, 𝑛=1

Para resolver, sustituimos en las fórmulas:

 2−1
ℎ= =1
1

𝑓(𝑎) = 𝑓(1) = 0.5

𝑓(𝑏) = 𝑓(2) = 3.313708

Reemplazamos en la fórmula:

𝑏
ℎ
𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 = [𝑓(𝑎) + 𝑓(𝑏)]
𝑎 2

1
𝐼= [0.5 + 3.313708]
2

𝐼 = 1.906854

6.- integrar la siguiente función entre los limites a=-1 y b=1 utilizando 6 intervalos y
por newton-Cotes (método del trapecio Compuesto).

Formulas:

1 𝑥2
1 −
∫ 𝑒 2 𝑑𝑥
−1 √2𝜋

𝑏−𝑎
ℎ=
𝑛

como tenemos un numero de intervalos en este caso n=6 sustituimos para obtener
h.

1 − (−1) 1
ℎ= = = 0.333333
6 3

I Xi Xi xi
0 𝑎 = 𝑥0 −1 −1
1 𝑥1 = 𝑥0 + ℎ -2/3 -0.66666
2 𝑥2 = 𝑥1 + ℎ -1/3 -0.33333
3 𝑥3 = 𝑥2 + ℎ 0 0
4 𝑥4 = 𝑥3 + ℎ 1/3 0.33333
 5 𝑥5 = 𝑥4 + ℎ 2/3 0.6666
6 𝑥6 = 𝑥5 + ℎ 1 1

De esta manera la integral dividida en los intervalos quedaría:

𝑏=𝑥6
𝐼=∫ 𝑓(𝑥)𝑑𝑥
𝑎=𝑥0
𝑥1 𝑥2 𝑥3 𝑥4 𝑥5
= ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥
𝑥0 𝑥1 𝑥2 𝑥3 𝑥4
𝑥6
+ ∫ 𝑓(𝑥)𝑑𝑥
𝑥5

Que se puede sintetizar en la siguiente formula:

𝑏=𝑥6 𝑛
ℎ
𝐼=∫ 𝑓(𝑥)𝑑𝑥 = [𝑓(𝑥0 + 2 ∑ 𝑓(𝑥𝑖 + 𝑓(𝑥𝑛 )]
𝑎=𝑥0 2
𝑖=0

Para hallar los valores de 𝑓(𝑥𝑖 ), se evalúa la función en los puntos como se
muestra en la siguiente tabla:

𝑖 𝑥𝑖 𝑓(𝑥𝑖 )
0 -1 0.241971
1 -0.666667 0.319448
2 -0.333333 0.377383
3 0 0.398942
4 0.3333333 0.377383
5 0 0.319448
6 0.6666667 0.241971

Reemplazando en la fórmula:
𝑛

∑ 𝑓(𝑥𝑖 ) = 0.319448 + 0.377383 + 0.398942 + 0.377383 + 0.319448 = 1.792604

𝑖=0
 0.333333
𝐼= [0.241971 + 2(1.792604) + 0.241971] = 0.678191
2

𝐼 = 0.678191

7.- evalúe la integral de los siguientes datos tabulados, integrando por el método
del Simpson 1/3 compuesto.

Puntos 0 1 2 3 4 5
𝑥 0 0.1 0.2 0.3 0.4 0.5
𝑓(𝑥) 1 7 4 3 5 2

En primer lugar se verifica que el valor de “n”(5), no es múltiplo de 2, por lo que no

se podría aplicar el método. Pero se puede dividir la integral y aplicar distintos
métodos en los intervalos.

Luego la formula seria:

 𝑏=𝑥5
𝐼=∫ 𝑓(𝑥)𝑑𝑥
𝑎=𝑥0

ℎ ℎ
= [𝐹(𝑥0 + 4𝑓(𝑥1 ) + 𝑓(𝑥2 )] + [𝑓(𝑥3 + 4𝑓(𝑥3 ) + 𝑓(𝑥4 ]
3 3
ℎ
+ [𝑓(𝑥4 ) + 𝑓(𝑥5 )]
2

Reemplazando la tabla de datos:

Puntos 0 1 2 3 4 5
𝑥 0 01 0.2 03 0.4 0.5
𝑓(𝑥) 1 7 4 3 5 2

0.1 0.1 0.1

𝐼= [1 + 4(7) + 4] + [4 + 4(3) + 5] + [5 + 2]
3 3 2

𝐼 = 2.15

8.-evalue la siguiente integral, utilizando las fórmulas de cuadratura de gauss-

Legendre.

3
𝑒 𝑥 𝑠𝑒𝑛 𝑥
∫ 2
𝑑𝑥
0 1+𝑥

𝑏=3
𝑒 𝑥 𝑠𝑒𝑛 𝑥
𝐼=∫ 𝑑𝑥
𝑎=0 1 + 𝑥2
 2𝑥 − (𝑎 − 𝑏) 2𝑥 − (0 + 3) 2𝑥 − 3
𝑧= = =
𝑏−𝑎 3−0 3

Luego cambiando la variable y el diferencial:

3−0 3−0 3 3
𝑥= 𝑧+ = 𝑧+
2 2 2 2

𝑏−𝑎 3−0 3
𝑑𝑥 = 𝑑𝑧 = 𝑑𝑧 = 𝑑𝑧
2 2 2

Aplicando en la formula

𝑏=3
𝑒 𝑥 𝑠𝑒𝑛 𝑥 3 1 3 3
𝐼=∫ 𝑑𝑥 = ∫ 𝑓 ( 𝑧 + ) 𝑑𝑧
𝑎=0 1 + 𝑥2 2 −1 2 2

𝑛
3 3 3
𝐼 = ∑ 𝑤𝑖 𝑓 ( 𝑧𝑖 + )
2 2 2
𝑖=1

Puntos 𝑊𝐼 𝑍𝐼
2 𝑊1 = 𝑊2 = 1.0 −𝑍1 = 𝑍2 = 0.277350269

Evaluando los coeficientes z1:

𝑒 𝑥 𝑠𝑒𝑛 𝑥
𝑓(𝑥) =
1 + 𝑥2
3 3 3 3
3 3 𝑒 (2𝑧𝑖 +2) 𝑠𝑒𝑛 (2 𝑧𝑖 + 2)
𝑓 ( 𝑧𝑖 + ) =
2 2 3 3 2
1 + (2 𝑧𝑖 + 2)

3 3 3 3
𝑓 ( 𝑧𝑖 + ) = 𝑓 ( (−0.577350269) + ) = 0.796501031
2 2 2 2

3 3 3 3
𝑓 ( 𝑧2 + ) = 𝑓 ( (−0.577350269) + ) = 1.130596636
2 2 2 2

Remplazando en la fórmula:
 2
3 3 3 3 3 3 3 3
𝐼 = ∑ 𝑊𝐼 𝑓 ( 𝑧𝑖 + ) = [𝑊𝐼 𝑓 ( 𝑧𝑖 + ) + 𝑊2 𝑓 ( 𝑧2 + )]
2 2 2 2 2 2 2 2
𝐼=

3
𝐼= [(1.0)(0.796501031) + (1.0)(1.133059636)]
2

𝐼 = 2.890646501

9.- evalúe la siguiente integral doble, con 4 intervalos por el método del trapecio
compuesto:

2 4
∫ ∫ (𝑥 2 − 3𝑦 2 + 𝑥𝑦 3 )𝑑𝑥𝑑𝑦
−2 0

Primero se aplica a formula del trapecio compuesto en la integral interna para X,

con la variable Y como constante:

𝑦1 𝑥1 2 4
𝐼 = ∫ ∫ 𝑓(𝑥, 𝑦)𝑑𝑥𝑑𝑦 = ∫ [∫ ((𝑥 2 − 3𝑦 2 + 𝑥𝑦 3 )𝑑𝑥] 𝑑𝑦
𝑦0 𝑥0 −2 0

2 3
ℎ𝑥
𝐼 = ∫ [ [𝑓(𝑥0 , 𝑦) + 2 ∑ 𝑓(𝑥𝑖 , 𝑦) + 𝑓(𝑥4 , 𝑦)]]
−2 2 𝑖=1

Dónde:

𝑏 − 𝑎 𝑥4 − 𝑥0 (4 − 0)
ℎ𝑥 = = = =1
𝑛 4 4

Los puntos xi se muestran en la siguiente tabla:

𝑖 (𝑥𝑖 , 𝑦) (𝑥𝑖 , 𝑦)
0 (𝑥0 , 𝑦) (0,y)
1 (𝑥1 , 𝑦) = (𝑥0 + ℎ𝑥 , 𝑦) (1,y)
 2 (𝑥2 , 𝑦) = (𝑥1 + ℎ𝑥 , 𝑦) (2,y)
3 (𝑥3 , 𝑦) = (𝑥2 + ℎ𝑥 , 𝑦) (3,y)
4 (𝑥4 , 𝑦) = (𝑥3 + ℎ𝑥 , 𝑦) (4,y)

Desarrollando la sumatoria:

𝑛−1

∑ 𝑓(𝑥𝑖 , 𝑦) = 𝑓(𝑥1 , 𝑦) + 𝑓(𝑥2 , 𝑦) + 𝑓(𝑥3 , 𝑦)

𝑖=1

Luego se evalúan los puntos en la función f (xi,y):

𝑓(𝑥0 , 𝑦) = 𝑓(0, 𝑦) = 𝑥02 − 3𝑦 2 + 𝑥0 𝑦 3 = −3𝑦 2

𝑓(𝑥1 , 𝑦) = 𝑓(1, 𝑦) = 𝑥12 − 3𝑦 2 + 𝑥1 𝑦 3 = 𝑦 3 − 3𝑦 2 + 1

𝑓(𝑥2 , 𝑦) = 𝑓(2, 𝑦) = 𝑥22 − 3𝑦 2 + 𝑥2 𝑦 3 = 2𝑦 3 − 3𝑦 2 + 4

𝑓(𝑥3 , 𝑦) = 𝑓(3, 𝑦) = 𝑥32 − 3𝑦 2 + 𝑥3 𝑦 3 = 3𝑦 3 − 3𝑦 2 + 9

𝑓(𝑥4 , 𝑦) = 𝑓(4, 𝑦) = 4 − 3𝑦 2 + 𝑥4 𝑦 3 = 4𝑦 3 − 3𝑦 2 + 16

Reemplazando en la fórmula:

2
ℎ𝑥
𝐼=∫ [ [𝑓(𝑥0 , 𝑦) + 2[(𝑥1 , 𝑦) + 𝑓(𝑥2 , 𝑦) + 𝑓(𝑥3 , 𝑦)] + 𝑓(𝑥4 , 𝑦)]]𝑑𝑦
−2 2

2
1
𝐼 = ∫ [ [−3𝑦 2 + 2[𝑦 3 − 3𝑦 2 + 1 + 2𝑦 3 − 3𝑦 2 + 4 + 3𝑦 3 − 3𝑦 2 + 9] + 4𝑦 3 − 3𝑦 2
−2 2

+ 16]] 𝑑𝑦

2 2
𝐼=∫ [8𝑦 3 2
− 12𝑦 + 22]𝑑𝑦 = ∫ [𝑔(𝑦)]𝑑𝑦
−2 −2

Se tiene una función solo en función de y, que se puede definir como g(y), luego
aplicando la fórmula de un trapecio compuesta en la segunda integral, para la
variable y, se tiene:
 5
ℎ𝑦
𝐼= [𝑔(𝑦)] + 2 ∑ 𝑔(𝑦𝑖 ) + 𝑔(𝑦6 )]
2
𝑖=0

Dónde:

𝑦1 − 𝑦0 2 − (−2)
ℎ𝑦 = = =4
𝑛 1

Los puntos se muestran en la siguiente tabla:

𝑖 𝑦𝑖 𝑦𝑖 𝑔(𝑥𝑖 )
0 𝑦0 -2 -90
1 𝑦1 = 𝑦0 + ℎ𝑦 -1 2
2 𝑦2 = 𝑦1 + ℎ𝑦 0 22
3 𝑦3 = 𝑦2 + ℎ𝑦 1 18
4 𝑦4 = 𝑦4 + ℎ𝑦 2 38

Reemplazando en la forma:

ℎ𝑦
𝐼= [𝑔(𝑦0 ) + 2[𝑔(𝑦1 ) + 𝑔(𝑦2 ) + 𝑔(𝑦3 ) + 𝑔(𝑦4 )]
2

1
= [−90 + 2(2 + 22 + 18) + 38]
2

𝐼 = 16

10.- aproximar la siguiente integral, aplicando la regla de Simpson de 1/3y

subdividiendo en 5 intervalos.

1
2
∫ 𝑒 𝑥 𝑑𝑥
0
Tenemos que n= 5, y la partición que se genera es:

𝑃 = [0,0.2,0.4,0.6,0.8,1]

Los puntos medios de cada sub intervalos son

𝑃𝑀=[0.1,0.3,0.5,0.7,0.9]

Por lo tanto; sustituimos los datos en la fórmula para obtener:

1
2 𝑓(0) + 4[𝑓(0.1) + 𝑓(0.3) + 𝑓(0.5) + 𝑓(0.7) + 𝑓(0.9)] + 2[𝑓(0.2) + 𝑓(0.4) + ⋯ + 𝑓(0.8)] + 𝑓(1)
∫ 𝑒 𝑥 𝑑𝑥 ≈ (1 − 0) [ ]
0 6(5)

1 2 2 2 2 2
= [1 + 4[𝑒 (0.1) + 𝑒 (0.3) + 𝑒 (0.5) + 𝑒 (0.7) + 𝑒 (0.9) ]
30
2 2 2 2
+ 2[𝑒 (0.2) + 𝑒 (0.4) + 𝑒 (0.6) + 𝑒 (0.8) ] + 𝑒]

= 1.4626