Integración Numérica 2022-2

Universidad Nacional Agraria
La Molina
FA C U LTA D E D E C I E N C I A S
DEPARTAM E N TO ACADÉM ICO DE M ATEM ÁTICA
CURSO : CÁLCULO PARA CIENCIA S

CICLO : 2 0 2 2 - II PRESENC I A L - VIRTUA L
TEMA: INTEGRACIÓN NUMÉRICA

Capítulo II: INTEGRACIÓN NUMÉRICA
2.1 Definiciones básicas para Integración Numérica

Partición de un Intervalo
Sea 𝐼 = 𝑎, 𝑏 un intervalo cerrado, llamaremos una partición de 𝐼 a un conjunto de puntos

de dicho intervalo, como: 𝑃 = 𝑥0 , 𝑥1 , 𝑥2 , … , 𝑥𝑛 −1 , 𝑥𝑛 , tal que 𝑥0 = 𝑎; 𝑥𝑛 = 𝑏, donde
𝑎 = 𝑥0 < 𝑥1 < 𝑥2 < ⋯ < 𝑥𝑛 −1 < 𝑥𝑛 = 𝑏.
Se obtiene n sub intervalos de la partición P, 𝑎, 𝑥1 , 𝑥1 , 𝑥2 , 𝑥2 , 𝑥3 , … , 𝑥𝑛 −1 , 𝑏 ; en

forma simplificada, un sub intervalo cualquiera se expresa como:
𝑥𝑖−1 , 𝑥𝑖 , 𝑖 = 1, 2, 3, … , 𝑛
La longitud de cada sub intervalo es, 𝑥𝑖 − 𝑥𝑖−1 = ∆𝑥𝑖 y la sumatoria es

𝑛
෍ ∆𝑥𝑖 = 𝑏 − 𝑎,
𝑖=1
INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL
Teorema
Sea f una función continua en un intervalo cerrado 𝑎, 𝑏 , M y m son valores
máximos y mínimos absolutos respectivamente de f en 𝑎, 𝑏 , entonces se
cumple:
𝑏
𝑚 𝑏 − 𝑎 ≤ න 𝑓 𝑥 𝑑𝑥 ≤ 𝑀 𝑏 − 𝑎 .
𝑎
Ejemplo
Aplicando el teorema anterior, determine en qué intervalo se encuentra

el valor de la siguiente integral:
4
න 1 + 𝑥 4 𝑑𝑥
0

Solución
Notemos que 𝑓 𝑥 = 1 + 𝑥 4 , 0 ≤ 𝑥 ≤ 4, como la función 𝑓 es continua en
0,4 , los valores extremos absolutos de 𝑓 son:
𝑚 = 𝑓 0 = 1 y 𝑀 = 𝑓 4 = 257 = 16.0312.
Aplicando el teorema anterior obtenemos:
4
1 4−0 ≤න 1 + 𝑥 4 𝑑𝑥 ≤ 257 4 − 0
0
4
4.00 ≤ න 1 + 𝑥 4 𝑑𝑥 ≤ 64.1249
0

Teorema del Valor Medio para Integrales Definidas
Sea f una función continua y en un intervalo cerrado 𝑎, 𝑏 . Entonces, existe 𝑐 ∈ 𝑎, 𝑏

tal que:
𝑏
1
𝑓 𝑐 = න 𝑓 𝑥 𝑑𝑥
𝑏−𝑎 𝑎
o
𝑏
න 𝑓 𝑥 𝑑𝑥 = 𝑏 − 𝑎 𝑓 𝑐 .
𝑎
Al valor f(c) se le denomina “valor medio” de la función f en 𝑎, 𝑏 .
Para la interpretación geométrica se requiere además que 𝑓 𝑥 > 0.
Ejemplo 1
10
Establezca un intervalo en el cual se encuentra el valor de la integral ‫׬‬3 𝑓 𝑥 𝑑𝑥, si
𝑓 𝑥 = 𝑙𝑛 𝑥 luego determine el valor de la altura media de 𝑓.

Solución
1
La función 𝑦 = 𝑓 𝑥 = 𝑙𝑛 𝑥 es continua para, 𝑥 > 0; 𝑓 ′ 𝑥 = > 0, 𝑥 > 0;
𝑥
máximo absoluto 𝑀 = 𝑓 10 = 𝑙𝑛 10 y mínimo absoluto 𝑚 = 𝑓 3 = 𝑙𝑛 3 .
10
10 − 3 𝑙𝑛3 ≤ න 𝑙𝑛 𝑥 𝑑𝑥 ≤ 10 − 3 𝑙𝑛10
3
Operando se tiene;
10
7.69 ≤ න 𝑙𝑛 𝑥 𝑑𝑥 ≤ 16.12
3
Entonces el intervalo donde se encuentra el valor de la integral es 7.69,16.12
Por el teorema del valor medio de la integral definida para 𝑐 ∈ 3,10 se tiene:
10
1 1
𝑦ത = 𝑓 𝑐 = න 𝑙𝑛 𝑥 𝑑𝑥 = 10𝑙𝑛10 − 3𝑙𝑛3 − 7 = 1.82
10 − 3 3 7

geométricamente, el valor de 𝑓 𝑐 = 𝑦ത = 1.82, es la altura media de la función
𝑦 = 𝑓 𝑥 en el intervalo 3,10 , como se muestra en la siguiente figura.

Ejemplo 2. Calcule el valor medio de la función 𝑓 𝑥 = 2 + 3𝑥 3 , en el
intervalo −1,2 . Luego determine el(los) número(s) que satisfaga el teorema
de valor medio.
Solución.
Usando el Teorema del valor medio de la integral definida, se tiene:
2
1 1 69 23
𝑦ത = 𝑓 𝑐 = න 2 + 3𝑥 3 𝑑𝑥 = = = 5.75
2 − (−1) −1 3 4 4
69
4
Como 𝑓 𝑐 = 2 + 3𝑐 3, calculamos
1 2 23
න 2 + 3𝑥 3 𝑑𝑥 = 𝑓 𝑐 = 2 + 3𝑐 3 =
3 −1 4

3
Obteniendo el valor que satisface el teorema de valor medio 𝑐 = 5/4.
Ejemplo 3. El crecimiento de una determinada planta en los primeros 30 días

