Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
Actividad 2: Laboratorio: Diferencia de población
Objetivos
En esta actividad pondremos en práctica diferentes técnicas para aproximar
numéricamente las integrales.
Pautas de elaboración
Queremos analizar el crecimiento de un tipo de bacterias que estamos estudiando.
Sea 𝐵(𝑡) la población de bacterias que varía a lo largo del tiempo. Si conocemos la
velocidad a la que la población varía (𝐵′ (𝑡)), podemos hallar la diferencia de
población entre dos instantes aplicando la regla de Barrow
𝑏
𝐵(𝑏) − 𝐵(𝑎) = ∫ 𝐵′ (𝑡)𝑑𝑡
𝑎
La velocidad de variación de tipo de bacterias que estamos estudiando se describe
como:
′ (𝑡)
30 𝑒 −0.1𝑡
𝐵 =
(1 + 3𝑒 −0.1𝑡 )2
Donde 𝑡 es el tiempo en horas.
Calcular la variación de la población de bacterias entre t=0 y t=20: obtener el
resultado exacto y el resultado aproximado utilizando la regla del trapecio y la
regla de Simpson con 2 subintervalos. Compare los resultados obtenidos.
Solución:
Para poder implementar la regla de trapecios y de Simpson con dos subintervalos, se
parte de la expresión general para valores pares de 𝑛 que denota los puntos en un
𝑥𝑛 −𝑥0 𝑏−𝑎
intervalo de análisis [𝑥0 , 𝑥𝑛 ] con ℎ = = . Si 𝑓 𝜖 𝐶 𝑛+2 [𝑎, 𝑏] se cumple:
𝑛 𝑛
1
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
𝑏 𝑛 𝑛
ℎ𝑛+3
∫ 𝑓(𝑥)𝑑𝑥 = ∑ 𝑎𝑖 𝑓(𝑥𝑖 ) + 𝑓 𝑛+2 (𝜉) ∫ 𝜏 2 (𝜏 − 1) ∙∙∙ (𝜏 − 𝑛)𝑑𝜏
(𝑛 + 2)!
𝑎 𝑖=0 0
𝑛+1
Sin embargo, para valores impares de 𝑛 y si 𝑓 𝜖 𝐶 [𝑎, 𝑏] se cumple que:
𝑏 𝑛 𝑛
ℎ𝑛+2
∫ 𝑓(𝑥)𝑑𝑥 = ∑ 𝑎𝑖 𝑓(𝑥𝑖 ) + 𝑓 𝑛+1 (𝜉) ∫ 𝜏(𝜏 − 1) ∙∙∙ (𝜏 − 𝑛)𝑑𝜏
(𝑛 + 1)!
𝑎 𝑖=0 0
Para la selección de dos subintervalos entre [0,20] se puede dividirlos como se
muestra a continuación:
Es decir que 𝑥0 = 0, 𝑥1 = 10, 𝑥2 = 20 para utilizar el método de Simpson, mientras
que para utilizar el método de Trapecios la selección 𝑥0 = 0, 𝑥1 = 20 es la que servirá
para aproximar la diferencia de población entre dos instantes de tiempo aplicando la
regla de Barrow.
Para el método de Trapecios, partimos de la fórmula general:
𝑛−1
ℎ
𝐼 ≈ ∑(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 ))
2
𝑖=0
ℎ
= (𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑥𝑛 ))
2
Con el error:
ℎ2
𝜀=− (𝑏 − 𝑎)𝑓"(𝜉(𝑥))
12
2
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
Particularizando para solamente dos subintervalos, es decir, 𝑛 = 2, se obtiene la
siguiente expresión:
2
ℎ ℎ
𝐼 ≈ ∑(𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑖+1 )) = (𝑓(𝑥0 ) + 2𝑓(𝑥1 ) + 𝑓(𝑥2 ))
2 2
𝑖=0
Para el método de Simpson con 𝑛 = 2 se obtiene la siguiente expresión:
ℎ 𝑎+𝑏
𝐼≈ (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏))
3 2
𝑎+𝑏
En donde 𝑥0 = 𝑎, 𝑥1 = , 𝑥2 = 𝑏, con su respectivo error definido por:
2
ℎ5 𝑖𝑣
𝜀=− 𝑓 (𝜉(𝑥))
90
Implementando el código proporcionado en las sesiones exclusivas sobre el tema en
particular para aproximar una integral por método de Simpson y Trapecios, se
obtuvieron los siguientes resultados:
𝟐𝟎
30 𝑒 −0.1𝑡
∫ = 𝟒𝟓. 𝟒𝟒𝟖
(1 + 3𝑒 −0.1𝑡 )2
𝟎
𝑛 Simpson Trapecio Error Simpson 𝐸𝑠 Error Trapecio 𝐸𝑡
𝟐 46.348 44.583 0.22497 1.5401
Tabla 1 . Tabla de resultados.
Como se puede observar en la figura 1, existe similitud entre el valor obtenido de la
integral utilizando el método de Simpson con el valor de la integral exacta entre sus
puntos de interés, sin embargo, con el método de Trapecios existe una variación
considerable de 1.5401 unidades de crecimiento del tipo de bacteria. Cabe recalcar
que el valor exacto respecto al crecimiento de bacterias entre 0 ℎ𝑜𝑟𝑎𝑠 y 20 ℎ𝑜𝑟𝑎𝑠 es
de 46.123 unidades de crecimiento de bacteria.
3
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
Calcular la variación de la población de bacterias entre t=0 y t=20: obtener el
resultado aproximado utilizando la regla del trapecio y la regla de Simpson con
3,4,5,6,7,8,9 y 10 subintervalos siempre que sea posible. Compare los resultados
obtenidos.
Solución:
Para obtener los distintos valores del crecimiento de este tipo de bacteria para
una serie de intervalos, se puede utilizar de igual manera el mismo código
mostrado en la figura 1, de tal forma que se puede variar 𝑛1 = 3,4,5,6,7,8,9,10.
𝟐𝟎
30 𝑒 −0.1𝑡
∫ = 𝟒𝟓. 𝟒𝟒𝟖
(1 + 3𝑒 −0.1𝑡 )2
𝟎
𝑛 Simpson Trapecio Error Simpson 𝐸𝑠 Error Trapecio 𝐸𝑡
𝟑 - 45.4480 - 0.6754
𝟒 46.133 45.745 0.0090788 0.37823
𝟓 - 45.8819 - 0.2416
𝟔 46.125 45.956 0.0017275 0.16757
𝟕 - 46.0004 - 0.1230
𝟖 46.124 46.124 0.00053927 0.094152
𝟗 - 46.0491 - 0.07440
𝟏𝟎 46.124 46.063 0.00021952 0.060227
Tabla 2. Tabla de resultados
Como se puede observar en la tabla 2, no hay resultados para el método de
Simpson con un número de intervalos impares. Esto sucede debido a que la
fórmula de Simpson utiliza tres puntos de interés, por ende, el número de
intervalos para poder utilizar el método debe ser par. En la siguiente figura se
ilustra como se distribuyen los intervalos para número par e impar.
4
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
Figura 1. Distribución de intervalos
5
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
Si la población de bacterias inicial es de 25 (𝑩(𝟎) = 𝟐𝟓) ¿cuál es el tamaño de la
población al cabo de 20 horas?
Solución:
Lo que se necesita básicamente es encontrar el valor del tamaño de la población
exactamente cuando 𝑡 = 20 es decir:
20
𝐵(20) = 25 + ∫ 𝐵′ (𝑡)𝑑𝑡
0
20
𝐵(20) = 25 + ∫ 𝐵′ (𝑡)𝑑𝑡 = 25 + 45.448
0
𝐵(20) = 70.448
Es decir que el tamaño de la población de este tipo de bacterias exactamente después
de 20 ℎ𝑜𝑟𝑎𝑠 es 70.448 unidades.
6
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
ANEXOS
7
Actividades
Asignatura Datos del alumno Fecha
Métodos Numéricos Apellidos: AVILA AGUILAR
14/01/2025
Aplicados I Nombre: WILLIAM ALEXANDER
8
Actividades
