1

Actividad Eje #4

Autores
Kelly Johana Landazuri Guanga
Eliana Andrea Perez Rosado
Salvador Velandia Montero

Docente
Richard Fabián Reyes Ramos

Facultad de Ingenierías
Fundación Universitaria del Área Andina
Análisis Numérico NRC 1476
2 de octubre de 2023
 2

Desarrollo de la actividad

1. Utilice la regla de la diferencia hacia adelante entre dos puntos para aproximar 𝑓′(1)
donde 𝑓(𝑥) = 𝐼𝑛𝑋, para:

a. ℎ = 0.1

b. ℎ = 0.01

c. ℎ = 0.001

Solución:

a. 𝒉 = 𝟎. 𝟏 es ⇒ incremento

se determina el punto de interés y la función:

𝑿𝟎 = 𝟏

𝒇(𝒙) = 𝑰𝒏𝑿

Se remplazan los valores en la fórmula:

𝒇(𝑿𝟎 + 𝒉) − 𝒇(𝑿𝟎 − 𝒉)
𝒇′(𝑿𝟎 ) =
𝟐𝒉

𝒇(𝟏 + 𝟎. 𝟏) − 𝒇(𝟏 − 𝟎. 𝟏) 𝒇(𝟏. 𝟏) − 𝒇(𝟎. 𝟗)

𝒇′(𝟏) = =
𝟐(𝟎. 𝟏) 𝟎. 𝟐

Se remplaza la función inicial:

𝒇(𝟏. 𝟏) − 𝒇(𝟎. 𝟗) 𝑰𝒏(𝟏. 𝟏) − 𝑰𝒏(𝟎. 𝟗)

𝒇′(𝟏) = =
𝟎. 𝟐 𝟎. 𝟐

𝟎. 𝟎𝟗𝟓𝟑𝟏𝟎𝟏𝟕𝟗𝟖𝟎𝟒𝟑𝟐𝟓 − (−𝟎. 𝟏𝟎𝟓𝟑𝟔𝟎𝟓𝟏𝟓𝟔𝟓𝟕𝟖𝟑) 𝟎. 𝟐𝟎𝟎𝟔𝟕𝟏

𝒇′(𝟏) = = = 𝟏, 𝟎𝟎𝟑𝟑𝟓𝟑
𝟎. 𝟐 𝟎. 𝟐
 3

Comprobamos:

Al momento de comprobar se deriva la función 𝒇(𝒙) = 𝑰𝒏𝑿, y luego remplazamos el punto de

observación.

𝟏
𝑰𝒏𝒙 =
𝒙

𝟏
= 𝟏 este valor se aproxima al obtenido por el método
𝒙

Error relativo:

𝟏−𝟏.𝟎𝟎𝟑𝟑𝟓𝟑
𝑬𝒓 = ∗ 𝟏𝟎𝟎 = 𝟎. 𝟑𝟑% este es error relativo que se halla
𝟏

b. 𝒉 = 𝟎. 𝟎𝟏 es ⇒ incremento

se determina el punto de interés y la función:

𝑿𝟎 = 𝟏
 4

𝒇(𝒙) = 𝑰𝒏𝑿

Se remplazan los valores en la fórmula:

𝒇(𝑿𝟎 + 𝒉) − 𝒇(𝑿𝟎 − 𝒉)
𝒇′(𝑿𝟎 ) =
𝟐𝒉

𝒇(𝟏 + 𝟎. 𝟎𝟏) − 𝒇(𝟏 − 𝟎. 𝟎𝟏) 𝒇(𝟏. 𝟏) − 𝒇(𝟎. 𝟗𝟗)

𝒇′(𝟏) = =
𝟐(𝟎. 𝟎𝟏) 𝟎. 𝟎𝟐

Se remplaza la función inicial:

𝒇(𝟏. 𝟎𝟏) − 𝒇(𝟎. 𝟗𝟗) 𝑰𝒏(𝟏. 𝟎𝟏) − 𝑰𝒏(𝟎. 𝟗𝟗)

