“AÑO DEL BICENTENARIO DEL PERÚ:200 AÑOS DE INDEPENDENCIA”

UNIVERSIDAD NACIONAL
AGRARIA LA MOLINA

Curso: Cálculo para Ciencias

Profesores: Esaul Obregon Parraga / Ivan Solis Saavedra Mendoza

Título: Aplicaciones de ecuaciones diferenciales a la Biología

INTEGRANTES:

● Mogollón Ampuero, Gabriel Alejandro 20200142

● Ramos Céspedes, Joaquín Matías Sebastián 20200146
● Salcedo Montoya, Cristopher Alexander 20200148
● Tasayco Flores, Alvaro Luis 20200150

Octubre- 2021
 ÍNDICE

Aplicación 1: Datación por radiocarbono 2

Aplicación 2: Dinámica de crecimiento en un individuo: modelo de Bertalanffy 5

Aplicación 3: Modelo epidemiológico 7

Aplicación 4: Leyes de Kirchhoff en circuitos RLC 11

Referencias Bibliográficas 15

1
Aplicación 1: Datación por radiocarbono

Esta técnica se utiliza para determinar la edad de objetos compuestos de sustancias

orgánicas, la Datación por radiocarbono está basada en la Ley de crecimiento y
decrecimiento exponencial.

Esta técnica puede utilizar el fechado por carbono-14, su fundamento radica en una
reacción química de gran importancia, que ocurre en la atmósfera terrestre, tal como plantea
Chang (2003): “Una de las reacciones importantes entre la atmósfera y los rayos cósmicos es
la captura de neutrones por el nitrógeno atmosférico (el isótopo nitrógeno-14) para producir
el isótopo radiactivo carbono-14 e hidrógeno…”.

El carbono-14 es producido de forma continua en la atmósfera que luego formará parte

del dióxido de carbono (CO2), al combinarse con el Oxígeno, para ser asimilado por las
plantas que a su vez son ingeridas por los animales.

Los átomos de carbono-14 presentes en los seres vivos están en constante desintegración,
pero simultáneamente, son reemplazados por nuevos átomos a un ritmo constante, de modo
que la proporción de carbono-14 presente en cualquier organismo es similar a la de la
atmósfera, sin embargo cuando un organismo muere deja de absorber carbono-14, debido a
esto al transcurrir los años la proporción de carbono-14 presente en las sustancias orgánicas
tiende a disminuir.

La técnica llamada del 𝐶14 , para datar una sustancia consiste en medir la cantidad de 𝐶14
que queda en la actualidad en dicho objeto, y utilizar la forma de las soluciones de la
ecuación de decrecimiento exponencial para calcular el tiempo que ha transcurrido.

Se sabe que la velocidad de desintegración de una sustancia radiactiva (es decir, el

número de átomos que se desintegran por unidad de tiempo) en un instante dado es
proporcional al número de átomos existentes en ese instante. En consecuencia, si se denota
por A(t) el número de átomos de la sustancia original presentes en el instante t, se puede
escribir:

A´(t) = −λA(t),
donde el signo menos indica que la velocidad es negativa (el número de átomos disminuye) y
la constante de proporcionalidad, λ > 0, se llama constante de descomposición o de
decrecimiento, y es propia de cada sustancia radiactiva.

Si se conoce el número de átomos presentes en un instante dado, por ejemplo se sabe que
en t = 0 es A(0) = 𝐴0, y se conoce también la constante de decaimiento, λ, entonces se puede
predecir el número de átomos presentes en cualquier instante posterior, ya que A(t) es la
solución del problema de valor inicial:

2
Podemos resolver esta ecuación lineal de primer orden como una ecuación de Bernoulli.

Una ecuación diferencial lineal de primer orden es de la siguiente forma:

Para llegar a la solución general explícita resolvemos la ecuación de la siguiente manera:

Aplicando estas fórmulas a la ecuación diferencial para la desintegración radiactiva se

obtiene:

3
Ejemplo
Si se encuentra un fósil que tiene un 1% de la cantidad original presente en la atmósfera
de C-14. Determinar la edad de dicho fósil si se sabe que el C-14 tiene una vida media de
5730 años.

1. Calculamos el valor de la constante de descomposición (λ) para el C-14 utilizando el

dato de la vida media (Vm).

2. Como la proporción del C-14 es el 1% de la concentración original, eso significa que:

4
Esto quiere decir que el fósil tiene una edad aproximada de 38 059 años

Aplicación 2: Dinámica de crecimiento en un individuo: modelo de Bertalanffy

"En los años 50 del siglo XX, el biólogo austriaco L. von

