INTITUTO TECNOLOGICO DE MINATITLAN

METODS NUMERICOS

PROFESOR ASIGNADO: ALVARADO ALVARADO

ANDRES

UNIDAD 6: SOLUCION DE ECUCIONES

DIFERENCIALES

INTEGRANTES DE EQUIPO:
CRUZ SEBA ANA GUADALUPE
TORRES REYES YARA NAYROBI
Unidad 6 Solución de ecuaciones diferenciales (Valor inicial y valor en la frontera)

6.1 Fundamentos
6.2 Métodos de un paso
6.3 Métodos rígidos y de paso múltiples
6.4 Métodos multipaso
6.5 Métodos de tamaño de paso variable
6.6 Sistemas de ecuaciones diferenciales ordinarias
6.7 Solución de ecuaciones diferenciales ordinarias de orden n
6.8 Métodos generales para problemas con valores en la frontera (Lineal y no lineal)
6.9 Clasificación de ecuaciones diferenciales parciales
6.10 Aplicaciones
 6.1 Fundamentos

Una ecuación diferencial es una ecuación en la que intervienen derivadas de una o

más funciones desconocidas dependiendo del número de variables independientes
respecto de las que se deriva, las ecuaciones diferenciales se dividen en:
Ecuaciones diferenciales ordinarias: Aquellas que contienen derivadas respecto
a una sola variable independiente.
Ecuaciones en derivadas parciales: Aquellas que contienen derivadas respecto
a dos o más variables.
Una ecuación diferencial es una ecuación que incluye expresiones o términos que
involucran a una función matemáticas con incógnita y sus derivadas. Algunos
ejemplos de ecuaciones diferenciales son:
𝒀’ = 𝟐𝒙𝒚 + 𝟏
Es una ecuación diferencial ordinaria, donde representa una función no especificada
𝒅𝒚
de la variable independiente “𝑥”, es decir, 𝒀 = 𝒇(𝒙), 𝒚´ = 𝒅𝒙 , donde es la derivada
de “𝑦” con respecto a “𝑥”.
𝒂𝒖 𝒂𝒖
La expresión 𝒂𝒙 + 𝒂𝒚 = 𝟎 es una ecuación en derivadas parciales

A la variable dependiente también se le llama función incógnita (desconocida). La

resolución de ecuaciones diferenciales es un tipo de problema matemático que
consiste en buscar una función que cumpla una determinada ecuación diferencial.
Se puede llevar a cabo mediante un método específico para la ecuación diferencial
en cuestión o mediante una transformada (Por ejemplo, la transformada de
Laplace).
Orden de la ecuación:
El orden de la derivada más alta en una ecuación diferencial se denomina orden de
la ecuación.
Grado de la ecuación:
Es la potencia de la derivada de mayor orden que aparece en la ecuación, siempre
y cuando la ecuación este en forma polinómica, de no ser así se considera que no
tiene grado.
Ecuación diferencial lineal:
Se dice que una ecuación es lineal si tiene la forma: 𝒂𝒏(𝒙)𝒚(𝒏−𝟏) +, … , +𝒂𝟏(𝒙)𝒚´ +
𝒂𝟎(𝒙)𝒚 = 𝒈(𝒙) , es decir:
- Ni la función, ni sus derivadas están elevadas a ninguna potencia distinta de uno
o cero.
- En cada coeficiente que aparece multiplicándolas solo interviene la variable
independiente.
- Una combinación lineal de sus soluciones es también solución de la ecuación.
Ejemplos
 𝒚´ = 𝒚
Es una ecuación diferencial ordinaria lineal de primer orden, tiene como soluciones
𝒚 = 𝒇(𝒙) = 𝒌. 𝒆𝒙 , Con “k“, Siendo un número real cualquiera.
 −𝒀𝒏 + 𝒚 = 𝟎
Es una ecuación diferencial ordinaria lineal de segundo orden, tiene como
soluciones:
𝒀 = 𝒇(𝒙) = 𝒂𝒄𝒐𝒔(𝒙) + 𝒃𝒔𝒊𝒏(𝒙), con a y b reales.
 𝒀𝒏 − 𝒚 = 𝟎
Es una ecuación diferencial ordinaria lineal de segundo orden, tiene como
soluciones:
𝟏
𝒂. 𝒆𝒙 + 𝒃. 𝒆𝒙 , con a y b reales.

6.2 Métodos de un paso

Método de Euler
El método de Euler, Es un procedimiento de integración numérica para resolver
ecuaciones diferenciales ordinarias a partir de un valor inicial dado.
El método de Euler es el más simple de los métodos numéricos para resolver un
problema del siguiente tipo:

𝒅𝒚
= 𝒇(𝒙, 𝒚)
𝒅𝒙

PVI 𝒚(𝒙𝟎) = 𝒚𝟎
𝒚(𝒙𝒊) = ?
Consiste en multiplicar los intervalos que va de “𝒙𝟎” a “𝒙𝒇” en “n” sub-intervalos de
ancho h; ósea:
𝒙𝒇−𝒙𝟎
𝒉= , De manera que se obtiene un conjunto discreto de 𝒏 + 𝟏 puntos:
𝒏
𝒙𝟎, 𝒙𝟏, 𝒙𝟐, … 𝒙𝒏 del intervalo de interés [𝒙𝟎, 𝒙𝒇].
Para cualquiera de estos puntos se cumple que:
𝑿𝒊 = 𝒙𝒐 + 𝒊𝒉, 𝟎 ≤ + ≤ 𝒏
La condición inicial 𝒚(𝒙𝟎) = 𝒚𝟎, representa el punto 𝑷𝟎 = (𝒙𝟎, 𝒚𝟎), por donde pasa
la curva.
La solución de la ecuación del planteamiento inicial, se denotara como 𝒇(𝒙) = 𝒚.
Teniendo el punto 𝒑𝟎, se puede evaluar la primera derivada de 𝒇(𝒙) en ese punto;
por lo tanto:
𝒅𝒚
𝒇´(𝒙) = 𝒑𝒐 = 𝒇(𝒙𝟎, 𝒚𝟎)
𝒅𝒙

Ejemplo
𝒀´ = 𝟐𝒙 − 𝟑𝒚 + 𝟏, 𝒚(𝟏) = 𝟓, 𝒚(𝟏. 𝟐)
Primer caso 𝒉 = 𝟎, 𝟏
𝒅𝒚
Primero. Escribimos la Ec. Dif. En la forma = 𝒇(𝒙, 𝒚), para extraer su segundo
𝒅𝒙
𝒅𝒚
miembro 𝒅𝒙 = 𝟐𝒙 + 𝟑𝒚 + 𝟏