de trasplantado sigue el modelo 𝑓 𝑡 = 1.2 ln 𝑡 + 1 + 20, en donde 𝑡 se
mide en días y 𝑓 𝑡 es el tamaño de la planta, en centímetros a los 𝑡 días que
se trasplantó.
a) ¿Cuál es el tamaño promedio de la planta en los 30 días?
b) ¿A los cuantos días del trasplante alcanza el tamaño promedio?
Solución.
a) Por el TVM de la integral definida se tiene
30
1
𝑦ത = 𝑓 𝑐 = න 1.2 ln 𝑡 + 1 + 20 𝑑𝑡
30 − 0 0

1 6 6 94 30 691.74
𝑦ത = 𝑡𝑙𝑛 𝑡 + 1 + ln 𝑡 + 1 + 𝑡൨ = = 23.06 𝑐𝑚.
30 5 5 5 0 30
El tamaño promedio de la planta en los 30 días es de 𝑦ത = 23.06 𝑐𝑚, como se

muestra en la siguiente figura.

a) Para encontrar el día se calcula el valor de 𝑡 = 𝑐 donde ocurre el tamaño medio,
será en (𝑐, 𝑓 𝑐 ) como se muestra en la figura anterior:
30
1
𝑓 𝑐 = 1.2 ln 𝑡 + 1 + 20 = න 1.2 ln 𝑡 + 1 + 20 𝑑𝑡 = 23.06
30 − 0 0
resolviendo se tiene 𝑐 = 11.81 días, entonces el día aproximado donde la planta

tiene una altura media de 23.06 cm es a los 11 días con 19.44 horas
aproximadamente.

Integración Numérica
𝑏
Sea f una función continua en 𝑎, 𝑏 . Para calcular ‫ 𝑥𝑑 𝑥 𝑓 𝑎׬‬se usa el Teorema
Fundamental del Cálculo, pero cuando no es posible encontrar la antiderivada de f,
debemos aproximar el valor de esta integral por medio de integración numérica,
usando sumas de Riemann de 𝑓 para una partición regular. La elección de n depende
de la precisión que nosotros queramos para nuestra integral. Es decir la integral se
aproxima como:
𝑏 𝑛
න 𝑓 (𝑥 )𝑑𝑥 ≈ ෍ 𝑓 (𝑤𝑖 )Δ𝑥

𝑎 𝑖=1
donde 𝑤𝑖 es un número arbitrario en el i-ésimo subintervalo 𝑥𝑖−1 , 𝑥𝑖 de la partición

𝑃 = 𝑥0 , 𝑥1 , 𝑥2 , … , 𝑥𝑖−1 , 𝑥𝑖 , … , 𝑥𝑛 .

En la siguiente figura.
Si 𝑤𝑖 = 𝑥𝑖−1, calculamos el valor de f en
el extremo izquierdo del sub intervalo
𝑥𝑖−1 , 𝑥𝑖 , entonces la integral definida se
aproxima como:
𝑏 𝑛
න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖−1 Δ𝑥
𝑎 𝑖=1
Si 𝑤𝑖 = 𝑥𝑖 calculamos f en el extremo derecho del sub intervalo 𝑥𝑖−1 , 𝑥𝑖 , entonces la

integral definida se aproxima como:
𝑏 𝑛
න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖 Δ𝑥
𝑎 𝑖=1

Una mejor aproximación se obtiene cuando se usa el promedio de las sumas
de la derecha e izquierda según:
𝑏 𝑛 𝑛 𝑛 𝑛
1 Δ𝑥
න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖−1 Δ𝑥 + ෍ 𝑓 𝑥𝑖 Δ𝑥 = ෍ 𝑓 𝑥𝑖−1 + ෍ 𝑓 𝑥𝑖
𝑎 2 2
𝑖=1 𝑖=1 𝑖=1 𝑖=1
Ejemplo
Calcule el valor aproximado de las siguientes integrales:
3 3
𝑖) න 1 + 𝑥2 𝑑𝑥 , para 𝑛 = 6 𝑖𝑖) න 10 − 𝑥 2 𝑑𝑥 , para 𝑛 = 6
0 1

Solución
3−0
para 𝑛 = 6; 𝑎 = 0 y 𝑏 = 3; se tiene Δ𝑥 = ; 𝑓 𝑥 = 1 + 𝑥2
6
a) Aproximamos la integral evaluando f en el

extremo izquierdo:
3 6
න 1 + 𝑥 2 𝑑𝑥 ≈ ෍ f xi−1 Δx
0 i=1
෍ 𝑓 𝑥𝑖−1 Δ𝑥 = 𝑓 𝑥0 + 𝑓 𝑥1 + 𝑓 𝑥2 + 𝑓 𝑥3 + 𝑓 𝑥4 + 𝑓 𝑥5 Δ𝑥
𝑖=1

En la siguiente tabla se muestra los valores de 𝑓 para cada 𝑥𝑖−1 y la suma total.
𝑖−1 𝑥𝑖−1 𝑓 𝑥𝑖−1
0 0.0 1.0000
1 0.5 1.2500
2 1.0 2.0000
3 1.5 3.2500
4 2.0 5.0000
5 2.5 7.2500
6
෍ 𝑓 𝑥𝑖−1 = 19.7500
𝑖=1
Por lo tanto la integral aproximada por la izquierda es:

3
න 1 + 𝑥 2 𝑑𝑥 ≈ 19.75 0.5 = 9.875
0

b) Aproximamos la integral evaluando f en el extremo derecho:
3 6
න 1 + 𝑥 2 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖 Δ𝑥
0 𝑖=1
En la siguiente tabla se muestra los valores de f para cada 𝑥𝑖 y la suma total.

:
i xi f(xi)
1 0.5 1.2500
2 1.0 2.0000
3 1.5 3.2500
4 2.0 5.0000
5 2.5 7.2500
6 3.0 10.0000
6
෍ 𝑓 𝑥𝑖 = 28.7500
𝑖=1

Por lo tanto la integral aproximada por la derecha es:
3
න 1 + 𝑥 2 𝑑𝑥 ≈ 28.75 0.5 = 14.375
0
Una buena aproximación es:
3
14.375 + 9.875
න 1+ 𝑥2 𝑑𝑥 ≈ = 𝟏𝟐. 𝟏𝟐𝟓
0 2

