1
UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
Universidad del Perú, Decana de América
FACULTAD DE QUÍMICA E INGENIERÍA QUÍMICA
Escuela Profesional de Ingeniería Química
“Métodos numéricos”
Informe Individual
“Ejercicio:”
Alumno:
Ramos Cano Caín Fernando (22070040)
Docente:
Jose Orlando Calvay Castillo
Lima - Perú
2025
2
TABLA DE CONTENIDO
MARCO TEÓRICO......................................................................................................................3
a) Dióxido de carbono...............................................................................................................3
b) Método de Euler...................................................................................................................3
c) Método de Heun...................................................................................................................3
PLANTEAMIENTO DEL CASO................................................................................................4
PSEUDOCÓDIGO........................................................................................................................5
a) Método de Euler...................................................................................................................5
b) Método de Heun...................................................................................................................5
DIAGRAMA DE FLUJO..............................................................................................................6
a) Método de Euler...................................................................................................................6
b) Método de Heun...................................................................................................................6
ALGORITMO................................................................................................................................7
a) Método de Euler...................................................................................................................7
b) Método de Heun...................................................................................................................9
INTERPRETACIÓN DE RESULTADOS................................................................................11
BIBLIOGRAFÍA.........................................................................................................................11
3
MARCO TEÓRICO
a) Dióxido de carbono
El dióxido de carbono (CO2) es un importante gas que atrapa el calor, o gas de efecto
invernadero, que proviene de la extracción y quema de combustibles fósiles (como
carbón, petróleo y gas natural), de incendios forestales y de procesos naturales como
erupciones volcánicas.
b) Método de Euler
Según de la Universidad Nacional Mar de la Plata (UNMdP), el método de Euler
consiste en encontrar iterativamente la solución de una ecuación diferencial de
primer orden y valores iniciales conocidos para un rango de valores. Partiendo de un
valor inicial x0 y avanzando con un paso h, se pueden obtener los valores de la
solución de la siguiente manera:
Donde Y es solución de la ecuación diferencial y f es la ecuación diferencial en
función de las variables independientes.
c) Método de Heun
Según de la Universidad Nacional Mar de la Plata (UNMdP), el método de Heun
consiste en una mejora del método de Euler para resolver ecuaciones diferenciales de
primer orden y conocido el valor inicial. En este caso, lo que se realiza es un
promedio entre el valor obtenido por Euler y otro obtenido a partir de la
aproximación del valor de la función en el punto siguiente, también por Euler.
Donde Y es solución de la ecuación diferencial, f es la ecuación diferencial en
función de las variables independientes y la solución de Y ˙ k+1 es una aproximación
de Euler.
4
PLANTEAMIENTO DEL CASO
El sitio web de la NASA, Worldview, ofrece una plataforma interactiva que permite
visualizar datos satelitales en tiempo casi real sobre diversos parámetros atmosféricos,
incluyendo la concentración de dióxido de carbono (CO₂). El CO₂ es un gas
contaminante emitido principalmente por fuentes industriales y el tráfico vehicular. Su
acumulación en la atmósfera contribuye al efecto invernadero y al cambio climático.
En Perú, las concentraciones de CO₂ varían según la actividad industrial y las
condiciones meteorológicas, lo que hace necesario modelar su evolución en el tiempo.
Para ello, se emplea un modelo matemático que describe la concentración de CO₂ (C)
en función del tiempo (t), considerando la tasa de emisión (E), el coeficiente de
dispersión atmosférica (D) y la tasa de absorción por sumideros naturales como los
océanos y la vegetación (R):
dC
=E−DC −RC
dt
Dónde:
E= 8ug/m3 /día es la tasa de emisión
D= 0.15 es el coeficiente de dispersión atmosférica
R= 0.07 es la tasa de absorción del CO2 debido a procesos naturales
C0= 80 ug/m3 es la concentración inicial
El modelo matemático permite analizar la evolución de la concentración de CO₂ a lo largo de un
período de 30 días. Se busca determinar cómo varía la concentración con el tiempo y si alcanza
un equilibrio.
5
PSEUDOCÓDIGO
a) Método de Euler
- Inicio
- Definir los parámetros del problema
- Asignar los valores iniciales
- Resolver por el método de Euler
- Mostrar resultados
- Fin
b) Método de Heun
- Inicio
- Definir los parámetros del problema
- Asignar los valores iniciales
- Resolver por el método de Heun
- Mostrar resultados
- Fin
DIAGRAMA DE FLUJO
a) Método de Euler
6
b) Método de Heun
ALGORITMO
a) Método de Euler
7
// DEFINIR PARÁMETROS DEL PROBLEMA
E = 8; // Tasa de emisión (µg/m³/día)
D = 0.15; // Coeficiente de dispersión
R = 0.07; // Tasa de absorción
C0 = 80; // Concentración inicial (µg/m³)
t_final = 30; // Tiempo total en días
h = 1; // Paso de tiempo en días
k = D + R; // Parámetro combinado
N = t_final / h + 1; // Número de pasos
// INICIALIZAR VECTORES PARA ALMACENAR RESULTADOS
C_euler = zeros(N, 1);
t = (0:h:t_final)'; // Vector de tiempo
// ASIGNAR CONDICIÓN INICIAL
C_euler(1) = C0;
// MÉTODO DE EULER
for i = 2:N
dC = E - k * C_euler(i-1); // Evaluar derivada en el punto actual
C_euler(i) = C_euler(i-1) + h * dC; // Actualizar concentración
end
// MOSTRAR RESULTADOS EN CONSOLA
disp("Día | Euler (µg/m³)");
for i = 1:N
disp([t(i), C_euler(i)]);
end
// GRAFICAR RESULTADOS
clf;
plot(t, C_euler, 'r-o', "LineWidth", 2);
xlabel("Tiempo (días)");
ylabel("Concentración de CO2 (µg/m³)");
title("Evolución de la Concentración de CO2 (Método de Euler)");
legend("Euler", "Location", "best");
xgrid();
8
b) Método de Heun
9
// DEFINIR PARÁMETROS DEL PROBLEMA
E = 8; // Tasa de emisión (µg/m³/día)
D = 0.15; // Coeficiente de dispersión
R = 0.07; // Tasa de absorción
C0 = 80; // Concentración inicial (µg/m³)
t_final = 30; // Tiempo total en días
h = 1; // Paso de tiempo en días
k = D + R; // Parámetro combinado
N = t_final / h + 1; // Número de pasos
// INICIALIZAR VECTORES PARA ALMACENAR RESULTADOS
C_heun = zeros(N, 1);
t = (0:h:t_final)'; // Vector de tiempo
// ASIGNAR CONDICIÓN INICIAL
C_heun(1) = C0;
// MÉTODO DE HEUN
for i = 2:N
dC_initial = E - k * C_heun(i-1); // Evaluar derivada en el punto actual
C_pred = C_heun(i-1) + h * dC_initial; // Paso predictor (Euler)
dC_pred = E - k * C_pred; // Evaluar derivada en el punto predictor
C_heun(i) = C_heun(i-1) + (h / 2) * (dC_initial + dC_pred); // Paso corrector
end
// MOSTRAR RESULTADOS EN CONSOLA
disp("Día | Heun (µg/m³)");
for i = 1:N
disp([t(i), C_heun(i)]);
end
// GRAFICAR RESULTADOS
clf;
plot(t, C_heun, 'b-x', "LineWidth", 2);
xlabel("Tiempo (días)");
ylabel("Concentración de CO2 (µg/m³)");
title("Evolución de la Concentración de CO2 (Método de Heun)");
legend("Heun", "Location", "best");
xgrid();
10
11
INTERPRETACIÓN DE RESULTADOS
Los resultados obtenidos con los métodos de Euler y Heun muestran que la concentración
de CO₂ en la atmósfera disminuye con el tiempo debido a la dispersión y absorción,
estabilizándose alrededor de 36,36 µg/m³. El método de Euler, al basarse solo en la
pendiente inicial, presenta una mayor acumulación de errores y subestima ligeramente la
concentración, mientras que el método de Heun, al corregir con un promedio de
pendientes, ofrece una estimación más precisa y estable. Ambos métodos confirman que la
concentración de CO₂ no crece indefinidamente, sino que tiende a un valor de equilibrio
determinado por la relación entre la emisión, dispersión y absorción del gas en la
atmósfera.
Desde una perspectiva ambiental, estos modelos pueden utilizarse para predecir la
evolución de contaminantes y evaluar el impacto de las estrategias de reducción de
emisiones. Si la tasa de emisión supera la capacidad de dispersión y absorción, la
concentración de CO₂ aumentará progresivamente, afectando la calidad del aire y
contribuyendo al cambio climático. Por otro lado, políticas que reduzcan las emisiones o
mejoren la absorción atmosférica pueden ayudar a mantener niveles más bajos de
contaminación. En conclusión, el método de Heun es más recomendable para estudios
ambientales debido a su mayor precisión, mientras que Euler puede ser útil para
aproximaciones rápidas en cálculos más simples.
BIBLIOGRAFÍA
Dióxido de carbono | Signos vitales – Climate Change: Vital Signs of the Planet. (s. f.).
Climate Change: Vital Signs Of The Planet.
[Link]
Método de Euler (s/f). Universidad Nacional de Mar de la Plata. [Link]. Recuperado el 8 de
febrero de 2025, de [Link]
%[Link]