𝒇′(𝟏) = =
𝟎. 𝟎𝟐 𝟎. 𝟎𝟐

𝟎. 𝟎𝟎𝟗𝟗𝟓𝟎𝟑𝟑𝟎𝟖𝟓𝟑𝟏𝟔𝟖𝟏 − (−𝟎. 𝟎𝟏𝟎𝟎𝟓𝟎𝟑𝟑𝟓𝟖𝟓𝟑𝟓𝟎𝟏)

𝒇′(𝟏) =
𝟎. 𝟎𝟐

𝟎. 𝟎𝟐𝟎𝟎𝟎𝟎𝟔𝟔𝟔𝟕𝟎𝟔𝟔𝟔𝟗𝟏
= = 𝟏, 𝟎𝟎𝟎𝟎𝟑𝟑𝟑𝟑
𝟎. 𝟎𝟐

Comprobamos:

Al momento de comprobar se deriva la función 𝒇(𝒙) = 𝑰𝒏𝑿, y luego remplazamos el punto de

observación.

𝟏
𝑰𝒏𝒙 =
𝒙

𝟏
= 𝟏 este valor se aproxima al obtenido por el método
𝒙

Error relativo:

𝟏−𝟏.𝟎𝟎𝟎𝟎𝟑𝟑𝟑𝟑
𝑬𝒓 = ∗ 𝟏𝟎𝟎 = 𝟎. 𝟎𝟎𝟑𝟑% este es error relativo que se halla
𝟏

c. 𝒉 = 𝟎. 𝟎𝟎𝟏 es ⇒ incremento
 5

se determina el punto de interés y la función:

𝑿𝟎 = 𝟏

𝒇(𝒙) = 𝑰𝒏𝑿

Se remplazan los valores en la fórmula:

𝒇(𝑿𝟎 + 𝒉) − 𝒇(𝑿𝟎 − 𝒉)
𝒇′(𝑿𝟎 ) =
𝟐𝒉

𝒇(𝟏 + 𝟎. 𝟏) − 𝒇(𝟏 − 𝟎. 𝟏) 𝒇(𝟏. 𝟏) − 𝒇(𝟎. 𝟗)

𝒇′(𝟏) = =
𝟐(𝟎. 𝟎𝟎𝟏) 𝟎. 𝟎𝟎𝟐

Se reemplaza la función inicial:

𝒇(𝟏. 𝟎𝟏) − 𝒇(𝟎. 𝟗𝟗) 𝑰𝒏(𝟏. 𝟎𝟏) − 𝑰𝒏(𝟎. 𝟗𝟗)

𝒇′(𝟏) = =
𝟎. 𝟎𝟎𝟐 𝟎. 𝟎𝟎𝟐

𝟎. 𝟎𝟎𝟗𝟗𝟓𝟎𝟑𝟑𝟎𝟖𝟓𝟑𝟏𝟔𝟖𝟏 − (−𝟎. 𝟎𝟏𝟎𝟎𝟓𝟎𝟑𝟑𝟓𝟖𝟓𝟑𝟓𝟎𝟏)

𝒇′(𝟏) =
𝟎. 𝟎𝟐

𝟎. 𝟎𝟐𝟎𝟎𝟎𝟎𝟔𝟔𝟔𝟕𝟎𝟔𝟔𝟔𝟗𝟏
= = 𝟏, 𝟎𝟎𝟎𝟎𝟑𝟑𝟑𝟑
𝟎. 𝟎𝟐

Comprobamos:

Al momento de comprobar se deriva la función 𝒇(𝒙) = 𝑰𝒏𝑿, y luego remplazamos el punto de

observación.

𝟏
𝑰𝒏𝒙 =
𝒙

𝟏
= 𝟏 este valor se aproxima al obtenido por el método
𝒙

Error relativo:
6

𝟏−𝟏.𝟎𝟎𝟎𝟎𝟑𝟑𝟑𝟑
𝑬𝒓 = ∗ 𝟏𝟎𝟎𝟎 = 𝟎. 𝟎𝟎𝟑𝟑% este es error relativo que se halla
𝟏
7

