Modelación matemática

Una de las aplicaciones más importantes del Cálculo son las Ecuaciones Diferenciales.

Muchos procesos biológicos ocurren de manera continua con respecto al tiempo, por ejemplo
la concentración de un medicamento en el torrente sanguı́neo, el crecimiento de la masa en un in-
dividuo o las poblaciones siempre y cuando estas consideren un número de individuos muy grande.

Definición. Una ecuación diferencial es una ecuación (igualdad) que contiene una función
desconocida y una o más de sus derivadas.

Ejemplo: Sea y = y(t). Su derivada es y ′ = y ′ (t)

La igualdad

y ′ + 2y = 3
o
y ′ (t) + 2y(t) = 3

es una ecuación diferencial donde no se conoce la regla de correspondencia explı́cita de y con

respecto a t.

Definición El orden de una ecuación diferencial es el mayor oden de la derivada de y que

aparezca en la ecuación.

Ejemplos:

y = ay + bt + c es de primer orden

5y ′′ (t) − y ′ (t) = y(t) es de segundo orden

Fin del ejemplo.

Definición. La solución a una ecuación diferencial es una función tal que, al sustituirla en la
ecuación diferencial, produce una igualdad.

dy
Ejemplo: La ecuación diferencial = 2 + y(t) o y ′ = 2 + y tiene como solución la función
dt
y(t) = et − 2. Veamos:

Sustituimos y(t) = et − 2 en los miembros izquierdo y derecho de la ecuación diferencial:

1
 dy d t
Miembro izquierdo: = (e − 2)=et 1
dt dt
Miembro derecho: 2 + y=2 + (et − 2)=et 2

como 1 y 2 son iguales, y(t) = et − 2 es solución de la ecuación diferencial.

Fin del ejemplo.

Usualmente una ecuación diferencial tendrá varias soluciones y la elección de la que resuelva
el problema del que se derivó dependerá de una condición inicial que es una pareja de números
(t0 , y0 ) que cumplan con y(t0 ) = y0 .

Interpretación geométrica de las ecuaciones diferenciales

Antes de presentar los tipos y métodos de soluciones de ecuaciones diferenciales se verá la

interpretación geométrica de estas y se introducirá una téctnica llamada campo de direcciones.

Campo de direcciones

Para usar esta técnica recordaremos que y ′ = dydt

representa la pendiente de la recta tangente
a la gráfica de la función y = y(t) en el punto (x, y) dibujada en el plano t − y.

Se dibuja en el plano t − y un conjunto de pequeños segmentos de rectas tangentes determi-

nadas por la ecuación diferencial

Ejemplo: Sea y ′ = t + y

Resolución
En el miembro izquierdo de la ecuación anterior, se tiene y ′ , y el valor de esta estará determinada
por la suma de los valores t y y del mimebro derecho. Notemos que para cada pareja de valores
(t, y) se obtendrá un valor para y ′ . Por otro lado (t, y) nos indica un punto (o posición) en el
plano t − y. Esto nos permite dibujar pequeños segmentos de recta situados en la coordenada
(t, y) cuya pendiente esté determinada por t + y.

Veamos:

2
 y ′ (t, y) =t+y
y ′ (0, 0) =0+0 =0
y ′ (0, 1) =0+1 =1
y ′ (0, −1) =0−1 = −1
y ′ (0, 2) =0+2 =2
y ′ (0, −2) =0−2 = −2
y ′ (1, 0) =1+0 =1
y ′ (1, −1) =1−1 =0
y ′ (1, 1) =1+1 =2
′
y (1, −2) =1−2 = −1
y ′ (1, 2) =1+2 =2

El resultado es un campo de direcciones que se muestra a continuación:

Este campo de direcciones se asemeja a un flujo o una corriente. Las lı́neas de corriente o
flujo representarán las curvas solución y = y(t) a la ecuación diferencial. Se observa que hay
muchas soluciones que cumplen con la ecuación diferencial. La elección de una solución particular
se realiza seleccionando una condición inicial y siguiendo la única lı́nea que tiene como pendiente
uno de los segmentos de las rectas tangentes del campo de direcciones. En la el campo de
direcciones de este ejemplo, se toma la solución particular que pasa por el punto (0, 1).