c) Calculando 𝑓 en el punto medio de cada subintervalo:
3 𝑛
𝑥𝑖−1 + 𝑥𝑖
න 1+ 𝑥2 𝑑𝑥 ≈ ෍ 𝑓 𝑤𝑖 Δ𝑥, 𝑤𝑖 =
0 2
𝑖=1
𝑖 𝑥𝑖 𝑥𝑖−1 + 𝑥𝑖 𝑓 (𝑤𝑖 )
𝑤𝑖 =
2
0 0.0
1 0.5 0.25 1.0625
2 1.0 0.75 1.5625
3 1.5 1.25 2.5625
4 2.0 1.75 4.0625
5 2.5 2.25 6.0625
6 3.0 2.75 8.5625
6
෍ 𝑓 (𝑤𝑖 ) = 23.87500
𝑖=1
Con los valores de la tabla obtenemos:

3
න 1 + 𝑥 2 𝑑𝑥 ≈ 23.875 0.5 = 11.94
0

3
𝑖𝑖) න 10 − 𝑥 2 𝑑𝑥 , para 𝑛 = 6
1
Como: 𝑛 = 6; 𝑎 = 1; 𝑏 = 3; 𝑓 𝑥 = 10 − 𝑥 2 ,
1
entonces Δ𝑥 =
3
a) Aproximamos la integral evaluando f en el extremo izquierdo:
3 𝑛
න 10 − 𝑥 2 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖−1 Δ𝑥
1 𝑖=1
෍ 𝑓 𝑥𝑖−1 Δ𝑥 = 𝑓 𝑥0 + 𝑓 𝑥1 + 𝑓 𝑥2 + 𝑓 𝑥3 + 𝑓 𝑥4 + 𝑓 𝑥5 Δ𝑥
𝑖=1

En la siguiente tabla se muestra los valores de 𝑓 para cada 𝑥𝑖−1 y la suma total.
𝑖−1 𝑥𝑖−1 𝑓 𝑥𝑖−1
0 1.00 3.0000
1 1.33 2.8674
2 1.67 2.6874
3 2.00 2.4495
4 2.33 2.1344
5 2.67 1.6997
6
14.8384
෍ 𝑓 𝑥𝑖−1 =
𝑖=1
Por lo tanto la integral aproximada por la izquierda es:

3
න 10 − 𝑥 2 𝑑𝑥 ≈ 4.95
1

b) Aproximamos la integral evaluando f en el extremo derecho:
3 𝑛
න 10 − 𝑥 2 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖 Δ𝑥
1 𝑖=1
En la siguiente tabla se muestra los valores de f para cada 𝑥𝑖 y la suma total.
:
i xi f(xi)
t 1.33 2.8674
2 1.67 2.6874
3 2.00 2.4495
4 2.33 2.1344
5 2.67 1.6997
6 3.00 1.0000
6
෍ 𝑓 𝑥𝑖 = 12.8384
𝑖=1

Por lo tanto la integral aproximada por la derecha es:
3
1
න 10 − 𝑥 2 𝑑𝑥 ≈ 12.8384 ( ) = 4.28
1 3
Una mejor aproximación es:
3
4.95 + 4.28
න 10 − 𝑥 2 𝑑𝑥 ≈ = 4.61
1 2

2.2 Método trapezoidal y método trapezoidal compuesto.
Estimación del error
Método trapezoidal o Regla del trapecio.
Es una técnica de integración numérica para calcular cierta integral definida de una
función 𝑓 en un intervalo cerrado 𝑎, 𝑏 , mediante una aproximación al área de un
trapecio formado por los puntos:
𝐴 = 𝑎, 0 ; 𝑃 = 𝑎, 𝑓 𝑎 ; 𝑄 = 𝑏, 𝑓 𝑏 y 𝐵 = (𝑏, 0) ; como se muestra en la
siguiente figura:
y
𝑦=𝑓 𝑥
𝐴 𝐵
0 𝑎 x
𝑏

2.2 Método trapezoidal y método trapezoidal compuesto. Estimación del error
Teorema
Sea 𝑓 una función positiva continua en un intervalo cerrado 𝑎, 𝑏 , 𝑃𝑄 un
segmento cuyos extremos pertenecen a la gráfica de y = 𝑓(𝑥) mostrada en la
figura anterior. La integral definida de 𝑓 en 𝑎, 𝑏 es aproximadamente igual
al área del trapecio APQB, es decir:
𝑏
𝑓 𝑎 +𝑓 𝑏
න 𝑓 𝑥 𝑑𝑥 ≈ 𝑏 − 𝑎 . … (1)
𝑎 2
donde las bases son 𝐴𝑃 = 𝑓 𝑎 y 𝐵𝑄 = 𝑓(𝑏) y de altura (𝑏 − 𝑎).
Nota 𝑏
𝑓 𝑎 + 𝑓(𝑏)
න 𝑓 𝑥 𝑑𝑥 = 𝑏 − 𝑎 +𝐸
𝑎 2
Si 𝑓 tiene segunda derivada continua en un intervalo cerrado a, b , el error (E)

estimado por la aproximación está dado por:
ℎ3
𝐸 ≤ 𝑚á𝑥 𝑓 ′′ 𝑥 , 𝑥 ∈ 𝑎, 𝑏
12
Ejemplo
Use la Regla del Trapecio para:
4
1
a) Aproximar el valor de la siguiente integral, න 𝑑𝑥.
2 𝑙𝑜𝑔2 𝑥
b) Estimar el error máximo que se comete por la aproximación.
Solución
1 𝑙𝑛2
Tenemos: 𝑎 = 2, 𝑏 = 4, 𝑓 𝑥 = = ,𝑥 ∈ 2,4 además 𝑓 2 = 1,
𝑙𝑜𝑔2 𝑥 𝑙𝑛𝑥
1
𝑓 4 = , reemplazando estos valores en la aproximación (1).
2

4
1 𝑓 2 +𝑓 4
න 𝑑𝑥 ≈ (4 − 2)
2 𝑙𝑜𝑔2 𝑥 2
1
4
1 1+ 4
1
න 𝑑𝑥 ≈ (2) 2 → න 𝑑𝑥 ≈ 1.50
2 𝑙𝑜𝑔2 𝑥 2 2 𝑙𝑜𝑔2 𝑥
Cálculo del error máximo (E).
𝑙𝑛2 ′
𝑙𝑛2 2 + 𝑙𝑛𝑥
𝑓 𝑥 = , 𝑥 ∈ 2,4 ; 𝑓 𝑥 = − 2; 𝑓 ′′ 𝑥 = 𝑙𝑛2 2 3
𝑙𝑛𝑥 𝑥 𝑙𝑛𝑥 𝑥 𝑙𝑛 𝑥