2. Utilice la fórmula de la diferencia centrada de tres puntos para la segunda derivada y

aproxime 𝑓′′(𝑥) = 𝐶𝑜𝑠𝑥, para:

a. ℎ = 0.1

b. ℎ = 0.01

c. ℎ = 0.001

Solución

a. h=0.1
se determina el punto de interés y la función:

𝑿𝟎 = 𝟏

𝒇(𝒙) = 𝑪𝑶𝑺𝑿

Se remplazan los valores en la fórmula:

′ (𝒙)
𝒇(𝒙 + 𝒉) − 𝒇(𝒙) 𝒉𝟐 ′′′′
𝒇 = 𝐥𝐢𝐦 − 𝒇 (𝒄)
𝒉→𝟎 𝒉𝟐 𝟏𝟐

En este caso no tendremos en cuenta el coeficiente de error:

8

𝒄𝒐𝒔(𝒙 + 𝒉) − 𝒄𝒐𝒔(𝒙) 𝒉𝟐
𝒇′ (𝟏) = 𝐥𝐢𝐦 −
𝒉→𝟎 𝒉𝟐 𝟏𝟐

𝒄𝒐𝒔(𝟏, 𝟏) − 𝒄𝒐𝒔(𝟏)
𝒇′ (𝟏) =
(𝟎. 𝟏)

𝒄𝒐𝒔(𝟏 + 𝟎, 𝟏) − 𝒄𝒐𝒔(𝟏)
𝒇′ (𝟏) =
(𝟎, 𝟏)𝟐

𝟎, 𝟗𝟗𝟗𝟖𝟏𝟓𝟕𝟏 − 𝟎, 𝟗𝟗𝟗𝟖𝟒𝟕𝟔𝟗
𝒇′ (𝟏) =
𝟎. 𝟎𝟏

−𝟎, 𝟎𝟎𝟎𝟎𝟑𝟏𝟗𝟖
𝒇′ (𝟏) =
𝟎, 𝟎𝟏

𝒇′ (𝟏) = −𝟎, 𝟎𝟎𝟑𝟏𝟗𝟖

b. h=0.01
se determina el punto de interés y la función:

𝑿𝟎 = 𝟏

𝒇(𝒙) = 𝑪𝑶𝑺𝑿

Se remplazan los valores en la fórmula:

𝒇(𝒙 + 𝒉) − 𝒇(𝒙) 𝒉𝟐 ′′′′

𝒇′ (𝒙) = 𝐥𝐢𝐦 − 𝒇 (𝒄)
𝒉→𝟎 𝒉𝟐 𝟏𝟐

En este caso no tendremos en cuenta el coeficiente de error:

𝒄𝒐𝒔(𝒙 + 𝒉) − 𝒄𝒐𝒔(𝒙)
𝒇′ (𝟏) = 𝐥𝐢𝐦
𝒉→𝟎 𝒉𝟐

𝒄𝒐𝒔(𝟏 + 𝟎, 𝟎𝟏) − 𝒄𝒐𝒔(𝟏)

𝒇′ (𝟏) =
(𝟎, 𝟎𝟏)𝟐
9

𝒄𝒐𝒔(𝟏, 𝟎𝟏) − 𝒄𝒐𝒔(𝟏)

𝒇′ (𝟏) =
𝟎, 𝟎𝟎𝟎𝟏

𝟎, 𝟗𝟗𝟗𝟖𝟒𝟒𝟔𝟑 − 𝟎, 𝟗𝟗𝟗𝟖𝟒𝟕𝟔𝟗
𝒇′ (𝟏) =
𝟎. 𝟎𝟎𝟎𝟏

−𝟎, 𝟎𝟎𝟎𝟎𝟎𝟑𝟏𝟗𝟖
𝒇′ (𝟏) =
𝟎, 𝟎𝟎𝟎𝟏