Fin del ejemplo.

Tipos de ecuaciones diferenciales y soluciones

Ecuación diferencial dependiente solo de t

3
 dy
Estas ecuaciones tienen la forma = f (t) donde la derivada de y es solo una función de t
dt
pero no de y.

La solución de estas ecuaciones se obtiene por integración directa (en caso de que f (t) se
pueda integrar. Ası́, de forma general:

Z Z
dy
y= dt = f (t)dt.
dt

dy
Ejemplo: Sea la ecuación diferencial = t2 + 2. Hallar la solución más general.
dt
Resolución:

dy
= t2 + 2
Z dt Z
dy
dt = t2 + 2 dt
dt
t3
y(t) = + 2t + C
3

donde C se determina por condiciones iniciales.

Fin del ejemplo.

Ecuación diferencial autónoma.

Estas ecuaciones involucran tanto la deriavda de la función y ′ como la función misma y, pero
dy
no hay dependencia explı́cita en la variable independiente t. Toman la forma = g(y)
dt

Ejemplos:

dy
1. = 7y + 9
dt
2. Sea a un paciente al que se le está suministrando un medicamento por vı́a intravenosa a
una tasa constante A (mg/mLh). Supongamos que el cuerpo metaboliza el medicamento
a una tasa proporcional al medicamento y presente en el cuerpo en cada momento t, es
decir −ky k > 0. Bajo estas hipótesis, la situación planteada genera la siguiente ecuación
diferencial para la tasa total de medicamento está dada por:

4
 dy
= A − ky
dt

El positivo de A indica un aumento en el medicamento debido al suminstro de este y el

signo negativo de −ky indica que el cuerpo, al metabolizarlo, reduce la cantidad presente.

Esta es una ecuación diferencial autónoma ya que el miembro derecho de la ecuación solo
depende de y.

Fin del ejemplo.

Ecuación diferencial no autónoma

Estas ecuaciones combinan las dos anteriores:

dy
Ejemplo: = 1 + sen t − y
dt
Se observa que en el miembro derecho de la igualdad combina una función de y y otra de t,
es por eso que es una ecuación diferencial no atónoma.

Fin del ejemplo.

Ecuaciones diferenciales separables

dy
Es una ecuación diferencial de primer orden en que la expresión puede ser factorizada
dt
como una función de t por una función de y, es decir
dy
= f (t)g(t)
dt
Para mostrar la técnica general de resolución usaremos la notación de diferenciales:

5
 dy
= f (t)g(y)
dt
1 dy
= f (t)) si g(y) ̸= 0
g(y) dt
1
dy = f (t)dt donde se usó la notación de diferenciales
g(y)
1
h(y)dy = f (t)dt renombrando h(y) =
g(y)
Z Z
h(y)dy = f (t)dt integrando de ambos lados

El miembro izquierdo de la última igualdad solo depende de y mientras que el miembro dere-
cho solo depende de t. Si son funciones integrables, se realiza la integral y, de ser posible, se
despeja y en función de t.

Caso particulares:
dy
Ecuación del tipo: = ky k =constante.
dt
Si y = 0 para toda t entonces la solución de la ecuación es y = C donde C es una constate.

Si y ̸= 0 entonces

dy
= kdt
dt
dy
=k dividiendo entre y
y
ln |y| = kt + C integrando respecto a las variables indicadas

|y| = ekt+C = eC ekt despejando y

y = Aekt donde A = ±eC

A continuación se muestran los campos de direcciones para valores de k > 0 y k < 0.

6
Aplicación: Modelo de crecimiento poblacional de Malthus

7
 Este modelo de crecimiento poblacional considera las siguientes hipótesis:

ˆ No hay limitantes de espacio.

ˆ No hay limitantes de alimento.

ˆ No hay predadores.

ˆ Hay inmunidad a enfermedades.

ˆ Se supone que la población crece con una tasa que es proporcinoal al tamaño de la población
presente al tiempo t.

La última hipótesis es la que contiene la información que generará la ecuación diferencial de