Bertalanffy (1901-1972) desarrolló un modelo matemático para
la talla de un individuo en función de su edad, que se utiliza
con frecuencia para predecir el tamaño de los peces" .
(Ecuaciones diferenciales aplicadas a la Biología, s. f. , p23)

Para ello se estableció lo siguiente:

1. La velocidad con la que crecen los individuos no es la misma a lo largo del tiempo
2. Al inicio la velocidad de crecimiento es muy grande; sin embargo, con el aumento del
tamaño y la madurez sexual, su rapidez va disminuyendo hasta ser prácticamente
cero. Dejando de crecer.

Figura 2: Curva de crecimiento típica de los peces.

Se puede imaginar el comportamiento de la gráfica
gracias a los principios 1. y 2

Fuente: Tomado de El modelo de crecimiento de Von

Bertalanffy, 2011.

5
 Entonces la ley de crecimiento, nos dice que la velocidad de crecimiento es proporcional
a la diferencia entre la longitud maxima y la actual.
𝐿´(𝑡) = 𝑘(𝐴 − 𝐿(𝑡))

- “L´(t)” la velocidad de crecimiento en el instante t

- k > 0, la constante de proporcionalidad de cada especie.
- "L(t)" la longitud del individuo en la edad t
- "A" la talla máxima de la especie
Ahora le sumamos una condición. "Si en el instante inicial, t = 0, la longitud del individuo es
0 < L0 < A, entonces L0, es la talla en el instante 0, será solución del siguiente problema de
valor inicial":
𝐿´ = 𝑘(𝐴 − 𝐿) .... α
𝐿(0) = 𝐿0 .... β
Entonces vemos que la ecuación diferencial α es de variables separables:

𝑑𝐿 𝑑𝐿
𝑑𝑡
= 𝑘(𝐴 − 𝐿) ⇔ ∫ 𝐴−𝐿
= ∫ 𝑘𝑑𝑡 ⇔ − 𝑙𝑛|𝐴 − 𝐿| = 𝑘𝑡 + 𝐶

Como sabemos 0 < 𝐿(𝑡) ≤ 𝐴

−𝑘𝑡
− 𝑙𝑛(𝐴 − 𝐿) = 𝑘𝑡 + 𝐶 ⇔ 𝐴 − 𝐿 = 𝐶𝑒
Por lo que la solución general es:
−𝑘𝑡
𝐿 = 𝐴 + 𝐶𝑒 , 𝐶 ∈ 𝑅
Con la condición inicial β:
−𝑘0
𝐿0 = 𝐿(0) = 𝐴 + 𝐶𝑒 ⇔ 𝐿0 − 𝐴 = 𝐶
Así la solución del problema del valor inicial. Que viene a ser la longitud de la especie en un
instante “t”

−𝑘𝑡
𝐿(𝑡) = 𝐴 + (𝐿0 − 𝐴)𝑒

Ejemplo
A partir de la ecuación podemos plantearnos un caso:
A = 50 (cm) , K = 0.5 (1/año) , L(0) = 0 cm
Reemplazando
−0.5𝑡
𝐿(𝑡) = 50 + (0 − 50)𝑒
¿Como se comporta la ecuación cuando el tiempo tiende al ∞?
−0.5𝑡 50
lim 50 − 50𝑒 = lim 50 − ∞ = 50 = 𝐴
𝑡→∞ 𝑡→∞ 𝑒

Entonces en la gráfica vemos una asíntota “A” que como indica la máxima longitud de la
especie nunca se pasará de este.

6
 Figura 3: Curva de crecimiento de acuerdo
a los parámetros A = 50 (cm) ,
K = 0.5 (1/año) , L(0) = 0 cm

Fuente: Tomado de Ecuaciones diferenciales

aplicadas a la Biología, s. f.

Aplicación 3: Modelo epidemiológico

El siguiente modelo es aplicado para poder hallar la difusión de una enfermedad

contagiosa y como va aumentado el número de infectados con el tiempo.Para el siguiente
proceso debe establecerse o plantearse algunas hipótesis:
● La población será un número fijo el cual será representado por P y sabes que cada
miembro de la población será susceptible a la enfermedad.
● La enfermedad tendrá una duración larga, la cual no se podrá curar mientras esté en
proceso o en periodo de estudio.
● Todos los individuos en el área infectada serán contagiosos y circulan libremente
entre la población.
● Durante pase el tiempo cada persona que esté infectada tendrá contactos el cual se
denotará por la letra C y cada contacto con una persona que no está infectada
redundará en la trasmisión de la enfermedad.
Una vez se tiene las hipótesis o aclaraciones ,vamos a considerar un periodo corto de
tiempo que irá desde t hasta t+h.Y cada persona infectada tendrá ch contactos de ello
¿Cuántos de estos contactos son personas no infectadas?. Si establecemos a f(t) como en
número de personas infectadas en un tiempo t, entonces se sabrá que P- f(t) es el número de
personas que no estarán infectadas, y que (P-f(t))/P será la fracción de la población que no
está infectada.Entonces podemos llegar que la cantidad de ch contactos hechos por alguien
infectado sería, los cuales habrán sido con personas no infectadas.