𝒇′ (𝟏) = −𝟎, 𝟎𝟑𝟎𝟔

c. h=0.001
se determina el punto de interés y la función:

𝑿𝟎 = 𝟏

𝒇(𝒙) = 𝑪𝑶𝑺𝑿

Se remplazan los valores en la fórmula:

′ (𝒙)
𝒇(𝒙 + 𝒉) − 𝒇(𝒙) 𝒉𝟐 ′′′′
𝒇 = 𝐥𝐢𝐦 − 𝒇 (𝒄)
𝒉→𝟎 𝒉𝟐 𝟏𝟐

En este caso no tendremos en cuenta el coeficiente de error:

𝒄𝒐𝒔(𝒙 + 𝒉) − 𝒄𝒐𝒔(𝒙)
𝒇′ (𝟏) = 𝐥𝐢𝐦
𝒉→𝟎 𝒉𝟐

𝒄𝒐𝒔(𝟏 + 𝟎, 𝟎𝟎𝟏) − 𝒄𝒐𝒔(𝟏)

𝒇′ (𝟏) =
(𝟎, 𝟎𝟎𝟏)𝟐

𝒄𝒐𝒔(𝟏, 𝟎𝟎𝟏) − 𝒄𝒐𝒔(𝟏)

𝒇′ (𝟏) =
𝟎, 𝟎𝟎𝟎𝟎𝟎𝟏
10

𝟎, 𝟗𝟗𝟗𝟖𝟒𝟕𝟑𝟗 − 𝟎, 𝟗𝟗𝟗𝟖𝟒𝟕𝟔𝟗
𝒇′ (𝟏) =
𝟎. 𝟎𝟎𝟎𝟎𝟎𝟏

−𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟑
𝒇′ (𝟏) =
𝟎, 𝟎𝟎𝟎𝟎𝟎𝟏

𝒇′ (𝟏) = −𝟎, 𝟑

3. Aplique la regla del trapecio para integrar numéricamente la función:

solución

𝒇(𝒙) = 𝟎, 𝟐 + 𝟐𝟓𝒙 − 𝟐𝟎𝟎𝒙𝟐 + 𝟔𝟕𝟓𝒙𝟑 − 𝟗𝟎𝟎𝒙𝟒 + 𝟒𝟎𝟎𝒙𝟓

Empezamos eligiendo los puntos 2 y 4

𝒃 𝟒

∫ 𝒇(𝒙) 𝒅𝒙 = ∫ 𝟎, 𝟐 + 𝟐𝟓𝒙 − 𝟐𝟎𝟎𝒙𝟐 + 𝟔𝟕𝟓𝒙𝟑 − 𝟗𝟎𝟎𝒙𝟒 + 𝟒𝟎𝟎𝒙𝟓 𝒅𝒙

𝒂 𝟐

Reemplazamos X por el límite inferior en la fórmula:

𝒇(𝒙) = 𝟎, 𝟐 + 𝟐𝟓𝒙 − 𝟐𝟎𝟎𝒙𝟐 + 𝟔𝟕𝟓𝒙𝟑 − 𝟗𝟎𝟎𝒙𝟒 + 𝟒𝟎𝟎𝒙𝟓

𝒇(𝒂) = 𝟎, 𝟐 + 𝟐𝟓(𝟐) − 𝟐𝟎𝟎(𝟐)𝟐 + 𝟔𝟕𝟓(𝟐)𝟑 − 𝟗𝟎𝟎(𝟐)𝟒 + 𝟒𝟎𝟎(𝟐)𝟓

𝒇(𝒂) = 𝟎, 𝟐 + 𝟐𝟎 − 𝟖𝟎𝟎 + 𝟓𝟒𝟎𝟎 − 𝟏𝟒𝟒𝟎𝟎 + 𝟏𝟐𝟖𝟎𝟎

𝒇(𝒂) = 𝟑𝟎𝟓. 𝟐

Reemplazamos X por el límite superior en la fórmula:

11