El valor máximo de 𝑓 ′′ 𝑥 (como la función 𝑓 ′′ es decreciente en 2,4 )

ocurre en x=2, entonces 𝑓 ′′ 2 = 1.4014. El error máximo es:
𝑓 ′′ 𝑐 3
1.4014 3
|𝐸| = 𝑏−𝑎 = 4−2 = 0.9343
12 12

Método Trapezoidal compuesto o Regla de los Trapecios.
y
Consideremos la función f positiva, 𝑃2
𝑃1 𝑦=𝑓 𝑥
continua en un intervalo cerrado 𝑎, 𝑏 y
𝑃0 𝑃𝑖−2
𝑃𝑛
una partición regular 𝑃= 𝑃𝑖−1 𝑃𝑖
𝑥0 , 𝑥1 , 𝑥2 , … , 𝑥𝑖−1 , 𝑥𝑖 , … , 𝑥𝑛 en 𝑎, 𝑏 .
Sea A𝑖 el área de un trapecio formado por
los puntos (𝑥𝑖−1 , 0), 𝑥𝑖 , 0 ,
0 a=𝑥0 𝑥1 𝑥2 x
𝑥𝑖−2 𝑥𝑖−1 𝑥𝑖 𝑥𝑛 = 𝑏
𝑃𝑖−1 = (𝑥𝑖−1 , 𝑓 𝑥𝑖−1 y 𝑃𝑖 = (𝑥𝑖 , 𝑓 𝑥𝑖 ) =b
entonces:
Δ𝑥
A𝑖 = 𝑓 𝑥𝑖−1 + 𝑓 𝑥𝑖
2

Del gráfico anterior, el área 𝐴 bajo la curva de f, es aproximadamente igual a la

suma de todos los trapecios formados en la partición, es de decir,
𝑏 𝑛
න 𝑓 𝑥 𝑑𝑥 = 𝐴 ≈ ෍ 𝐴𝑖
𝑎 𝑖=1
𝑏 𝑛−1
Δ𝑥
න 𝑓 𝑥 𝑑𝑥 = 𝐴 ≈ 𝑓 𝑥0 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 𝑥𝑛
𝑎 2
𝑖=1
𝑏 𝑛−1
𝑏−𝑎
න 𝑓 𝑥 𝑑𝑥 ≈ 𝑓 𝑥0 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 𝑥𝑛
𝑎 2𝑛
𝑖=1
A esta aproximación se le conoce como la regla compuesta del trapecio o regla de

los trapecios, cuyo valor es aproximadamente la integral definida de f en 𝑎, 𝑏 .

Estimación del Error en la Regla de los Trapecios
Si f tiene segunda derivada continua en un intervalo cerrado a, b , el error (E)

𝑏
cuando se aproxima ‫ 𝑥𝑑 𝑥 𝑓 𝑎׬‬usando la regla compuesta del trapecio es:
𝑏−𝑎 3
𝐸 ≤ 𝑀, para n subintervalos en 𝑎, 𝑏
12𝑛2
donde M, es un número real positivo tal que 𝑓 ′′ 𝑥 ≤ 𝑀, 𝑥 ∈ 𝑎, 𝑏 .

Ejemplo 1
Dada la integral:
4
1
𝐼 = න 𝑑𝑥; .
1 𝑥
Use la Regla de los Trapecios para:
a) aproximar el valor de la integral 𝐼, para 𝑛 = 6.
b) Estime el error máximo que se comete.
c) Determine 𝑛, que se requiere para aproximar la integral 𝐼, con una exactitud de 10-4.
Solución
a=1, b=4 y n=6

4−1 1
También tenemos ∆𝑥 = = 0.5, 𝑓 𝑥 = , 𝑥 ∈ 1,4 .
6 𝑥
a) Para estimar integral usamos la aproximación siguiente y construimos una

tabla:
4 5
1 4−1
න 𝑑𝑥 ≈ 𝑓 1 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 4
1 𝑥 2(6)
𝑖=1
i 𝒙𝒊 = 𝒂 + 𝒊∆𝒙 𝒇(𝒙𝒊 ) 𝒎 𝒎𝒇(𝒙𝒊 )

0 1.0 1.0000 1 1.0000
1 1.5 0.6667 2 1.3333
2 2.0 0.5000 2 1.0000
3 2.5 0.4000 2 0.8000
4 3.0 0.3333 2 0.6667
5 3.5 0.2857 2 0.5714
6 4.0 0.2500 1 0.2500
Suma Total 5.6214

Reemplazando valores obtenemos:
4
1 0.5
න 𝑑𝑥 ≈ 5.6214 = 1.4054
1 𝑥 2
b) Cálculo del error máximo (E).
1 1 2 2
𝑓 𝑥 = ; 𝑓′ 𝑥 = − ; 𝑓 ′′ 𝑥 = , el valor máximo de 𝑓 ′′ 𝑥 = , (como
𝑥 𝑥2 𝑥3 𝑥3
2
la función es decreciente en 1,4 ) ocurre en x=1, 𝑀 = 𝑓 ′′ 1 = = 2.
13
Entonces, el error (E) máximo es:
𝑀 3
2 3
E= 𝑏−𝑎 = 4−1 = 0.125
12𝑛2 12 6 2

c) Para asegurar el grado de exactitud de 10-4, necesitamos tener:
𝑀 3 −4
E= 𝑏−𝑎 < 10
12𝑛2
2 3 −4
E= 4−1 < 10
12𝑛2
𝑛2 > 450000 ⇒ 𝑛 > 212.13
Entonces el valor de n debe ser mayor o igual a 213.