7
El número total de de nuevas infecciones deberá ser,

ahora dividiéndolo entre h, y haciendo que h tienda a cero obtenemos

del cual notamos que es un EDO de primer orden y es el caso para aplicar la ecuación de
Bernoulli

Donde K=2 entonces aplicando un cambio de variable para transformar a una ecuación lineal
decimos que:

8
 De la ecuación lineal tenemos que P(x)=c y Q(x)=c/P,para luego hallar el factor de
integración y resolver como una ecuación lineal:

resolviendolo la función f verifica la ecuación la ecuación diferencial que dará lugar a la

ecuación logística de:

donde c y K se podrán determinar las características de la epidemia.

Ejemplo:
Los servicios de salud pública registran la difusión de una epidemia de gripe de duración
particularmente larga en una ciudad de 500.000 personas. Al inicio de la primera semana de
registro se habían contabilizado 200 casos; durante la primera semana aparecieron 300
nuevos casos. Nos proponemos estimar el número de individuos infectados después de 6
semanas.

De los datos ya mencionados podemos sacar que el valor de P que aparece es la

capacitancia del sistema, que será en nuestro caso el número de individuos que a largo plaza
se llegaran a infectar, P=500000.También tenemos que para el t = 0, f(0)=200 y sustituyendo
en la ecuación,

9
obtenemos que el valor de K=2449. Cómo se reportaron que en la primera semana
aparecieron 300 nuevos casos,tendríamos que en total habrían 500 casos la primera
semana,pudiendo ahora aplicar la ecuación para:

Obteniendo el valor de c en consecuencia:

De ello obtenemos que el número de personas infectadas al final de la sexta semana

será 44646 aproximadamente.

Representación gráfica de

10
Aplicación 4: Leyes de Kirchhoff en circuitos RLC

Supongamos tener el siguiente circuito RLC

Figura 4. Circuito RLC.

Fuente. Elaboración propia.

Observamos que el circuito está formado por un batería conectada en serie a una resistor
de resistencia (R) un condensador de capacitancia C y un inductor de inductancia L. La
corriente generada por la batería y por ende la carga variará con respecto al tiempo a
consecuencia de la descarga y carga del condensador. Por segunda Ley de Kirchoff sabemos
que: “La suma de las diferencias de potencial a través de todos los elementos alrededor de
cualquier espira de un circuito cerrado debe ser igual a cero” (Serway R. y Jewett J. ,2008 ,
p.785) , entonces:

V + V R + V L + VC = 0
Sabemos que la resistencia , el condensador y el inductor generan una caída de voltaje:

𝑑𝐼
VR = - I R VL = - L VC = - Q/C
𝑑𝑡

Reemplazando:
𝑑𝐼
V - I R - L( ) - Q/C = 0
𝑑𝑡

𝑑𝐼
V = I R + L( ) + Q/C ….( 1 )
𝑑𝑡

11
 2
𝑑𝑄 𝑑𝐼 𝑑𝑄
sabemos que I = entonces = 2 , reemplazando en (1)
𝑑𝑡 𝑑𝑡 𝑑𝑡

2
𝑑𝑄 𝑑𝑄 1
V=L 2 +R + Q -----> EDO 2do orden
𝑑𝑡 𝑑𝑡 𝐶

Si derivamos la ecuación 1 con respecto al tiempo:

2
𝑑𝑉 𝑑𝐼 𝑑𝐼 𝑑𝑄 1 𝑑𝑄
= R +L( 2 )+( ) ; reemplazando I =
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝐶 𝑑𝑡

2
𝑑𝑉 𝑑𝐼 𝑑𝐼 1
= L( 2 )+ R + I -----> EDO 2do orden
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝐶

Obtenemos una ecuación diferencial de segundo orden homogénea. Al

resolverla obtenemos una función que nos brinda la corriente en
cualquier instante t

EJEMPLO:
Se conecta en serie una fuente de voltaje V = 1.5 V, una resistencia R = 20 Ω, un
capacitor de 10-3 F y un inductor L = 0.1 H . Determinar la carga en el capacitor y la corriente
que circula por el circuito en todo tiempo, si inicialmente el capacitor está totalmente
descargado y no fluye corriente sobre el circuito. ( Universidad Autónoma Metropolitana ,
p.2)

Gráfica:

Figura 5. Circuito RLC.

Fuente. Elaboración propia.