Definimos
𝒙𝟎, 𝒚𝟎 𝒚 𝒉 De acuerdo a los datos del problema
𝑿𝟎 = 𝟏
𝒀𝟎 = 𝟓
Planteamos la ecuación de Euler utilizando los datos iniciales:
𝒀𝟎 + 𝟏 = 𝒚𝟎 + 𝒉𝒇(𝒙𝟎, 𝒚𝟎)
𝒀𝟏 = 𝒚𝟎 + 𝒉 ∗ (𝟐𝒙𝟎 − 𝟑𝒚𝟎 + 𝟏)
𝒀𝟏 = 𝟓 + (𝟎. 𝟏)(𝟐(𝟏) − 𝟑(𝟓) + 𝟏)
Desarrollamos hasta el valor buscado en x, en este caso:
𝑿 = 𝟏. 𝟐
𝒀𝟏 = 𝟓 + (𝟎. 𝟏)(𝟐 − 𝟏𝟓 + 𝟏)
= 𝟓 + (𝟎. 𝟏)(−𝟏𝟐)
= 𝟓 − 𝟏. 𝟐
𝒀𝟏 = 𝒚(𝟏. 𝟏) = 𝟑. 𝟖𝟎𝟎𝟎

𝒀𝟏 + 𝟏 = 𝒚𝟏 + 𝒉𝒇(𝒙𝟏, 𝒚𝟏)
𝒀𝟐 = 𝒚𝟏 + 𝒉 ∗ (𝟐𝒙𝟏 − 𝟑𝒚𝟏 + 𝟏)
𝒀𝟐 = 𝟑. 𝟖 + (𝟎. 𝟏)(𝟐(𝟏. 𝟏) − 𝟑(𝟑. 𝟖) + 𝟏)
= 𝟑 . 𝟗 + (𝒐. 𝟏)(𝟐. 𝟐 − 𝟏𝟏. 𝟒 + 𝟏)
= 𝟑. 𝟖 + (𝟎. 𝟏)(−𝟖. 𝟐)
= 𝟑. 𝟖 − 𝟎. 𝟖𝟐
𝒚𝟐 = 𝒚(𝟏. 𝟐) = 𝟐. 𝟗𝟖𝟎𝟎

Método de Euler mejorado

Este método se basa en la misma idea del método anterior, pero hace un
refinamiento en la aproximación, tomando un promedio entre ciertas pendientes.
La fórmula es la siguiente:
𝒇(𝒙𝒏, 𝒚𝒏) + 𝒇(𝒙𝒏 + 𝟏, 𝒚𝒏 + 𝟏)
𝒀𝒏 + 𝟏 = 𝒚𝒏 + 𝒉 [ ]
𝟐

Ejemplo
Aplicar el método de Euler mejorado para aproximar 𝒚(𝟏. 𝟑) si tenemos:
𝒀´ = 𝒙 − 𝒚 + 𝟓
𝒀(𝟏) = 𝟐
Datos en una primera iteración tenemos:
𝑯 = 𝟎. 𝟏
𝑭(𝒙, 𝒚) = 𝒙 − 𝒚 + 𝟓 𝒙𝟏 = 𝒙𝟎 + 𝒉 = 𝟏. 𝟏
𝑿𝟎 = 𝟏 𝒚𝟏 = 𝒚𝟎 + 𝒉. 𝒇(𝒙𝟎, 𝒚𝟎) = 𝟐. 𝟒
𝟒+(𝟏.𝟏−𝟐.𝟒+𝟓)
𝒀𝟎 = 𝟐 𝒚𝟏 = 𝟐 + 𝟎. 𝟏 ( ) = 𝟐. 𝟑𝟖𝟓
𝟐

Resumimos los resultados en la siguiente tabla

𝒏 𝒙𝒏 𝒚𝒏
0 1 2
1 1.1 2.385
2 1.2 2.742925
3 1.3 3.07635

Concluimos entonces que la aproximación buscada es:

𝒀(𝟏. 𝟑) = 𝟑. 𝟎𝟕𝟖𝟔𝟑𝟓

Método de Runge-kutta
Es un método genérico de resolución genérica de ecuaciones diferenciales. Son un
conjunto de métodos iterativos (explícitos e implícitos) para la aproximación de
soluciones de ecuaciones diferenciales ordinarias, concretamente del problema de
valor inicial.
Sea 𝒚´(𝒕) = 𝒇(𝒕, 𝒚(𝒕)).
Una ecuación diferencial ordinaria, con 𝐟: 𝛀 𝐂𝐑𝐱𝐑´´ → 𝐑´´, Donde omega es un
conjunto abierto, junto con la condición de que el valor inicial de “𝒇” sea (𝐭𝟎, 𝐲𝟎)𝛆𝛀
Entonces el método Rk (de orden s) tiene la siguiente expresión, en su forma más
general
𝑠

𝑦𝑛 + 1 = 𝑦𝑛 + ℎ ∑ 𝑏𝑖𝑘𝑖
𝑖=1

Donde “𝒉“es el paso por iteración, o lo que es lo mismo el incremento ∆𝒕𝒏 entre los
sucesivos puntos 𝒕𝒏 y 𝒕𝒏 + 𝟏. Los coeficientes ”𝒌𝒊“ son términos de aproximación
intermedios, evaluados en t de manera local.
𝒔