Ejemplo 2.
Dada la integral
2
න 𝑒 −4𝑥 + cos(𝜋𝑥 )𝑑𝑥
1
Utilice la regla del trapecio compuesto para:
a) Aproximar la integral dada, para 𝑛 = 13
b) Estime el error que se comete en la aproximación
Solución.
a)
2−1
𝑓 𝑥 = 𝑒 −4𝑥 + cos 𝜋𝑥 , 𝑎 = 1, 𝑏 = 2, ∆𝑥 = = 0.07692308
13
2 12
2−1
න 𝑓 𝑥 𝑑𝑥 ≈ 𝑓 1 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 2
1 2(13)
𝑖=1

2
2−1
න 𝑒 −4𝑥 + cos(𝜋𝑥 )𝑑𝑥 ≈ 0.11797326
1 2(13)
2
න 𝑒 −4𝑥 + cos(𝜋𝑥 )𝑑𝑥 ≈ 0.0045374
1
b)
𝑓 ′ 𝑥 = −4𝑒 −4𝑥 − 𝜋𝑠𝑒𝑛(𝜋𝑥)
𝑓 ′ ′ 𝑥 = 16𝑒 −4𝑥 − 𝜋 2 𝑐𝑜𝑠(𝜋𝑥)
Para 𝑥 ∈ 1,2
𝑓′′ 𝑥 = 16𝑒 −4𝑥 − 𝜋 2 𝑐𝑜𝑠(𝜋𝑥) ≤ 16𝑒 −4𝑥 + 𝜋 2 cos 𝜋𝑥 ≤ 16𝑒 −4 1 + 𝜋2 1 =𝑀
2−1 3 10.1627
𝐸 ≤ 2 16𝑒 −4 1 + 𝜋2 1 =
12 13 12 13 2
𝐸 ≤ 0.0050112

2.3 Método de Simpson y método de Simpson compuesto. Estimación
del error.
Regla o Método de Simpson
Dada la función f positiva, continua en un intervalo cerrado 𝑎, 𝑏 y una partición

regular 𝑃 = 𝑥0 , 𝑥1 , 𝑥2 , … , 𝑥𝑖−1 , 𝑥𝑖 , … , 𝑥𝑛 en 𝑎, 𝑏 con n par.
𝑃1
y 𝑦=𝑓 𝑥
Este método aproxima en cada sub
𝑃2
intervalo la función f mediante un 𝑃0 𝑃 𝑥 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐
polinomio de segundo grado

𝑃 𝑥 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐, 𝑎, 𝑏, 𝑐 ∈ ℝ ,
como se muestra en la figura
0 𝑥0 𝑥1 𝑥2 x
𝑥2 −𝑥0
En este caso se tiene, ∆𝑥 = ℎ = , 𝑥1 = 𝑥0 + ℎ, 𝑥2 = 𝑥0 + 2ℎ, la integral definida
2
en 𝑥0 , 𝑥2 se aproxima de la siguiente forma:

2.3 Método de Simpson y método de Simpson compuesto. Estimación del error.
𝑥2 𝑥2
න 𝑓(𝑥)𝑑𝑥 ≈ න 𝑃 𝑥 𝑑𝑥
𝑥0 𝑥0
Sean 𝑃0 = 𝑥0 , 𝑦0 , 𝑃1 = 𝑥1 , 𝑦1 y 𝑃2 = 𝑥2 , 𝑦2 los puntos de intersección

de la curva f con el polinomio, entonces se cumple:
𝑦0 = 𝑎𝑥02 + 𝑏𝑥0 + 𝑐 … (1)

𝑦1 = 𝑎𝑥12 + 𝑏𝑥1 + 𝑐 = 𝑎 𝑥0 + ℎ 2 + 𝑏(𝑥0 +ℎ) + 𝑐 … (2)
𝑦2 = 𝑎𝑥22 + 𝑏𝑥2 + 𝑐 = 𝑎 𝑥0 + 2ℎ 2 + 𝑏(𝑥0 + 2ℎ) + 𝑐 … (3)
Integrando el polinomio en el intervalo 𝑥0 , 𝑥2 ,
𝑥2 𝑥2 𝑥0 +2ℎ
න 𝑓(𝑥)𝑑𝑥 ≈ න 𝑃 𝑥 𝑑𝑥 = න 𝑎𝑥 2 + 𝑏𝑥 + 𝑐 𝑑𝑥
𝑥0 𝑥0 𝑥0

𝑥0 +2ℎ
𝑥0+2ℎ
2
𝑥3 𝑥2
න 𝑎𝑥 + 𝑏𝑥 + 𝑐 𝑑𝑥 = 𝑎 + 𝑏 + 𝑐𝑥อ
𝑥0 3 2
𝑥0
Evaluando y simplificando obtenemos
𝑥2 𝑥2
ℎ
න 𝑓(𝑥)𝑑𝑥 ≈ න 𝑃 𝑥 𝑑𝑥 = 𝑦 + 4𝑦1 + 𝑦2 .
𝑥0 𝑥0 3 0
𝑏−𝑎
Si f es una función positiva continua en un intervalo 𝑎, 𝑏 , ∆𝑥 = ℎ = ,
2
la integral definida por la regla de Simpson está dada por la siguiente

aproximación:
𝑏
𝑏−𝑎 𝑏−𝑎 𝑎+𝑏
න 𝑓 𝑥 𝑑𝑥 ≈ 𝑓 𝑥0 + 4𝑓 𝑥1 + 𝑓 𝑥2 = 𝑓 𝑎 + 4𝑓 +𝑓 𝑏
𝑎 6 6 2

Ejemplo
Usando el método de Simpson, calcule el valor aproximado de la siguiente

integral:
3
6
න 𝑑𝑥
1 𝑥2
Solución
Tenemos:
6 3−1
𝑓 𝑥 = 2 ; 𝑎 = 1, 𝑏 = 3 → ∆𝑥 = = 1; 𝑥1 = 𝑎 + ∆𝑥 = 2.
𝑥 2

Aplicando el método de Simpson, se tiene:
3
6 3−1
න 𝑑𝑥 ≈ 𝑓 𝑥0 + 4𝑓 𝑥1 + 𝑓 𝑥2
1 𝑥2 6
𝑓 𝑥0 = 𝑓 1 = 6.00
𝑓 𝑥1 = 𝑓 2 = 1.50
𝑓 𝑥2 = 𝑓 3 = 0.667
Reemplazando valores se tiene,
3
6 1 1
න 2 𝑑𝑥 ≈ 6.00 + 4 1.50 + 0.667 = 12.667 = 4.222
1 𝑥 3 3