este modelo:

Sea N el número de individuos de la población al tiempo t. Entonces

dN
= kN k > 0 constante de proporcionalidad
dt
La solución a esta ecuación diferencial es

N (t) = N0 ekt
con N0 = N (0) es la población inicial y k > 0 para indicar que la población N crece conforme
al tiempo, haciéndolo de forma exponencial.

Comprobación:

dN d
= N0 ekt = N0 kekt = kN0 ekt = kN
dt dt
dN
dN
La ecuación diferencial = k es equivalente a dt = k que se interpeta como la tasa de
dt N
crecimiento per cápita.

Este modelo de crecimiento tiene hipótesis muy restrictivas y que no modelan a largo plazo
el crecimiento de una población; sin embargo, modela bien los primeros momentos del establec-
imiento de un conjunto de individuos creciendo de forma exponencial.

Ecuación del tipo y ′ = ay + b con a, b =const.

Si b = 0 estamos en el caso anterior.

8
 dy
Sı́ b ̸= 0, entonces tenemos la ecuación = ay + b.
dt
Resolución:

dy b
= a(y + )
dt a

dy b
= adt dividiendo entre y +
y + ab a
Z Z
dy
ds = adt
y + ab

b
ln |y + | = at + C integrando respecto a las variables indicadas
a

b
|y + | = eat+C = eC eat despejando y
a

b
y+ = Aeat donde A = ±eC
a
Ası́ que la solución es:

b
y = Aeat −
a

Aplicación: Crecimiento restringido

Hipótesis:

ˆ Ningún organismo o población crece indefinidamente (por falta de alimento, condiciones

fı́sicas desfavorables, etc.)

ˆ Hay un lı́mite superior para N denominado B al cual se aproxima asintóticamente, es decir,

dN
dt
tiende a cero cuando B − N se hace pequeña (B > N ).

La ecuación diferencial que se genera debido a estas hipótesis se escribe como sigue:
dN
= k(B − N ) k>0
dt
Observaciones:

9
 ˆ Si N << B (N mucho menor que B), entones N ′ ≈ kB que es una constante, lo que
implica que N tiene un crecimiento que se aproxima a una lı́nea recta de pendiente kB.

ˆ Si N ≈ B entonces 0 < B − N << 1 que es positiva pero muy cercana a cero lo que
implica que 0 < N ′ << 1 que indica que N ′ ≈ 0 y queire decir que N , al acercarse al valor
de B tiene una tasa de crecimiento que se aproxima a cero.

Reescribiendo la ecuación diferencial en la forma

dN
= −kN + kB
dt
observamos que tiene la forma de y ′ = ay + b con y = N , a = −k y kB = b

Por lo tanto, la solución a esta ecuación diferencial es

N (t) = Ae−kt + B
como y < B entonces A debe ser negativa.

Este modelo puede ser usado para modelar tamaños de cosechas, crecimiento de un individuo
o un tejido, etc., cuyo crecimiento en etapas tempranas es lineal y se aproxima asintóticamente
a un valor fijo B para etapas posteriores.

Ecuación del tipo y ′ = ay 2 + by + c

Observemos que ay 2 + by + c es un polinomio de grado 2. Supondremos que tiene 2 raı́ces

reales distintas A y B (que se pueden hallar usando la fórmula general), ası́:

dy
= ay 2 + by + c
dt
dy
= a(y − A)(y − B)
Z dt Z
dy
= a dt
(y − A)(y − B)

La integral en el miembro izquierdo de la ecuación se resuelve por el método de las fracciones

parciales:

10
 1 C D
= +
(y − A)(y − B) (y − A) (y − B)

Cy − CB + Dy − AD (C + D)y − CB − AD
⇔ =
(y − A)(y − B) (y − A)(y − B)

Lo que genera el siguiente sistema de ecuaciones:

⇔C +D =0
−CB − AD = 1

cuya solución es

1 1
D= C=−
B−A B−A

De esta forma se encuentran las fracciones parciales

1 1 1
=− +
(y − A)(y − B) (B − A)(y − A) (B − A)(y − B)

 
1 1 1
= −
(B − A) y−B y−A