𝒇(𝒃) = 𝟎, 𝟐 + 𝟐𝟓(𝟒) − 𝟐𝟎𝟎(𝟒)𝟐 + 𝟔𝟕𝟓(𝟒)𝟑 − 𝟗𝟎𝟎(𝟒)𝟒 + 𝟒𝟎𝟎(𝟒)𝟓

𝒇(𝒃) = 𝟎, 𝟐 + 𝟏𝟎𝟎 − 𝟑𝟐𝟎𝟎 + 𝟒𝟑𝟐𝟎𝟎 − 𝟐𝟑𝟎, 𝟒𝟎𝟎 + 𝟒𝟎𝟗, 𝟔𝟎𝟎

𝒇(𝒃) = 𝟐𝟏𝟗. 𝟑𝟎𝟎, 𝟐

Aplicamos la siguiente fórmula:

𝒃
𝒇(𝒂) + 𝒇(𝒃)
∫ 𝒇(𝒙) 𝒅𝒙 = (𝒃 − 𝒂) ∗
𝟐
𝒂

𝟒
𝟑𝟎𝟓𝟎, 𝟐 + 𝟐𝟏𝟗𝟑𝟎𝟎, 𝟐
∫ 𝟎, 𝟐 + 𝟐𝟓𝒙 − 𝟐𝟎𝟎𝒙𝟐 + 𝟔𝟕𝟓𝒙𝟑 − 𝟗𝟎𝟎𝒙𝟒 + 𝟒𝟎𝟎𝒙𝟓 𝒅𝒙 = (𝟒 − 𝟐)
𝟐
𝟐

= 𝟐𝟐𝟐. 𝟑𝟓𝟎, 𝟒 𝒖𝒏𝒅𝒔𝟐

12

Comprobación:

Grafica
13

4. Evalué la siguiente integral con una sola aplicación de la regla compuesta del trapecio con
𝑛 = 2𝑦 4
solución
𝜋
2
∫ (6 + 3 𝑐𝑜𝑠𝑒𝑥) 𝑑𝑥
0

Primero identificamos los componentes del ejercicio n=2

𝒂=𝟎

𝝅
𝒃=
𝟐

𝒏=𝟐

Aplicamos la fórmula para hallar el intervalo de partición

𝝅
𝒃−𝒂 −𝟎 𝝅
∆𝒙 = = 𝟐 =
𝒏 𝟐 𝟒

Obtenemos los términos de partición

𝒙𝒊 (𝒙𝒊 − 𝟏) + ∆𝒙

𝒙𝟎 0

𝒙𝟏 𝜋
4

𝒙𝟐 𝜋
2
14

Luego reemplazos los valores de la ecuación de la siguiente forma:

𝒇(𝒙𝒊 ) 𝟔 + 𝟑 𝒄𝒐𝒔𝒙 𝑹𝒆𝒔𝒖𝒍𝒕𝒂𝒅𝒐

𝒇(𝒙𝟎 ) 𝟔 + 𝟑 𝒄𝒐𝒔 𝟎 9

𝒇(𝒙𝟏 ) 𝜋 8,121
𝟔 + 𝟑 𝒄𝒐𝒔
4

𝒇(𝒙𝟐 ) 𝜋 6
𝟔 + 𝟑 𝒄𝒐𝒔
2

Aplicamos los resultados en la siguiente formula;

𝝅
𝟐 ∆𝒙
∫ (𝟔 + 𝟑 𝒄𝒐𝒔𝒆𝒙) 𝒅𝒙 = (𝒇(𝒙𝟎 ) + 𝟐𝒇(𝒙𝟏 ) + 𝒇(𝒙𝟐 ))
𝟎 𝟐

𝝅 𝝅
𝟐
∫ (𝟔 + 𝟑 𝒄𝒐𝒔𝒆𝒙) 𝒅𝒙 = 𝟒 (𝟗 + 𝟐(𝟖, 𝟏𝟐𝟏) + 𝟔)
𝟎 𝟐