𝒌𝒊 = 𝒇(𝒕𝒏 + 𝒉𝒄𝒊, 𝒚𝒏 + 𝒉 ∑ 𝒂𝒊𝒋 𝑲𝒊

𝒋=𝟏
Con 𝒂𝒊𝒋, 𝒃𝒊, 𝒂 coeficientes propios del esquema numérico elegido dependiente de la
regla cuadrática utilizada. Los esquemas runge-kutta pueden ser explícitos o
implícitos dependiendo de las constantes ”𝒂𝒊𝒋“del esquema.

Ejemplo
𝟎. 𝟓
𝒚(𝟎. 𝟓) = 𝒚(𝟎) + [(𝟒𝒆𝟎.𝟖(𝟎) − 𝟎. 𝟓(𝟐)) + 𝒇(𝟎. 𝟓, 𝟐 + 𝟎. 𝟓𝒇(𝟎. 𝟐))]
𝟐
= 𝟐 + 𝟎. 𝟓[𝟑 + 𝒇(𝟎. 𝟓, 𝟑. 𝟓)] = 𝟑. 𝟖𝟎𝟒𝟑𝟐𝟓
𝟎. 𝟓
𝒀(𝟏. 𝟎) = 𝒚(𝟎. 𝟓) + [𝒇(𝟎. 𝟓, 𝟑. 𝟖𝟎𝟒𝟑𝟐𝟓) + 𝒇(𝟎. 𝟓 + 𝟎. 𝟓, 𝟑. 𝟖𝟎𝟒𝟑𝟐𝟓
𝟐
+ 𝟎. 𝟓𝒇(𝟎. 𝟓, 𝟑. 𝟖𝟎𝟒𝟑𝟐𝟓))] = 𝟔. 𝟑𝟏𝟔𝟓𝟑𝟖
𝟎. 𝟓
𝒀(𝟏. 𝟓) = 𝒚(𝟏. 𝟎) + [𝒇(𝟏. 𝟎, 𝟔. 𝟑𝟏𝟔𝟓𝟑𝟖) + 𝒇(𝟏. 𝟎 + 𝟎. 𝟓, 𝟔. 𝟑𝟏𝟔𝟓𝟑𝟖
𝟐
+ 𝟎. 𝟓𝒇(𝟏. 𝟎, 𝟔. 𝟑𝟏𝟔𝟓𝟑𝟖)) = 𝟗. 𝟗𝟐𝟒𝟎𝟔𝟖
𝒀(𝟐. 𝟎) = 𝟏𝟓. 𝟏𝟗𝟔𝟐𝟗𝟖
𝒀(𝟐. 𝟓) = 𝟐𝟐. 𝟗𝟕𝟓𝟗𝟑𝟖
𝒀(𝟑. 𝟎) = 𝟑𝟒. 𝟓𝟏𝟒𝟗𝟐𝟎
𝒀(𝟑. 𝟓) = 𝟓𝟏. 𝟔𝟕𝟔𝟖𝟏𝟏
𝒀(𝟒. 𝟎) = 𝟕𝟕. 𝟐𝟑𝟖𝟓𝟐𝟒

6.3 Métodos rígidos y de pasos múltiples

Método rígido
Las ecuaciones rígidas son aquellas cuyas soluciones contienen escalas diferentes
para la variable independiente cuando la escala más grande es la de interés, pero
la escala más pequeña dicta el tamaño de paso de un método con base en la
estabilidad.
EDO (Ecuaciones diferenciales ordinarias), tanto las EDO como los sistemas de
EDO pueden ser rígidas.
 Ejemplo
Una EDO rígida es la que se muestra a continuación:
𝒅𝒚
𝟏. = −𝟏𝟎𝟎𝟎𝒚 + 𝟑𝟎𝟎𝟎 − 𝟐𝟎𝟎𝟎𝒆−𝟏
𝒅𝒕
Si se considera la condición inicial 𝒚(𝟎) = 𝟎, la solución analítica que se obtiene
esta dada por:
𝟐. 𝒚 = 𝟑 = 𝟎. 𝟗𝟗𝟖𝒆−𝟏𝟎𝟎𝟎𝒕 − 𝟐. 𝟎𝟎𝟐𝒆−𝒕
La solución, al principio se encuentra denominada por el termino exponencial rápido
𝒆−𝟏𝟎𝟎𝟎𝒕 , después de un periodo muy corto (𝒕 < 𝟎. 𝟎𝟎𝟓), esta parte transitoria termina
y la solución se rige por el exponencial lento 𝒆−𝒕 .
Analizando la parte homogénea de la ecuación (1), se puede determinar el tamaño
de paso necesario para que la ecuación sea estable, sea entonces la ecuación
homogénea.
𝒅𝒚
𝟑. = −𝒂𝒚
𝒅𝒕
Con la condición 𝒚(𝟎) = 𝒚𝟎, se obtiene la solución que se muestra en (4), que
asintóticamente se aproxima a cero, comenzando en el valor y0.
𝟒. 𝒚 = 𝒚𝟎𝒆−𝒂𝒕
Si se aplica el método de Euler a la ecuación (3), se obtiene la formula:
𝟓. 𝒚𝒊 + 𝟏 = 𝒚𝒊 − 𝒂𝒚𝒊𝒉
Que trabajada algebraicamente queda:
𝟔. 𝒚𝒊 + 𝟏 = 𝒚𝒊(𝟏 − 𝒂𝒉)
Por lo que podemos decir que 𝒚𝒊 = 𝒚𝟎(𝟏 − 𝒂𝒉)𝒊 , para que esta solución numérica
𝟐
sea acotada deberá se [𝟏 − 𝒂𝒉] < 𝟏. De aquí surge que si 𝒉 > , entonces [𝒚𝒊]
𝒂
crece indefinidamente cuando i tiende a infinito.

Métodos de pasos múltiples

Se considera el problema de valores iniciales (P,V,I) 𝟖 <: 𝒚𝟎(𝒙) =
𝒇(𝒙, 𝒚(𝒙)); 𝒙𝟐 [𝒂; 𝒃] , 𝒚(𝒂) = 𝒚𝟎 dado el que supondremos tiene solución única y
[𝒂; 𝒃]
 Ejemplo
Dada una partición del intervalo [𝒂; 𝒃], 𝒂 = 𝒙𝟎 < 𝒙𝒏 = 𝒃; los métodos que hemos
visto hasta aquí solo usan la información del valor ”𝒚𝒊“ de la ecuación.

6.4 Métodos Multipaso

En este tipo de métodos, para calcular un valor en la malla de puntos se usan varios
de los valores previamente calculados y no solo el valor anterior.
Donde hemos usado integración numérica para la integral de 𝒇(𝒕, 𝒚). al despejar
𝒚𝒊 + 𝟏 se tiene
𝒎 𝒎

𝒚𝒊 = ∑ 𝒂𝒓𝒚𝒊 + 𝟏 − 𝒓 + 𝒉 ∑ 𝒃𝒓 + (𝒕𝒊 + 𝟏 − 𝒓, 𝒚𝒊 + 𝟏 − 𝒓)
𝒓=𝟏 𝒓=𝟎

Si 𝒃𝟎 = 𝟎, lo que queremos calcular (𝒚𝒊 + 𝟏) esta en ambos ladosm entonces el

método se llama implicitom en otro caso, explicito. Ahora no se toman puntos
intermedios 𝜃𝑟, sino puntos anteriores en los cuales tenemos la información que
necesitamos.
Cuando empezamos, en un PVI de primer orden, solo tenemos un valor, así que
para preparar el principio de tabla hasta los m valores que necesita el método
multipaso, se usan métodos unipaso.
Veamos el caso en que 𝒃𝒐 = 𝟎:
𝒎

𝒚𝒊 + 𝟏 = 𝒗𝟏, 𝒊 + 𝒉𝒃𝟎𝒇(𝒕𝒊 + 𝟏, 𝒚𝒊 + 𝟏) + 𝒉 ∑ 𝒃𝒓 + (𝒕𝒊 + 𝟏 − 𝒓, 𝒚𝒊 + 𝟏 − 𝒓)

𝒓=𝟏

𝒚𝒊 + 𝟏 = 𝒂 + 𝒉𝒃𝟎𝒇(𝒕𝒊 + 𝟏, 𝒚𝒊 + 𝟏)
Los métodos multipaso se pueden deducir por integración de la ecuación diferencial
y usando interpolación polinomial o utilizando el método de los coeficientes
interminados. Imponiendo que la ecuación en diferencia sea exacta para todo
problema de solución un polinomio de bajo grado.
Ejemplo
Resolver el ejemplo utilizando los métodos de Adams del mismo orden en la forma
predictor-corrector, suponiendo que calculáramos por un método unipaso una
estimación de 𝑢𝑖
𝝅𝒕
𝒕𝒚𝒏 + 𝒚𝒚´ + 𝒔𝒆𝒏 =𝟎
𝟔
𝒚(𝟏) = 𝟏, 𝒚´(𝟏) = 𝟐
Solución
Teniendo el valor aproximado de vi, del ejemplo mencionado y usando los métodos
explicito como predictor y el implícito como corrector tenemos que usar:
𝒉
𝒚𝒊 + 𝟏 = 𝒚𝒊 + (𝟑𝒇𝒊 − 𝟏)
𝟐
Por ejemplo, para un método de segundo orden, predecimos el valor de 𝑢2:

𝒚(𝟐) 𝒉
( 𝒚´(𝟐)= 𝒖𝒑𝟐 = 𝒖𝟏 + (𝟑𝒇(𝟏. 𝟓𝒖𝟏) − 𝒇(𝟏, . 𝒋𝟎))
𝟐
Vemos que el valor de la función crece, lo cual es un reflejo de que 𝒚´(𝟏. 𝟓) > 𝟎, tal
y como se ve eje el ejemplo anterior.
Ahora corregimos el valor

𝒚(𝟐) 𝟎. 𝟓
( 𝒚´(𝟐)= 𝒖𝟐−𝒄 𝒖𝟏 + (𝒇𝒄𝟐, 𝟐−𝒑 ) + 𝒇(𝟏. 𝟓, 𝒖𝟏))
𝟐
= (𝟐.𝟎𝟑𝟐𝟖𝟏
𝟎.𝟑𝟔𝟑𝟑 )

Si lo hacemos con el siguiente punto

𝒚(𝟐.𝟓) 𝟎. 𝟓
( 𝒚´(𝟐.𝟓)) = 𝒖𝟑−𝒑 = 𝒖𝟐−𝒄 + (𝟑𝒇((𝟐, 𝒖−𝒄 𝟐) − 𝒇(𝟏. 𝟓, 𝒖𝟏))
𝟐
Con lo que aplicando repetidamente, tenemos la solución.

6.5 Métodos de tamaño de paso variable

Los métodos desarrollados anteriormente tiene un gran tamaño de paso fijo (puntos
igualmente especificados) sin embargo, esto no permite tener un control sobre el
error de truncamiento local en cada paso, ya que puede darse el caso de que la
función tenga cambios buscos en el intervalo de integración.
Se han desarrollado técnicas numéricas para la estimación local del error que
permitan controlar el tamaño de paso óptimo para controlar el error global.

Se estudiara el error local con el objetivo de controlar indirectamente el error global

cometido en la resolución de la EDO.
Considérese un método de un paso de orden de consistencia P; El error local de
truncamiento es:
𝜺𝒏 + 𝟏 = 𝒚(𝒙𝒏 + 𝟏) − 𝒚𝒏 + 𝟏 = 𝒄(𝒙𝒏)𝒉𝒏𝒑+𝟏 + 𝟎(𝒉𝒏𝒑+𝟐 )
Siendo Y la solución exacta que paso por (𝑿𝒏, 𝒀𝒏). Remitiendose el análisis de la
evolución del error global, si el error local por unidad de longitud satisface
𝜺𝒏+𝟏
| |<𝟎
𝒉𝒏

Y las condiciones de unidad de la solución se cumplen, entonces el error global de

truncamiento satisface
𝒆𝟐(𝒙𝒏−𝒂) − 𝟏
|𝜺𝒏| < 𝒅 | | , ∀𝒙𝒏𝜺[𝒂, 𝒃]
𝟏

La estrategia a adoptar para la selección de cada paso deberá garantizar que

satisfaga la cota (2), la primera dificultad que aparece detrás de esta estrategia es
que el valor de 𝒉𝒏 no se puede calcular directamente de lo anterior porque
generalmente la función c(x) de (1) no se conoce, aunque se asumirá que la misma
varia suavemente.
Si 𝛆∗ 𝐧 + 𝟏 es el error local asociado con el paso H, entonces
𝛆∗ 𝐧 + 𝟏 = 𝐜(𝐱𝐧) + 𝑯𝒑+𝟏
Y el error local producido por un paso es asumiendo el mismo valor de
𝛆∗ 𝐧 + 𝟏 = 𝐜(𝐱𝐧) + (𝛉𝑯)𝒑+𝟏

La desigualdad (2) se puede satisfacer en caso que

|𝒄(𝒙𝒏)(𝜽𝑯)𝒑 | < 𝒐
Usando la estimación (3) para eliminar 𝒄(𝒙𝒏), se tiene que la elección
𝟎𝑯 𝟏
𝜽<| |
𝜺𝒏 + 𝟏 𝒑
Cumple con la cota del error local propuesta en (2)
En términos de costo de las operaciones, es conveniente tomar el mayor paso
posible, pero como rechazar un paso puede resultar inadecuado se toma:
𝜽𝑯
𝜽 = 𝟎. 𝟖 | |
𝝐𝒏 + 𝟏
 6.6 Sistemas de ecuaciones diferenciales ordinarias
En un sistema de ecuaciones diferenciales ordinarias de cualquier orden, pueden
ser reducidos a un sistema equivalente de primer orden si se introducen nuevas
variables y ecuaciones. Por esa razón en este artículo solo se consideran sistemas
de ecuaciones de primer orden. Un sistema de ecuaciones diferenciales ordinarias
de primer orden escrito en forma explícita es un sistema de la forma:
𝒅𝒙𝒊
= 𝒇𝟏(𝒙𝟏, 𝒙𝟐, … , 𝒙𝒏; 𝒕)
𝒅𝒕
𝒅𝒙𝟐
= 𝒇𝟐(𝒙𝟏, 𝒙𝟐, … , 𝒙𝒏; 𝒕)
𝒅𝒕
…
𝒅𝒙𝒏
= 𝒇𝒏(𝒙𝟏, 𝒙𝟐, … , 𝒙𝒏; 𝒕)
𝒅𝒕

Reducción a un sistema de primer orden

Dada un sistema de ecuaciones diferenciales de orden n con m ecuaciones:
𝒅𝒙𝒋 𝒅𝒏 𝒙𝒋
𝒇𝒊 (𝒙𝒋, , … , 𝒏 ; 𝒕) = 𝟎 𝒄𝒐𝒏 𝒊, 𝒋
𝒅𝒕 𝒅𝒕
Existe un sistema equivalente de primer orden con a la suma (𝒏 + 𝟏)𝒙𝒎 ecuaciones
para ver esto consideramos un sistema en que interviene m funciones incógnitas
“𝑥𝑖" y sus “𝒏“derivadas e introduzcamos un nuevo conjunto de variables ”𝒀𝒊”,”𝒌”
definidos de la siguiente manera:
𝒅𝒌 𝒙𝟏(𝒕)
𝒚𝒊, 𝒌(𝒕) =
𝒅𝒕𝒌
El sistema de primer orden equivalente en las variables “𝒚𝒊”, “𝒌” resulta ser:
𝒅𝒚𝒊,𝒌
𝒚𝒊, 𝒌 + 𝟏 = 𝒌𝜺{𝟎, 𝟐, . . }
𝒅𝒕

𝒇𝒊(𝒚𝟏, 𝟎, 𝒚𝒊, 𝟏, … , 𝒚𝒋, 𝒏; 𝒕) = 𝟎 𝟏, 𝒋𝜺{𝟏, 𝟐, … }

Como ejemplo de reducción de un sistema de ecuaciones diferenciales podemos

considerar las ecuaciones de movimiento de la mecánica newtoniana de una
partícula que es un sistema de segundo orden de tres ecuaciones:
 𝒅𝟐 𝒙 𝒅𝒙 𝒅𝒚 𝒅𝒛
𝒎 = 𝒇𝒙 (𝒙, 𝒚, 𝒛; , ; ; 𝒕)
𝒅𝒕𝟐 𝒅𝒕 𝒅𝒕 𝒅𝒕
𝒅𝟐 𝒚 𝒅𝒙 𝒅𝒚 𝒅𝒛
𝒎 𝟐
= 𝒇𝒚 (𝒙, 𝒚, 𝒛; , ; ; 𝒕)
𝒅𝒕 𝒅𝒕 𝒅𝒕 𝒅𝒕
𝒅𝟐 𝒛 𝒅𝒙 𝒅𝒚 𝒅𝒛
𝒎 𝟐 = 𝒇𝒛 (𝒙, 𝒚, 𝒛; , ; ; 𝒕)
𝒅𝒕 𝒅𝒕 𝒅𝒕 𝒅𝒕
Si se introducen tres funciones incógnitas nuevas que representan la velocidad, el
sistema anterior se puede reducir a un sistema de primer orden y seis ecuaciones.
𝒅𝒙 𝒅𝒗𝒙
𝒎 = 𝒗𝒙(𝒙), 𝒎; = 𝒇𝒙(𝒙, 𝒚, 𝒛, 𝒗𝒙)
𝒅𝒕 𝒅𝒕
𝒅𝒚 𝒅𝒗𝒚
𝒎 = 𝒗𝒚(𝒕), 𝒎; = 𝒇𝒚(𝒙, 𝒚, 𝒛, 𝒗𝒙)
𝒅𝒕 𝒅𝒕
𝒅𝒛 𝒅𝒗𝒛
𝒎 = 𝒗𝒙(𝒛), 𝒎; = 𝒇𝒛(𝒙, 𝒚, 𝒛, 𝒗𝒙)
𝒅𝒕 𝒅𝒕

Ejemplo
La población 𝒑(𝒕) de un suburbio de una gran ciudad en un instante cualquiera se
rige por
𝒅𝒑
= 𝑷(𝟏𝟎−𝟏 − 𝟏𝟎−𝟕 𝒑)
𝒅𝒕
𝒑(𝟎) = 𝟓𝟎𝟎𝟎

En donde “𝒕” se mide en meses. ¿Cuál es el valor límite de la población?, ¿en qué
momento será la población igual a la mitad de su valor límite?
Solución: calculamos el tamaño de la población, 𝒑(𝒕), resolviendo el problema de
valores iniciales. La ecuación diferencial tiene sus variables separadas:
𝒑´
=𝟏
𝒑(𝟏𝟎 − 𝟏 − 𝟏𝟎−𝟕 𝒑)

𝒅𝒑
Donde hemos denotado 𝒑𝟏 = , integrando los dos miembros de esta identidad
𝒅𝒕
entre “𝟎” y ”𝒕“obtenemos.
𝒑(𝒕)
𝒅𝑸
𝟏𝟎𝟕 ∫ 𝟔
=𝒕
𝟓𝟎𝟎𝟎 𝑸(𝟏𝟎 − 𝑸)
Donde hemos efectuado el cambio de variable 𝑸 = 𝒑(𝒕)
Teniendo en cuenta que
𝟏 −𝟔
𝟏 𝟏
= 𝟏𝟎 ( + )
𝑸(𝟏𝟎𝟔 − 𝑸) 𝑸 𝟏𝟎𝟔 − 𝑸
Concluimos tras una serie de cálculos simples que la única solución de nuestro
problema es:
𝟏
𝟏𝟎𝟔 𝒆𝟏𝟎
𝒑(𝒕) =
𝟏𝟗𝟗 + 𝒆𝟏/𝟏𝟎
Por lo tanto el valor límite de la población es:
𝐥𝐢𝐦 𝐩𝟖𝐭) = 𝟏𝟎𝟔
𝐭→∞

Como se desprende de una simple aplicación de la regla, tenemos que encontrar el

𝟏𝟎𝟔
valor “𝒕𝟎” para que 𝒑(𝒕𝟎)´ = , resolvemos la ecuación
𝟐

𝟏𝟎
𝟏𝟎−𝟔 𝒆 𝟏𝟎 𝟏𝟎𝟔 𝟏𝟎
= ↔ 𝒆𝟏𝟎 = 𝟏𝟗𝟗
𝟏𝟎 𝟐
𝟏𝟗𝟗 + 𝒆 𝟏𝟎

Concluimos que
𝒕𝟎 = 𝟏𝟎𝒍𝒐𝒈(𝟏𝟗𝟗)𝒎𝒆𝒔𝒆𝒔 = 𝟒. 𝟒𝟏 𝒂ñ𝒐𝒔

6.7 Solución de ecuaciones diferenciales ordinarias de orden n

Una ecuación ordinaria de orden n es una ecuación que liga la variable
independiente x, una función incógnita 𝒚 = 𝒚(𝒙) y sus derivadas sucesivas
𝒚´, 𝒚´´, … , 𝒚𝒏, es decir, es una expresión de la forma

𝒇(𝒙, 𝒚, 𝒚´, 𝒚𝒏 , . . , 𝒚(𝒏) ) = 𝟎 Forma implicita

O bien, si se puede despejar la derivada de mayor orden

𝒇(𝒙, 𝒚, 𝒚´, 𝒚𝒏 , . . , 𝒚(𝒏−𝟏) ) Forma explicita

Solución de una ecuación diferencial ordinaria (EDO) n
Dada la ecuación diferencial de orden m:

𝒚(𝒏) = 𝒇(𝒙, 𝒚, 𝒚´, … , 𝒚(𝒏−𝟏) 𝒇: 𝑫𝑪𝟏𝑹𝒏+𝟏 → 𝟏𝑹

Se dice que 𝒛 = 𝒛(𝒙) 𝒛: 𝟏𝒄 𝒍𝑹 – 𝒍𝑹 es la solución de la ecuación diferencial si
satisface:
1. z es n veces derivable en 1

2. (𝒙, 𝒛(𝒙), 𝒛´(𝒙), … . , 𝒛(𝒏+𝟏) (𝒙))𝝐𝑫∀𝒙𝜺𝑰

3. 𝒛𝒏 (𝒙) = 𝒇´(𝒙, 𝒛(𝒙), 𝒛´(𝒙), … , 𝒛(𝒏−𝟏) , (𝒙))

Es decir, solución de una EDO, Es toda función que sustituida justamente con sus
derivadas en la ecuación conduciendo a una identidad.

Tipos de soluciones de EDO

Las soluciones de una EDO pueden ser de tres tipos
1. Solución general Solución de la ecuación diferencial en la que aparecen tantas
constantes arbitrarias como orden de la ecuación, en nuestro caso al ser de orden
n la solución general será una familia de curvas.
2. solución particular Es una solución que se obtiene al fijar los valores de las
constantes arbitrarias de la solución general.
3. solución singular Es una solución que no esta incluida en la solución general;
es decir, no se puede obtener a partir de ella asignando valores convenientes alas
constantes arbitrarias.
Formas en las que pueden aparecer las soluciones de una EDO
La solución de una ecuación diferencial puede venir dada de 3 formas distintas:
1. forma explicita: si la incógnita y viene despejada en función de la variable
independiente x.
2. forma implícita: si la solución viene expresada por una ecuación que liga la
incógnita y por la variable x.
3. forma paramétrica: si la función viene dada en función de un parámetro.
 Ejemplo
problema de cauchy
Un problema de cauchy o problema de valores iniciales, asociados a una EDO, es
un problema (𝒑) de la forma:

𝒚(𝒏) = 𝒇(𝒙, 𝒚, 𝒚′ , … , 𝒚(𝒏−𝟏) )

𝒚𝟖𝒙𝟎) = 𝒚𝟎
𝒚´(𝒙𝟎) = 𝒚´𝟎

𝒚(𝒏−𝟏) (𝒙𝟎) = 𝒚𝟎𝒏−𝟏

Teorema de existencia y unidad de soluciones de cauchy
Sea (𝒑) un problema de cauchy, o de valores iniciales asociada a una EDO es un
problema (p) de la forma que:
1. existencia: si f es continua en un entorno de punto (𝒙𝟎, 𝒚𝟎, 𝒚´𝟎, . . , 𝒚𝟎𝒏−𝟏 )
Entonces (𝒑) posee solución.
𝝏𝒇 𝝏𝒇 𝝏𝒇
2. unidad: si además las derivadas parciales 𝝏𝒚 , 𝝏𝒚 ,.., 𝝏𝒚(𝒏−𝟏) .existen y son continuas
en un entorno del punto (𝒙𝟎, 𝒚𝟎, 𝒚´𝟎, … , 𝒚𝟎𝒏−𝟏 ) entonces (𝒑) posee solución única.
Es interesante hacer un estudio teórico previo de un problema antes de intentar
calcular la solución general de la ecuación diferencial (𝒑):

𝒚(𝒙 + 𝟏) = 𝒚 + 𝒙𝒚(𝒙𝟏)
𝒚(𝟎) = 𝟎
𝒚´(𝟎) = 𝟎
…
𝒚(𝒙𝟏)(𝟎) = 𝟎

6.8 Métodos generales para problemas con valores en la frontera

(Lineal y no lineal)
En el caso en que las condiciones para la variable dependiente estén definidas para
distintos valores en el rango de la variable independiente (normalmente en los
extremos), tenemos una ecuación diferencial con valores en la frontera.
Método de disparo
Este método se basa en la conversión de un problema en la frontera a su
equivalente de un problema de valores iniciales, implementándose un esquema de
prueba y error para alcanzar una solución adecuada.
Solución de una ecuación diferencial lineal
Una ecuación diferencial es línea; si en ella no aparecen potencias de la variable
dependiente y sus derivadas, ni productos de la variable por sus derivadas o
productos entre derivadas.
Ejemplo
𝟖𝒅𝟐 𝒚 𝒅𝒚
Resolver la ecuación diferencial + 𝟏𝟔 𝒅𝒙 − 𝟒𝒚 = 𝟐𝟎 con la condición de
𝒅𝒙𝟐
frontera 𝒚(𝟎)𝟓𝒚 𝒚 𝒚(𝟐𝟎) = 𝟐
Transformación
𝒅 𝒅𝒚𝒊 𝒅𝒚
Si 𝒚 = 𝒚𝟏 entonces 𝟖 𝒅𝒙 ( 𝒅𝒙 ) + 𝟏𝟔 𝒅𝒙 − 𝟒𝒚𝒊 = 𝟐𝟎

𝒅𝒚𝒊 𝒅𝟐 𝒚𝟐
Si entonces 𝟖 𝒅𝒙 + 𝟏𝟔 𝒅𝒙 − 𝟒𝒚𝒊 = 𝟐𝟎
𝒅𝒙

Solución de una ecuación diferencial no lineal

Una ecuación diferencial es no lineal, si en ella aparecen potencias de la variable
dependiente y sus derivadas, productos de la variable dependiente por sus
derivadas o productos entre derivadas.
Ejemplo
𝒅𝟐 𝒚 𝟏 𝒅𝒚
Resolver la ecuación diferencial = (𝟖)(𝟑𝟐 + 𝟐𝒙𝟑 − 𝒚𝟏 𝒅𝒙
𝒅𝒙𝟐

Transformación
𝒅 𝒅𝒚𝒊 𝟏 𝒅𝒚
Si 𝒚 = 𝒚𝟏 entonces 𝒅𝒙 ( 𝒅𝒙 ) = (𝟖)(𝟑𝟐 + 𝟐𝒙𝟑 − 𝒚𝟏 𝒅𝒙

𝒅𝒚𝒊 𝒅𝒚𝟐 𝟏
Si entonces = (𝟖)(𝟑𝟐 + 𝟐𝒙𝟑 − 𝒚𝟏𝒚𝟐)
𝒅𝒙 𝒅𝒙
Interpolación cuadrática
𝒇(𝒙𝟏) − 𝒇(𝒙𝟎) 𝒇(𝒙𝟐) − 𝒇(𝒙𝟏) − 𝒇(𝒙𝟏) − 𝒇(𝒙𝟎)
𝒇(𝒙) = 𝒇(𝒙𝟎) + (𝒙 − 𝒙𝟎) +
𝒙𝟏𝒙𝟎 (𝒙𝟐 − 𝒙𝟏) (𝒙𝟏 − 𝒙𝟎)
Método de diferencias finitas
La solución numérica más común para resolver ecuaciones diferenciales en la
frontera, se basa en el uso de ecuaciones diferenciales finitas para evaluar las
derivadas, ya que sustituyendo estas por su equivalente en la ecuación diferencial,
esta se transforma en una ecuación algebraica en diferencias, la técnica incluye la
construcción de una retícula donde se ubican los puntos que representan el
fenómeno a estudiar.
Solución de una ecuación diferencial lineal
Ejemplo
𝟖𝒅𝟐 𝒚 𝒅𝒚
Resolver la ecuación diferencial + 𝟏𝟔 𝒅𝒙 − 𝟒𝒚 = 𝟐𝟎
𝒅𝒙𝟐

Construcción de la retícula
𝒚(𝟎) = 𝟓 𝒚𝟏(𝟐) =? 𝒚𝟐(𝟒) =? 𝒚𝟑(𝟔) =? 𝒚𝟒(𝟖) =?
𝒚𝟓(𝟏𝟎) =? 𝒚𝟔(𝟏𝟐) =? 𝒚𝟕(𝟏𝟒) =? 𝒚𝟖(𝟏𝟔) =? 𝒚𝟗(𝟏𝟖) =? 𝒚𝟏𝟎(𝟐𝟎) = 2
𝒚𝒊 + 𝟏 − 𝟐𝒚𝒊 + 𝒚𝒊 − 𝟏 𝒚𝒊 + 𝟏 − 𝒚𝒊 + 𝟏
𝟖( ) + 𝟏𝟔 ( ) − 𝟒𝒚𝒊 = 𝟐𝟎
𝒉𝟐 𝟐𝒉
𝒙𝒏 − 𝒙𝟎 𝟐𝟎 − 𝟎
𝒉= = =𝟐
𝒏 𝟏𝟎
𝟐𝒚𝒊 + 𝟏 − 𝟒𝒚𝒊 + 𝟐𝒚𝒊 − 𝟏 + 𝟒𝒚𝒊 + 𝟏 − 𝟓𝒚𝒊 − 𝟏 = 𝟐𝟎
Solución de una ecuación diferencial no lineal
Ejemplo
𝒅𝟐 𝒚 𝟏 𝒅𝒚
= (𝟖) (𝟑𝟐 + 𝟐𝒙𝟑 − 𝒚𝟏 𝒅𝒙) Con la condición de frontera 𝒚(𝟏) = 𝟏𝟕 𝒚 𝒚(𝟑) = 𝟑
𝒅𝒙𝟐

Construcción de la retícula
𝒚𝟎(𝟏) = 𝟏𝟕 𝒚𝟏(𝟏. 𝟐) =? 𝒚𝟐(𝟏. 𝟒) =? 𝒚𝟑(𝟏. 𝟔) =? 𝒚𝟒(𝟏. 𝟖) =?
𝒚𝟓(𝟐. 𝟎) =? 𝒚𝟔(𝟐. 𝟐) =? 𝒚𝟕(𝟐. 𝟒) =? 𝒚𝟖(𝟐. 𝟔) =? 𝒚𝟗(𝟐. 𝟖) =? 𝒚𝟏𝟎(𝟑) = 14.333
𝒚𝒊 + 𝟏 − 𝟐𝒚𝒊 + 𝒚𝒊 − 𝟏 𝟏 𝒚𝒊 + 𝟏 − 𝒚𝒊 − 𝟐
𝟐
= (𝟑𝟐 + 𝟐𝒙𝟑 − 𝒚𝒊 )
𝒉 𝟖 𝒉𝟐
𝒙𝒏 − 𝒙𝟎 𝟑 − 𝟏
𝒉= = = 𝟎. 𝟐
𝒏 𝟏𝟎
 𝟐𝟓𝒚𝒊 + 𝟏 − 𝟓𝟎𝒚𝒊 + 𝟐𝟓𝒚𝒊 − 𝟏 = 𝟒 + 𝟎. 𝟐𝟓𝒙𝟑 − 𝟐. 𝟓𝒚𝒊𝒚𝒊 + 𝟏 + 𝟐𝟓𝒚𝒊𝒚𝒊 − 𝟏 = 𝟎

6.9 Clasificación de ecuaciones diferenciales parciales

Si una variable y depende de más de una variable independiente, las derivadas de
u con respecto de una o más variables independientes, se llaman derivadas
parciales.
𝝏𝟐 𝒖 𝝏𝟐 𝒖 𝝏𝟐 𝒖
𝑨 𝟐 +𝑩 𝟐 +𝑪 𝟐+𝑫=𝟎
𝝏𝒙 𝝏𝒙 𝒚 𝝏𝒚
La ecuación diferencial parcial de segundo orden para dos variables independientes
cuya fórmula es:
𝝏𝒖 𝝏𝒖
Donde 𝑨,𝑩 y 𝑪 son funciones de 𝒙,𝒚 y 𝑫, es función de 𝒙, 𝒚, 𝝏𝒙 𝒚 𝝏𝒚

Dependiendo de la relación entre los coeficientes 𝑨, 𝑩, y 𝑪 se clasifican en:

Elíptica si 𝐵 2 − 4𝐴𝐶 < 0 Ecuación de Laplace
Parabólica si 𝐵 2 − 4𝐴𝐶 = 0 Ecuación de conducción de calor
Hiperbólica 𝐵 2 − 4𝐴𝐶 > 0 Ecuación de onda
Ecuación diferenciales parciales parabólicas
Estas ecuaciones aparecen en ingeniería cuando se estudian los fenómenos de
conducción de calor en estado transitorio, así como en el estudio de la difusión
molecular en el seno de un fluido, etc.
Ejemplo
Conducción del calor de una varilla aislada, cuyos extremos libres se encuentran a
distintas temperaturas.
Haciendo un balance de energía, se encuentra la ecuación que gobierna el flujo de
calor
𝝏𝟐 𝑻 𝝏𝑻
𝒂 𝝏𝟐 𝒙 = Donde a se denomina difusividad térmica
𝝏𝒙

Para aplicar el método de diferencias finitas se construye una retícula.

El método explicito requiere el valor en (𝒊, 𝒋) a partir de (𝒊 − 𝟏, 𝒋 − 𝒊), (𝒊𝒋 − 𝟏), (𝒊 +
𝟏, 𝒋 − 𝟏).
𝑻(𝒙, 𝟎) = 𝒕𝒊
Condiciones iniciales en la frontera izquierda
𝑻(𝒐, 𝒇) = 𝒕𝟎
𝑻(𝒙, 𝒕) = 𝒕𝟎
En el método implícito predice el valor en las (𝒊, 𝒋) a partir de (𝒊, 𝒋) mediante la
generación de un sistema de ecuaciones obtenidas de los nodos.

6.10 Aplicaciones
Ley de newton del enfriamiento
Aplicando la primera ley de la termodinámica o la esfera, y suponiendo que el calor
fluye tan rápidamente en la esfera que la temperatura es prácticamente la misma
en todos los puntos de la misma, el calor descifrado por la esfera se puede expresar
analíticamente por medio de la ecuación diferencial homogénea.
𝒅𝑻 𝒉𝑨
+ (𝑻 − 𝑻∞) = 𝟎
𝒅𝒕 𝑷𝑪𝑽
Donde
H= coeficiente de transferencia de calor
A= área de la esfera para transferencia de calor
r= densidad de la esfera
V= calor especifico de la esfera
C= calor especifico de la esfera
Fenómenos de transporte ll
Ejemplo
Una esfera de aluminio de 3cm de diámetro se calienta hasta una temperatura de
200°c, entonces en el instante t00, se coloca en el aire que se mantiene a una
temperatura de 30°c. Si el coeficiente promedio de transferencia de calor es de
20W/M°C, calcule el tiempo necesario para que la esfera alcance una temperatura
de 150°C
Suponga las siguientes propiedades del aluminio:
𝑲 = 𝟐𝟏𝟎/𝑾/𝑴°𝑪
𝑪 = 𝟎. 𝟗𝟖𝟓𝒋/𝒈°𝑪
𝒓 = 𝟐. 𝟕𝟐𝒈/𝒄𝒎𝟒
Solución
𝒉𝑨 𝒉(𝟒𝝅𝑹𝟐 𝟑𝒉 𝟑 + 𝟐𝟎
= 𝟑
= = = 𝟏𝟔. 𝟒𝟑𝒙𝟏𝟎−𝟒 𝟓−𝟏
𝑷𝑪𝑽 𝑷𝑪(𝟒/𝟑𝝅𝑹 𝑷𝑪𝑹 𝟐. 𝟕𝟓𝒙𝟏𝟎 ∗ 𝟎. 𝟖𝟗𝒙𝟏𝟎𝟑 ∗ 𝟏. 𝟓𝒙𝟏𝟎−𝟐
𝟑

𝒅𝑻
= −𝟏𝟔. 𝟒𝟑𝒙𝟏𝟎−𝟒 (±𝟑𝟎) 𝒕 = 𝟎. 𝟕 = 𝟐𝟎𝟎°𝑪
𝒅𝒕
EJEMPLO:
Para un sistema dado, el sistema de ecuaciones se puede expresar como:
𝒅𝒚
= 𝟐𝒙 − 𝟏. 𝟐𝒙𝒚 𝒆𝒏 𝒕 = 𝟎, 𝒙 = 𝟐
𝒅𝒕

𝒅𝒚
= −𝒚 + 𝟎. 𝟗𝒙𝒚 𝒆𝒏 𝒕 = 𝟎, 𝒚 = 𝟏
𝒅𝒕
Evaluación sumativa
EJEMPLO: Sistema de reactores tipo tanque con agitación
Considere el siguiente sistema de reactores tipo tanque donde:
Q=Flujo volumétrico en metros cúbicos por minuto
C= Concentracion en miligramos por metro cubico.
𝒎𝟑 𝒎𝒈 𝒎𝒈
Flujo masico = QC = 𝒎𝒊𝒏 = 𝒎𝒊𝒏
𝒎𝟑

Balance de materia (ley de conservación de la materia)

Acumulación = entradas – salidas
Acumulación = v dc/dt
V = volumen del reactor
𝑑𝑦
V 𝑑𝑥 = 𝐸𝑛𝑡𝑟𝑎𝑑𝑎𝑠 − 𝑆𝑎𝑙𝑖𝑑𝑎𝑠; 𝑟𝑒𝑎𝑐𝑜𝑚𝑜𝑑𝑎𝑛𝑑𝑜

𝑑𝑦 𝐸𝑛𝑡𝑟𝑎𝑑𝑎𝑠 − 𝑆𝑎𝑙𝑖𝑑𝑎𝑠
=
𝑑𝑥 𝑣
BIBLIOGRAFIAS
1. BURGEN L. ANALISIS NUMERICOS IBEROAMERICA
2. CHAPRAS S.C Y CANALE R. NUMENCAL METHODS POR ENGINENNG ,
MC GRAW-HILL.
3. CONSTATINIDES, A APPLIED NUMERICAL METHODS WITH PERSONAL
COMPUTERS. MC GRAW HILL.
4. CONTE. S.D Y DE BOOR C. ANALISIS NUMERICO ELEMENTAL. MC
GRAW –MILL.
5. LUTHE OLIVERA Y SCHUITZ. METODOS NUMERICOS. LIMUSA.
6. MATHEWS, J Y FINK. C.D METODOS NUMERICOS CON MATLAD.
PRENTICE HALL.
7. NAKAMURA ,S, METODOS NUMERICOS ALICADA CON SOFTWARE.
PRENTICE-HALL.
8. NIEVE.A. DOMINGUEZ F.METODOS NUMERICOS APLICADOS A
INGENIERIA C.E.C.S.A.
9. R.E, METODOS NUMERICO BASICOS. MC GRAW- HILL.
10. WALFORD, METODOS NUMERICOS APLICADOS A LA COMPUTACION
DEGITAL REPRESENTACIONES DE SERVICIOS DE INGENIERIA DE
EDITORES.