11
Ası́, resolviendo le ecuación diferencial
Z Z
dy
= a dt
(y − A)(y − B)

Z   Z
1 1 1
− dy = a dt
(B − A) y−B y−A

Z   Z
1 1
− dy = a(B − A) dt
y−B y−A

Z Z
dy dy
− = a(B − A)t + C
y−B y−A

ln |y − B| − ln |y − A| = a(B − A)t + C

y−B
ln = a(B − A)t + C
y−A

y−B
= ea(B−A)t+C
y−A

Haciendo k = ±eC

y−B
= −kea(B−A)t
y−A

Resolviendo para y

y − B = (y − A)(−k)ea(B−A)t

y + ykea(B−A)t = B + Akea(B−A)t

y[1 + kea(B−A)t ] = B + Akea(B−A)t

B + Akea(B−A)t
y=
1 + kea(B−A)t

Donde B y A son las raı́ces del polinomio de segundo grado, a es el coeficiente del término
cuadrático de la ecuación diferencial y k se fija por condiciones iniciales.

12
 Aplicación: Ecuación logı́stica de crecimiento restringido

Esta ecuación combina el modelo de crecimiento de Malthus y el de crecimiento restringido

para conseguir un crecimiento exponencial en etapas tempranas del establecimiento de una
población y que, al paso del tiempo, el número de habitantes no rebase de cierto lı́mite. La
dN
ecuación diferencial que se obtiene de estas consideraciones se obtiene considerando que es
dt
proporcional tanto a N (Malthus) como a (B − N ) (crecimiento restringido). De esta forma se
tiene
dN
= λN (B − N ) λ>0
dt
Reescribiendo el miembro derecho de la igualdad:
dN
= λN (B − N ) = −λN 2 + λN B = −λ(N − 0)(N − B)
dt
Al reescribirla de esta forma, se observa que esta ecuación diferencial tiene una forma cuadrática
del tipo y ′ = ay 2 + by + c = a(y − A)(y − B). Identificando los valores A = 0, B = B y a = −λ
se tiene la solución dada por:

B + (0)ke−λ(B−0)
N=
1 + ke−λ(B−0)t

B
N=
1 + keλBt
Donde k se fija por condiciones iniciales.

Analicemos el lı́mite de N cuando t → ∞

B B
lim N = lim −λBt
= =B
t→∞ t→∞ 1 + ke 1 + k(0)
Ası́, para valores muy grandes de t, N se aproxima asintóticamente a B que es lo que es-
perábamos.

Para analizar la forma de las soluciones, se realizará un campo de direcciones para la ecuación
dN
diferencial = N (2 − N ) donde se ha fijado λ = 1 y B = 2.
dt

13
 En la gráfica se muestra la solución gráfica para un valor inicial de 0.1.

Esta forma de crecimiento modela muy bien el crecimiento de un cultivo de levadura como
es el que se muestra a continuación:

En la siguiente gráfica se muestran varias soluciones a la ecuación diferencial cuya población

inicial es menor que B. En todas ellas se marca un punto de inflexión

14
 El punto de inflexión, de forma simplificada, puede ser hallado al hallar los puntos para los
cuales la segunda derivada de N es igual a 0. Un método para hallar dónde ocurren los puntos
de inflexión es el siguiente:

dN
= λN (B − N )
dt

Derivando de ambos lados con respecto a t

d2 N
 
dN dN
=λ (B − N ) + N (− )
dt2 dt dt
dN
=λ (B − N − N )
dt
dN
=λ (B − 2N )
dt

Igualando a cero

dN
λ (B − 2N ) = 0
dt

15
Lo que ocurre si y solo si

B = 2N

dN
ya que dt
> 0 en el punto de inflexión.

De esta forma, el punto de inflexión ocurre cuando N = B/2, esto es, el punto de inflexión
se localiza cuando N alcanza la mitad de la población lı́mite B.

Sistema predador-presa (Modelo Lotka-Volterra)

En general, las poblaciones interaccionan con otras especies y el modelo que resulta de con-
siderar una interacción genera un sistema de ecuaciones diferenciales.