Regla compuesta de Simpson (regla 1/3)
Sea f una función positiva continua en un intervalo cerrado 𝑎, 𝑏 y 𝑃 =

𝑥0 , 𝑥1 , 𝑥2 , … , 𝑥𝑛−1 , 𝑥𝑛 una partición regular de 𝑎, 𝑏 con n un número entero
par; entonces la integral definida por la regla de Simpson dada por:
𝑛 𝑛
−1
𝑏 2 2
𝑏−𝑎
න 𝑓 𝑥 𝑑𝑥 ≈ 𝑓 𝑥𝑜 + 2 ෍ 𝑓 𝑥2𝑗 + 4 ෍ 𝑓 𝑥2𝑗−1 + 𝑓 𝑥𝑛
𝑎 3𝑛
𝑗=1 𝑗=1
La regla compuesta de Simpson es conocida como la regla 1/3 de Simpson.

Error máximo al estimar la integral por el método compuesto de Simpson:
𝑏
El error 𝐸 por el cálculo de la integral ‫ 𝑥𝑑 𝑥 𝑓 𝑎׬‬para n subintervalos es:
𝑏 𝑏
𝑏−𝑎 5 4
𝐸 = න 𝑓 𝑥 𝑑𝑥 − න 𝑃 𝑥 𝑑𝑥 ≤ 4 𝑚á𝑥 𝑓 𝑐 , 𝑐 ∈ 𝑎, 𝑏
𝑎 𝑎 180𝑛
Ejemplo 1
Del ejemplo de la Regla del Trapecio, ahora se pide aproximar la integral y
calcular el error, usando el método de Simpson para a) n = 2; b) n=6 y c) calcular
el valor de n para que la integral se aproxime con una exactitud de 10-4.

Solución
41 1 𝑏−𝑎
‫׬‬1 𝑥 𝑑𝑥: 𝑓(𝑥) = 𝑥 ; 𝑎 = 1; 𝑏 = 4 → 𝑥𝑖 = 𝑎 + 𝑖∆𝑥; ∆𝑥 = 𝑛
a) Para n=2; ∆𝑥 =1.5

𝑓 𝑥0 = 𝑓 1 = 1.00
𝑓 𝑥1 = 𝑓 2.5 = 0.40
𝑓 𝑥2 = 𝑓 4 = 0.25
Reemplazando valores en la regla de Simpson, la integral es aproximadamente:
4
1 4−1
න 𝑑𝑥 ≈ 𝑓 1 + 4𝑓 2.5 + 𝑓 4
1 𝑥 3(2)
4
1 1
න 𝑑𝑥 ≈ 1 + 4 0.4 + 0.25 = 1.425
1 𝑥 2

Para el cálculo, tenemos

1 ′ 1 24
𝑓 𝑥 = ; 𝑓 𝑥 = − 2 ;… , 𝑓 (4) 𝑥 =
𝑥 𝑥 𝑥5
24
El máximo ocurre en x=1 , entonces 𝑓 (4) 1 = = 24.
13
Entonces el error máximo es:
24 5
𝐸≤ 4 3 = 2.025
180 2

b) Para n=6, la integral es:

4
1 4−1
න 𝑑𝑥 ≈ 𝑓 1 + 4𝑓 𝑥1 + 2𝑓 𝑥2 + 4𝑓 𝑥3 + 2𝑓 𝑥4 + 4𝑓 𝑥5 + 𝑓 𝑥6
1 𝑥 3(6)
Construimos una tabla para estimar la integral:
i 𝒙𝒊 = 𝟏 + 𝒊∆𝒙 𝒇 𝒙𝒊 𝒎 𝒎𝒇 𝒙𝒊
0 1.0 1.0000 1 1.0000
1 1.5 0.6667 4 2.6667
2 2.0 0.5000 2 1.0000
3 2.5 0.4000 4 1.6000
4 3.0 0.3333 2 0.6667
5 3.5 0.2857 4 1.1429
6 4.0 0.2500 1 0.2500
suma 8.3262

Reemplazando valores, la integral aproximada es;
4
1 1
න 𝑑𝑥 ≈ 8.3262 = 1.39
1 𝑥 6
El cálculo del error máximo E, es como se realizó en la parte a)

24
𝑓 (4) 1 = = 24. Entonces el error máximo es:
13
24 5
𝐸≤ 4 3 = 0.025
180 6

Como se puede observar que el error encontrado para n=6 es menor al

estimado para n=2, entonces concluimos que cuando n es mayor, el error
es menor.
c) Para asegurar el grado de exactitud de 10-4, resolvemos:
𝑏−𝑎 5 4 −4
𝑚á𝑥 𝑓 𝑐 < 10
180𝑛4
4−1 5 −4
24 < 10
180𝑛4
𝑛4 > 324000 ⇒ 𝑛 = 23.86

Entonces el valor de n es mayor o igual que 24.
Se puede observar que para obtener la integral con la misma exactitud el número
de subintervalos, por el método de Simpson compuesto (n=24) es mucho menor
que para el método de los trapecios (n=213),
Ejemplo 2. Dada la integral

−3
𝑥6
න − 𝑒 −𝑥 𝑑𝑥
−5 180
Calcule el número de sub intervalos (𝑛), de manera que, al aproximar la integral

por la regla compuesta de Simpson, con una exactitud de 9 × 10−3.
Solución
𝑥6
Hallemos las derivadas de la función 𝑓 𝑥 = − 𝑒 −𝑥
180

𝑥5
𝑓′ 𝑥 = + 𝑒 −𝑥 ,…, 𝑓 ′𝑣 𝑥 = 2𝑥 2 − 𝑒 −𝑥 y 𝑓 𝑣 𝑥 = 4𝑥 + 𝑒 −𝑥 .
30
Como 𝑓 𝑣 𝑥 > 0 ⟹ 𝑓 ′𝑣 𝑥 es creciente para todo 𝑥 ∈ −5, −3 además