12
Datos : Capacitancia ( C ) = 10-3 F Batería (V) = 1.5 V Resistencia = 20 Ω
Inductancia (I) = 0.1 H

Condiciones iniciales : Q(0) = 0 y I (0 ) = 0

Primero hallamos la carga en función del tiempo usando la ecuación diferencial

obtenida al inicio.

2
𝑑𝑄 𝑑𝑄 1
V=L 2 +R +Q
𝑑𝑡 𝑑𝑡 𝐶

2
𝑑𝑄 𝑑𝑄
1.5 = 0.1 2 + 20 + 103Q
𝑑𝑡 𝑑𝑡

*Multiplicando por 10

2
𝑑𝑄 𝑑𝑄
15 = 2 + 200 + 104Q
𝑑𝑡 𝑑𝑡

Obtenemos una ecuación diferencial de segundo grado no homogénea , entonces su

solución general está compuesta por la solución asociada a la homogénea y la
solución particular.

Q(t) = Qh + Q p

Encontrando la solución asociada a la homogénea Qh:

Ecuación característica : λ2 + 200 λ + 104 = 0 ; λ1 = -100 y λ2 = -100

Cómo obtenemos raíces iguales la solución de la homogénea es:

−100𝑡 −100𝑡
Qh = C1 𝑒 + C2𝑡𝑒

Encontrando la solución particular Qp:

Como K(x) = 15 ; asumo una posible solución Qp = A , siendo A una constante.

Derivando Qp = A Q p’ = 0 Q p’’= 0 , reemplazando en la ecuación

diferencial.
15 = 0 + 200*0 + 104A

13
 15
Qp = A =
10000

Entonces la solución general es :

−100𝑡 −100𝑡 15
Q(t) = C1 𝑒 + C2𝑡𝑒 + 10000
𝑑𝑄
Si derivamos respecto al tiempo Q(t) obtenemos la ecuación de corriente porque I =
𝑑𝑡

𝑑𝑄 −100𝑡 −100𝑡 −100𝑡

𝑑𝑡
= -100C1 𝑒 + C2𝑒 - 100 C2 t 𝑒

−100𝑡 −100𝑡 −100𝑡

I(t) = -100C1 𝑒 + C2𝑒 - 100 C2 t 𝑒

Con las condiciones iniciales Q(0) = 0 y I (0 ) = 0 hallamos C1 y C2 :

15
Q(0) = 0 = C1 + 10000
I(0) = 0 = -100C1+ C2

Resolviendo : C1 = -0.0015 C2 = -0.15

Finalmente reemplazando:

−100𝑡 −100𝑡 −100𝑡

I(t) = 0.15 𝑒 - 0.15𝑒 +15 t 𝑒
−100𝑡 −100𝑡 15
Q(t) = - 0.0015 𝑒 - 0.15 𝑡𝑒 +
10000

Respuesta:

−100𝑡
I(t) = 15 t 𝑒
−100𝑡 −100𝑡 15
Q(t) = - 0.0015 𝑒 - 0.15 𝑡𝑒 +
10000

14
Referencias Bibliográficas

● Aplicación de las Ecuaciones Diferenciales Ordinarias en la Datación de Fósiles

mediante el Isótopo Radiactivo Carbono-14 y su resolución mediante el Método de
Runge-Kutta de 4° orden — Steemit. (2018). Steemit.com; Steemit. Aplicación de las
Ecuaciones Diferenciales Ordinarias en la Datación de Fósiles mediante el Isótopo
Radiactivo Carbono-14 y su resolución mediante el Método de Runge-Kutta de 4°
orden — Steemit

● ‌Chang , R. (2003) Química 7° Edición, Editorial McGraw-Hill pág. 530

● Ecuaciones diferenciales aplicadas a la Biología 2.1 Introducción. (s.f.).

http://departamento.us.es/edan/php/asig/LICBIO/LBII/Teoria2BIOII0910.pdf

● El modelo de crecimiento de Von Bertalanffy (2011) . Deducción del modelo de

crecimiento de Von Bertalanffy. Recuperado de
https://estadistica-dma.ulpgc.es/FCC/pdf/Bertalanffy.PDF

● Medina J . (2005) Tema 16. Ecuaciones diferenciales. Recuperado de:

C:\JUAN\...\OPENCALC\ecdif.DVI (upct.es)

● Modelos basados en E.D.O. (s. f.). . Recuperado 10 de octubre de 2021, de teoria

continuo tema2.pdf (ujaen.es)

● Serway R. y Jewett J. (2008) .Física para ciencias e ingeniería con física

moderna.Volumen 2. Séptima edición.

● Universidad Autónoma Metropolitana (s.f.) Aplicaciones de ED de segundo orden .

Recuperado de: ImpRL Continua.DVI (uam.mx)

15