𝝅
(𝟑𝟏, 𝟐𝟒𝟐) = 𝟏𝟐, 𝟐𝟔𝟖 𝒖𝒏𝒅𝒔𝟐
𝟖
15

Grafica

Ahora identificamos los componentes del ejercicio 𝑛 = 4

𝒂=𝟎

𝝅
𝒃=
𝟐

𝒏=𝟒

Aplicamos la fórmula para hallar el intervalo de partición

𝝅
𝒃−𝒂 −𝟎 𝝅
∆𝒙 = = 𝟐 =
𝒏 𝟒 𝟖
16

Obtenemos los términos de partición

𝒙𝒊 (𝒙𝒊 − 𝟏) + ∆𝒙

𝒙𝟎 0

𝒙𝟏 𝜋
8

𝒙𝟐 𝜋
4

𝒙𝟑 3𝜋
8

𝒙𝟒 𝜋
2

Luego reemplazos los valores de la ecuación de la siguiente forma:

𝒇(𝒙𝒊 ) 𝟔 + 𝟑 𝒄𝒐𝒔𝒙 𝑹𝒆𝒔𝒖𝒍𝒕𝒂𝒅𝒐

𝒇(𝒙𝟎 ) 𝟔 + 𝟑 𝒄𝒐𝒔 𝟎 9

𝒇(𝒙𝟏 ) 𝜋 8,771
𝟔 + 𝟑 𝒄𝒐𝒔
8

𝒇(𝒙𝟐 ) 𝜋 8,121
𝟔 + 𝟑 𝒄𝒐𝒔
4

𝒇(𝒙𝟑 ) 3𝜋 7,148
𝟔 + 𝟑 𝒄𝒐𝒔
8

𝒇(𝒙𝟒 ) 𝜋 6
𝟔 + 𝟑 𝒄𝒐𝒔
2

Aplicamos los resultados en la siguiente formula:

17

𝝅
𝟐 ∆𝒙
∫ (𝟔 + 𝟑 𝒄𝒐𝒔𝒆𝒙) 𝒅𝒙 = (𝒇(𝒙𝟎 ) + 𝟐𝒇(𝒙𝟏 ) + 𝟐𝒇(𝒙𝟐 ) + 𝟐𝒇(𝒙𝟑 ) + 𝒇(𝒙𝟒 ))
𝟎 𝟐

𝝅 𝝅
𝟐
∫ (𝟔 + 𝟑 𝒄𝒐𝒔𝒆𝒙) 𝒅𝒙 = 𝟖 (𝟗 + 𝟐(𝟖, 𝟕𝟕𝟏) + 𝟐(𝟖, 𝟏𝟐𝟏) + 𝟐(𝟕, 𝟏𝟒𝟖) + 𝟔)
𝟎 𝟐

𝝅
(𝟔𝟑, 𝟎𝟖) = 𝟏𝟐, 𝟑𝟖𝟓 𝒖𝒏𝒅𝒔𝟐
𝟏𝟔

Grafica

5. Evalúe la siguiente integral con la regla de Simpson 1/3 𝑐𝑜𝑛 𝑛 = 4

18

3
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥
0

Solución

Calculamos el ancho del rectángulo

𝑎−𝑏 3−0
∆𝑥 = = = 0,75
𝑛 4

Determinamos las coordenadas de cada uno de los puntos

𝑥𝑖 = 𝑎 + 𝑖∆𝑥

𝑥0 0

𝑥1 0,75

𝑥2 1,5

𝑥3 2,25

𝑥4 3

Calculamos las alturas

𝑥2𝑒 𝑥

(0)2 𝑒 0 0,0000000000

(0,75)2 𝑒 0,75 1,1908125093

(1,5)2 𝑒 1,5 10,0838004083

(2,25)2 𝑒 2,25 48,0316626716

(3)2 𝑒 3 180,7698323087
19

Aplicamos la siguiente formula:

𝑥𝑛
𝑏−𝑎
∫ 𝑓(𝑥) 𝑑𝑥 = (𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + 𝑓(𝑥4 ))
𝑥0 3𝑛