El sistema modela dos especies cuya cantidad de individuos depende una de la otra. El sistema
tendrá entonces dos funciones desconocidas dependientes del tiempo, en el caso que veremos,
una será el predador y el otro la presa.
Hipótesis del modelo:

ˆ Si no hay predadores presentes, las presas se reproducen con una tasa proporcional a su
población y no les afecta la sobrepoblación.

ˆ Los predadores se comen a las presas y la razón a la que las presas son devoradas es
proporcional a la interacción de presas y predadores.

ˆ Sin presas qué comer, la población de predadores declina a una razón proporcional a ella
misma.

ˆ La tasa de nacimiento de predadores va en proporción al número de presas comidas

por predador que, por la segunda hipótesis, es proporcional a la interacción del presas
y predadores.

Denotamos C a las presas (conejos) y L a los predadores (lobos). Entonces, por la primera
hipótesis, en ausencia de predadores se tendrı́a que la tasa de crecimiento para los conejos es
dC
|
dt L=0
= kC con k > 0 y de la tercera hipótesis, en ausencia de presas la tasa de crecimiento
para los lobos es dL |
dt C=0
= −rL con r > 0 donde el signo menos indica que la población de lobos
decrece en ausencia de presas.

Con la presencia de las dos epecies, la interacción se mide con el producto de número de
presas por el número de predadores, es decir CL.

16
 De lo anterior, tenemos dos ecuaciones:
dC
= kC − aCL k, a > 0 es la tasa de cambio total para las presas
dt
dL
= −rL + bCL r, b > 0 es la tasa de cambio total para los predadores
dt
Las dos ecuaciones anteriores forman un sistema de ecuaciones diferenciales no lineales de
primer orden y son llamadas ecuación predador-presa o Lotka-Volterra.

Una solución a este tipo de sistemas de ecuaciones diferenciales es una pareja de funciones
C(t) y L(t) que describe la población de presa-predador en función del tiempo. Pero como es
un sistema acoplado, no se puede hallar primero una ecuación y luego otra, se deben resolver si-
multáneamente. Desafortunadamente generalmente es imposible hallar una solución con fórmulas
para C y L pero se pueden usar métodos gráficos para su estudio.

Ejemplo: Tomaremos los valores k = 0.08, a = 0.001, r = 0.02 y b = 0.00002 para el sistema
de ecuaciones anterior, de esta forma

dC
= 0.08C − 0.001CL
dt
dL
= −0.02L + 0.00002CL
dt
Primero analizaresmos las soluciones de equilibrio que son la pareja C y L para las cuales
dC
dt
= 0 y dL
dt
= 0 que nos deja el siguiente sistema de ecuaciones:

C(0.08 − 0.001L) = 0
L(−0.02 + 0.00002C) = 0

La cual se soluciona si R = 0 y L = 0 (que es la solución trivial) o si

0.08
L= = 80
0.001

0.02
C= = 1000
0.00002
Este resultado quiere decir que si hay 80 lobos y 1000 conejos, la polación se encuentra en
equilibrio y ası́ se queda al transcurrir el tiempo.

El siguiente paso es analizar gráficamente el comportamiento de un campo de direcciones en

el plano L − C. Para esto, usaremos la regla de la cadena para hallar una expresión para dL
dt
, ası́:

17
 dL dL dC
=
dC dC dL

dL
dL dt −0.02L + 0.00002LC
=⇒ = dC
=
dC dt
0.08C − 0.001LC

El campo de direcciones se muestra a continuación:

las cuvras solución parten de un punto y realizan un ciclo regresando al mismo punto del que
partieron. Por ejemplo, una curva puede pasar por el punto (1000,40).

dC
Para determinar la dirección de movimiento de la curva, se analiza el signo de dt
en el punto
(1000,40):

dC
= 0.08(1000) − 0.001(1000)(40) = 80 − 40 = 40
dt

18
 dC
ası́ que = 40 > 0 en el punto (1000,40), por lo que indica que la curva se mueve
dt
en dirección de las manecillas del rejo, porque, para ese punto, la población de conejos debe
incrementarse.

Fin del ejemplo.

19