𝑓 ′𝑣 −5 = 50 − 𝑒 5 y 𝑓 ′𝑣 −3 = 18 − 𝑒 3,
por tanto, el máx 𝑓 ′𝑣 𝑥 = 𝑒 5 − 50 = 98.4131591 (como 𝑓 𝑣 𝑥 = 4𝑥 + 𝑒 −𝑥

y 0 < 𝑒 3 − 12 ≤ 4𝑥 + 𝑒 −𝑥 ≤ 𝑒 5 − 20, se tiene 𝑓 𝑣 > 0 entonces, 𝑓 𝑖𝑣 es
creciente en −5, −3 ). Luego por la fórmula de error máximo se tiene:
−3 + 5 5
𝑒 5 − 50 < 9 × 10−3 ,
180𝑛4
resolviendo 𝑛4 > 1943.963637 ⇒ 𝑛 > 6.640060466
De este modo,
𝒏≥𝟖

Ejemplo 3
Un sólido de revolución se obtiene al girar la región limitada por la curva

𝑥2
𝑦=𝑓 𝑥 =1+ ; y las rectas 𝑥 = 0; 𝑥 = 3; 𝑦 = 0, alrededor de eje X. Si el
9
volumen generado del sólido es:

3
𝑉 = 𝜋 න 𝑓 2 (𝑥) 𝑑𝑥
0
Calcule dicho volumen utilizando integración numérica; para n=10:
a) Por el método de los trapecios y el error por la aproximación.

b) Por el método compuesto de Simpson y el error por la aproximación.
Solución.

2
𝑥2
Como 𝑓 𝑥 = 1 + , el volumen generado (V) es:
9
3 3 2
𝑥2
𝑉 = 𝜋 න 𝑓 2 (𝑥)𝑑𝑥 = 𝜋 න 1+ 𝑑𝑥
0 0 9
3−0
Tenemos: 𝑎 = 0; 𝑏 = 3; 𝑛 = 10 → Δ𝑥 = = 0.30
(10)
Con esta información construimos la siguiente tabla, que sirve para aproximar la
integral para los dos métodos:

i 𝒙𝒊 = 𝟎 + 𝒊∆𝒙 𝑓 𝑥𝑖 𝑚 𝑚𝑓 𝑥𝑖 𝑚 𝑚𝑓 𝑥𝑖
Método del Trapecio Método de Simpson

0 0.0 1.0000 1 1.0000 1 1.0000
1 0.3 1.0201 2 2.0402 4 4.0804
2 0.6 1.0816 2 2.1632 2 2.1632
3 0.9 1.1881 2 2.3762 4 4.7524
4 1.2 1.3456 2 2.6912 2 2.6912
5 1.5 1.5625 2 3.1250 4 6.2500
6 1.8 1.8496 2 3.6992 2 3.6992
7 2.1 2.2201 2 4.4402 4 8.8804
8 2.4 2.6896 2 5.3792 2 5.3792
9 2.7 3.2761 2 6.5522 4 13.1044
10 3.0 4.0000 1 4.0000 1 4.0000
Suma Total 37.4666 56.00

Reemplazando valores para cada método:
a) Aproximación por el método de los trapecios.
3 2 9
𝑥2 3−0
𝑉 = 𝜋න 1+ 𝑑𝑥 ≈ 𝜋 𝑓 0 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 3
0 9 2 10
𝑖=1
3 2
𝑥2
𝑉 = 𝜋න 1+ 𝑑𝑥 ≈ 𝜋(0.15) 37.4666 = 5.62𝜋 𝑢3
0 9

Error máximo por la aproximación:
4 4 2 16
Como 𝑓 ′′ 𝑥 = + 𝑥 , entonces 𝑀 = 𝑓 ′′ 3 = , el error máximo es:
9 27 9
𝑀 3
16/9 3
𝐸= 𝑏−𝑎 = 3−0 = 0.04
12𝑛2 12 10 2
b) Aproximación por el Método de Simpson.
3 2
𝑥2
𝑉 = 𝜋න 1+ 𝑑𝑥
0 9
3−0
𝑉≈ 𝑓 0 + 4𝑓 𝑥1 + 2𝑓 𝑥2 + ⋯ + 2𝑓 𝑥8 + 4𝑓 𝑥9 + 𝑓 3 𝜋
3 10

Reemplazando los valores de la tabla, se tiene:
3 2
𝑥2
𝑉 = 𝜋න 1+ 𝑑𝑥 ≈ 0.1 56.00 𝜋 = 5.6𝜋 𝑢3
0 9
Error máximo para esta aproximación:
8 8
Como 𝑓 ′𝑣 𝑥 = , se tiene 𝑓 (4) 3 = = 0.2963. Entonces el error máximo
27 27
es:
0.26963 5
𝐸≤ 3 = 4.0𝑥10−5
180 10 4

Ejemplo 4. Si 𝑓 es una función, tal que −6 ≤ 𝑓 4 𝑥 ≤ 6, 2 ≤ 𝑥 ≤
48.4
48.4. Si se quiere aproximar la integral ‫׬‬2 𝑓 𝑥 𝑑𝑥 con la regla
compuesta de Simpson con una exactitud de 10−1 . Determine el número
de subintervalos (𝑛)
Solución.
Para asegurar el grado de exactitud de 10-1, resolvemos
𝑏−𝑎 5 4 −1
𝑚á𝑥 𝑓 𝑐 < 10
180𝑛4
48.4 − 2 5
𝑚á𝑥 𝑓 4 𝑐 =6⇒ 6 < 10 −1
180𝑛4
𝑛4 > 71691661.91
𝑛 = 92.02 ⇒ 𝑛 = 94

Ejemplo 5. Entre las 4:00 p.m. y 6:00 p.m. se midió el caudal Q (en litros por
minuto) que sale de un depósito de agua, para uso doméstico de un determinado
pueblo, cuyo resultado se muestra en la siguiente tabla.
T (tiempo) Q (litros/minuto) T (tiempo) Q (litros/minuto)
04:00 17 05:12 29
04:12 19 05:24 38
04:24 18 05:36 43
04:36 21 05:48 41
04:48 28 06:00 39
05:00 30
i) Qué volumen de agua fue utilizada en este período de 2 horas (use la regla
compuesta de Simpson);
ii) cuál es caudal medio en dicho intervalo de tiempo.
Solución.