Ahora realizamos el siguiente cuadro con el cálculo

𝑖 𝑥𝑖 Coeficiente 𝑓(𝑥𝑖 ) Coeficiente

*𝑓(𝑥𝑖 )

0 0 1 0,0000000000 0

1 0,75 4 1,1908125093 4,76325004

2 1,5 2 10,0838004083 20,1676008

3 2,25 4 48,0316626716 192,126651

4 3 1 180,7698323087 180,769832

Sumatoria total 397,827334

Luego hacemos la sumatoria y la multiplicamos para que nos del resultado:

3
∆𝑥
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥 = 𝑠𝑢𝑚𝑎𝑡𝑜𝑟𝑖𝑎 𝑡𝑜𝑡𝑎𝑙 ∗
0 3

3
0,75
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥 = 397,827334 ∗ = 99,45683346
0 3
 20

Grafica

6. Evalué siguiente integral mediante la regla de Simpson 1/3 𝑦 3/8 𝑐𝑜𝑛 𝑛 = 5

3
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥
0

Solución

Seguimos con la fórmula 1/3

𝑏−𝑎 3−0
∆𝑥 = = = 0,6
𝑛 5

(𝑥𝑛= 0 𝑓(𝑥0 = 𝑓(0) = 5 + 3 cos 0 = 5 + 3(1) = 8

(𝑥1= 0,6 𝑓(𝑥1 = 𝑓(0,6) = 5 + 3 cos 0,6 = 5 + 3(0,82533) = 7,476006

(𝑥2= 1,2𝑓(𝑥2 = 𝑓(1,2) = 5 + 3 cos 1,2 = 5 + 3(0,36235) = 6,087073

 21

(𝑥3= 1,8 𝑓(𝑥3 = 𝑓(1,8) = 5 + 3 cos 1,8 = 5 + 3(−0,22720) = 4,318393

(𝑥4= 2,4 𝑓(𝑥4 = 𝑓(2,4) = 5 + 3 cos 2,4 = 5 + 3(−0,73739) = 2,787818

(𝑥5= 3 𝑓(𝑥5 = 𝑓(3) = 5 + 3 cos 3 = 5 + 3(−0,98999) = 2,03002

3
𝑏−𝑎
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + 2𝑓(𝑥4 ) + 𝑓(𝑥5 )]
0 3∗𝑛

3
3−0
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [8 + 4(7,476006) + 2(6,087073) + 4(4,318393) + 2(2,787818) + (2,03002)]
0 3∗5

3
3−0
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [8 + 29,90402 + 12,17414 + 17,27357 + 5,57563 + 2,03002]
0 15

3
3−0
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [74,957408]
0 15

3
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ 14,99148
0

Seguimos con la fórmula 3/8

𝑏−𝑎 3−0
∆𝑥 = = = 0,6
𝑛 5

3
3∆𝑥
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [𝑓(𝑥0 ) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 3𝑓(𝑥3 ) + 3𝑓(𝑥4 ) + 𝑓(𝑥5 )]
0 8

3
3(0,6)
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [8 + 3(7,476006) + 3(6,087073) + 3(4,318393) + 3(2,787818) + (2,03002)]
0 8

3
1,8
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ ( ) [8 + 22,428018 + 18,261219 + 12,955179 + 8,363454 + 2,03002]
0 8

3
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ (0,225)[72,03789]
0
22

3
∫ (5 + 3 𝑐𝑜𝑠𝑥) 𝑑𝑥 ⇔ 16,2085
0

Grafica
23

Referencias

✓ Tarazona Oscar – Vega Bladimir Análisis numérico Digitk. Areandina.edu.co

https://digitk.areandina.edu.co/

✓ Rivera Diego – Diferenciación e integración numérica

https://www.youtube.com/watch?v=oM0sdG2QS3k

✓ Cedillo Hernández Antonio – Integración numérica – Regla de Simpson

https://www.youtube.com/watch?v=ypaNlJTPf9c