120
Tenemos: 𝑎 = 0; 𝑏 = 120; 𝑛 = 10 → Δt = = 12 y para Simpson
10
compuesto aproximamos el volumen
120
120
𝑉𝑜𝑙 = න 𝑄 𝑡 𝑑𝑡 ≈ 𝑄 0 + 4𝑄 12 + 2𝑄 24 + ⋯ + 4𝑄 108 + 𝑄 120
0 3 10
𝒊 𝒕𝒊 𝑸(𝒕𝒊 ) m 𝒎𝑸(𝒕𝒊 )
0 0 17 1 17
1 12 19 4 76
2 24 18 2 36
3 36 21 4 84
4 48 28 2 56
5 60 30 4 120
6 72 29 2 58
7 84 38 4 152
8 96 43 2 86
9 108 41 4 164
10 120 39 1 39
suma 888

Con los datos y la información de la tabla calculamos:
i) Volumen de agua:
120
120 − 0
𝑉𝑜𝑙 = න 𝑄 𝑡 𝑑𝑡 ≈ 888 = 3552
0 3 10
El volumen de agua utilizada es de 3552 litros.
ii) En base al teorema del valor medio para las integrales definidas, obtenemos el
caudal medio 𝑄ത :
1 120
𝑄ത = ‫׬‬ 𝑄 𝑡 𝑑𝑡 ≈ 29.60
120−0 0
litros
El caudal medio para uso doméstico es aproximadamente de 29.60 .
min

Ejemplo 7. En Lima, el día lunes 27 de diciembre, entre las 02:00 y las
23:00 horas, se registraron las rachas máximas de viento (en kilómetros por
hora), el registro se muestra en la siguiente tabla:
Utilice la regla compuesta de Simpson para determinar el valor promedio

aproximado de las rachas de viento en este periodo.

Solución.
Valor promedio (𝑉𝑃):
21
1
𝑉𝑃 = න 𝑉 𝑡 𝑑𝑡
21 − 0 0
Aproximación de la integral: 𝑎 = 0, 𝑏 = 21, 𝑛 = 7; ∆𝑥 = 3
Como 𝑛 es impar y pide el Método de Simpson se debe calcular hasta (𝑛 − 1)

sub intervalos con Simpson y el último intervalo por el método del trapecio.
21 18 21
න 𝑉 𝑡 𝑑𝑡 = න 𝑉 𝑡 𝑑𝑡 + න 𝑉 𝑡 𝑑𝑡
0 0 18

Para Simpson
18
18 − 0
න 𝑉 𝑡 𝑑𝑡 ≈ 392 = 392
0 3 6

Trapecio
𝑉 18 = 24; 𝑉 21 = 23
21
3
න 𝑉 𝑡 𝑑𝑡 ≈ 24 + 23 = 70.5
18 2
Reemplazando valores la integral se aproxima según:

21
න 𝑉 𝑡 𝑑𝑡 ≈ 392 + 70.5 = 462.5
0
21
1 1
𝑉𝑃 = න 𝑉 𝑡 𝑑𝑡 = (462.5)
21 − 0 0 21 − 0
𝑉𝑃 = 22.02
En valor promedio es de 22.02 𝑘𝑚/ℎ


Integración Numérica 2022-2

Cargado por

Copyright:

Formatos disponibles

Integración Numérica 2022-2

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Integración Numérica 2022-2

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional Agraria

CURSO : CÁLCULO PARA CIENCIA S

TEMA: INTEGRACIÓN NUMÉRICA

2.1 Definiciones básicas para Integración Numérica

Sea 𝐼 = 𝑎, 𝑏 un intervalo cerrado, llamaremos una partición de 𝐼 a un conjunto de puntos

Se obtiene n sub intervalos de la partición P, 𝑎, 𝑥1 , 𝑥1 , 𝑥2 , 𝑥2 , 𝑥3 , … , 𝑥𝑛 −1 , 𝑏 ; en

La longitud de cada sub intervalo es, 𝑥𝑖 − 𝑥𝑖−1 = ∆𝑥𝑖 y la sumatoria es

Aplicando el teorema anterior, determine en qué intervalo se encuentra

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Aplicando el teorema anterior obtenemos:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Sea f una función continua y en un intervalo cerrado 𝑎, 𝑏 . Entonces, existe 𝑐 ∈ 𝑎, 𝑏

Al valor f(c) se le denomina “valor medio” de la función f en 𝑎, 𝑏 .

Para la interpretación geométrica se requiere además que 𝑓 𝑥 > 0.

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

máximo absoluto 𝑀 = 𝑓 10 = 𝑙𝑛 10 y mínimo absoluto 𝑚 = 𝑓 3 = 𝑙𝑛 3 .

Entonces el intervalo donde se encuentra el valor de la integral es 7.69,16.12

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Usando el Teorema del valor medio de la integral definida, se tiene:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Ejemplo 3. El crecimiento de una determinada planta en los primeros 30 días

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

El tamaño promedio de la planta en los 30 días es de 𝑦ത = 23.06 𝑐𝑚, como se

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

resolviendo se tiene 𝑐 = 11.81 días, entonces el día aproximado donde la planta

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

න 𝑓 (𝑥 )𝑑𝑥 ≈ ෍ 𝑓 (𝑤𝑖 )Δ𝑥

donde 𝑤𝑖 es un número arbitrario en el i-ésimo subintervalo 𝑥𝑖−1 , 𝑥𝑖 de la partición

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Si 𝑤𝑖 = 𝑥𝑖 calculamos f en el extremo derecho del sub intervalo 𝑥𝑖−1 , 𝑥𝑖 , entonces la

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

a) Aproximamos la integral evaluando f en el

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

𝑖−1 𝑥𝑖−1 𝑓 𝑥𝑖−1

Por lo tanto la integral aproximada por la izquierda es:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

En la siguiente tabla se muestra los valores de f para cada 𝑥𝑖 y la suma total.

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Una buena aproximación es:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Con los valores de la tabla obtenemos:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

a) Aproximamos la integral evaluando f en el extremo izquierdo:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

𝑖−1 𝑥𝑖−1 𝑓 𝑥𝑖−1

Por lo tanto la integral aproximada por la izquierda es:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Una mejor aproximación es:

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Si 𝑓 tiene segunda derivada continua en un intervalo cerrado a, b , el error (E)

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

Cálculo del error máximo (E).

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL

El valor máximo de 𝑓 ′′ 𝑥 (como la función 𝑓 ′′ es decreciente en 2,4 )

INTEGRACIÓN NUMÉRICA UNALM 2022 - II PRESENCIAL - VIRTUAL