448 CAPÍTULO 8.

ECUACIONES DIFERENCIALES ORDINARIAS

∂f1


 J11,j = = l1 (m1 + m2 )
∂θ1,j+1







∂f1




 J12,j = = l2 m2 cos(θ1,j − θ2,j )
∂θ

2,j+1






+2m2 l2 (θ2,j+1 − θ2,j ) sin(θ1,j − θ2,j )





 ∂f2
J21,j = = l1 cos(θ1,j − θ2,j ) − 2l1 (θ1,j+1 − θ1,j ) sin(θ1,j − θ2,j )


∂θ1,j+1









 ∂f2
 J22,j = = l2


∂θ2,j+1

Finalmente, se pide:

a) Realizar un programa que calcule las trayectorias θ1 (t) y θ2 (t) del doble
péndulo en un intervalo de tiempo de 60s con un paso ∆t = 0,01s
para todo valor fı́sicamente factible de m1 , m2 , l1 , l2 , α1 , α2 , V1 y V2 .
Como tolerancia para el método de Newton-Raphson utilice el valor
de 1E −14 .

b) Realizar la gráfica de convergencia del método de Newton-Raphson

para los tiempos 30s y 60s. Verificar que dicha convergencia sea cua-
drática.

c) Realizar una animación de las trayectorias de las masas m1 y m2 del

doble péndulo en donde que se pinte con una lı́nea de color la trayec-
toria correspondiente a m2 a lo largo de todo el tiempo de simulación,
es decir, 60s (ver Figura 8.25).

8.12. Aplicación a la Medicina (Brote Epidémico)

El modelo SIR (Susceptibles-Infectados-Recuperados) fue utilizado por
primera vez por Kermack y McKendrick en 1927, desde entonces ha sido
aplicado a una variedad de enfermedades, en especial, enfermedades conta-
giadas por vı́a aérea en niños a los que se considera con inmunidad de por
vida luego de la recuperación, tales como: sarampión, paperas, rubéola, tos
ferina, y, en nuestros dı́as, coronavirus-19. S, I, y R representan el número de
8.12. APLICACIÓN A LA MEDICINA (BROTE EPIDÉMICO) 449

Figura 8.25: Trayectoria de la masa m2 del doble péndulo de la Figura 8,23).

individuos susceptibles, infectados, y recuperados, ası́ como N = S + I + R

es la población total de individuos a ser estudiada.

En el modelo SIR, los individuos que se han recuperados adquieren in-

munidad total con respecto al patógeno. Por otra parte, en el modelo SIRS
(Susceptibles-Infectados-Recuperados-Susceptibles), esta inmunidad se des-
vanece con el tiempo, por lo que los individuos pueden volverse a infectar.

El diagrama SIR/SIRS mostrado en la Figura 8.26 muestra como los

individuos se mueven a través de cada fase del modelo. La lı́nea punteada
muestra como el modelo SIR se convierte en un modelo SIRS, en el que
recuperarse no otorga inmunidad de por vida, por lo que los individuos
pueden volerse susceptibles nuevamente.
La tasa de infección, β, controla la velocidad de dispersión que representa
la probabilidad de transmitir una enfermedad de un individio infectado a un
individio susceptible. La tasa de recuperación, γ = 1/D, está determinada
450 CAPÍTULO 8. ECUACIONES DIFERENCIALES ORDINARIAS

β γ
Susceptibles Infectados Recuperados

Figura 8.26: Diagrama SIR/SIRS del flujo de los individuos a través de cada
fase del modelo. La lı́nea punteada muestra como el modelo SIR se convierte
en un modelo SIRS a través de ξ, la tasaa la cual los individuos recuperados
vuelven a ser susceptibles debido a la pérdida de inmunidad.

por el tiempo promedio de duración, D, de la infección. Para el modelo SIRS,

ξ es la tasa a la cual los individuos recuperados vuelven a ser suceptibles
debido a la pérdida de inmunidad.
Por otro lado, se puede agregar también la Dinámica Vital al modelo
SIRS, donde µ y ν representan las tasas de natalidad y mortalidad, respec-
tivamente. Con la finalidad de mantener una población constante, se puede
asumir µ = ν. Finalmente, note que en estado estacionario:

dI
=0
dt
Ası́, el sistema de ecuaciones diferenciales ordinarias al que se llega en el
modelo SIRS con dinámica vital, está dado por:

dS SI


 = µN − β + ξR − νS



 dt N


dI SI

=β − (γ − ν)I (8.121)


 dt N


 dR


 = γI − (ξ + ν)R
dt
donde, N = S + I + R es la población total de individuos en estudio.
8.12. APLICACIÓN A LA MEDICINA (BROTE EPIDÉMICO) 451

8.12.1. Discretización del Modelo SIRS

A continuación, se utiliza el esquema de la primera diferencia finita pro-
gresiva o delantera para la primera derivada (ver subsección 7.1.1.), con la
finalidad de discretizar el lado izquierdo del sistema de ecuaciones diferencia-
les ordinarias 8.121. A la vez, se usa un método de Euler Hacia Adelante (ver
subsección 8.4.1) para la discretización del lado derecho de las ecuaciones
8.121. Ası́, se llega al siguiente sistema de ecuaciones:

 Si+1 − Si = µNi − β Ii Si + ξRi − νSi



∆t Ni







Ii+1 − Ii Si Ii

=β − (γ + ν)Ii (8.122)


 ∆t N i


R − Ri


 i+1

 = γIi − (ξ + ν)Ri
∆t
Luego, operando en el sistema 8.122, se obtiene:

    
Ii

 Si+1 = (µNi + ξRi )∆t + 1 − β − ν ∆t Si
Ni






    
Si (8.123)
I
 i+1 = 1 + β − γ − ν ∆t Ii


 Ni




Ri+1 = [1 − (ξ + ν)∆t]Ri + γ∆tIi


que además, debe satisfacer la ecuación de equilibrio poblacional, que dis-

cretizada se re-escribe de la siguiente manera:

Ni+1 = Si+1 + Ii+1 + Ri+1 (8.124)

Note que, el sistema de ecuaciones 8.123 y 8.124 da origen a un método
explı́cito de diferencias finitas, en el cual, dado que se conoce la solución en
el tiempo ti , se puede encontrar la solución en el tiempo ti+1 .

8.12.2. Condiciones Iniciales

El sistema de ecuaciones diferenciales ordinarias 8.121 describe el com-
portamiento de un problema de valores iniciales. Es por esto, que para su
solución, es fundamental conocer ciertos valores de sus variables al inicio de
452 CAPÍTULO 8. ECUACIONES DIFERENCIALES ORDINARIAS

la epidemia.

A manera de ejercicio considere las siguientes condiciones iniciales:



 N (t = 0) = 100 personas




 S(t = 0) = 99 personas


(8.125)
I(t = 0) = 1 personas








R(t = 0) = 0 personas


Además de las condiciones iniciales del fenómeno, es necesario conocer

también los valores de las diferentes tasas involucradas en el modelo. Ası́,
para este ejercicio asumimos los siguientes valores de dichas tasa:

β = 2,2 1/dı́a










 γ = 0,1176 1/dı́a



ξ = 0,0014 1/dı́a (8.126)




µ = 0,005 1/dı́a









ν = 0,0021 1/dı́a
En la Figura 8.27, se muestra la evolución del número de individuos
susceptibles, infectados, y recuperados obtenidos al resolver el sistema de
ecuaciones diferenciales ordinarias 8.121, utilizando el esquema explı́cito de
diferencias finitas 8.123, 8.124 y con ∆t = 0,01 dı́a; bajo las condiciones
iniciales 8.125, y asumiendo los parámetros 8.126 del modelo, para un lapso
de tiempo de 90 dı́as.
Como ejercicio práctico, se pide obtener la evolución del número de in-
dividuos susceptibles, infectados, y recuperados para los mismos 90 dı́as,
asumiendo:

a) µ = ν = 0

b) µ = ν = ξ = 0

Comparar los resultados de los literales a) y b) entre estos y con respecto

a los mostrados en la Figura 8.27. Explique los cambios vistos en las curvas
8.13. EJERCICIOS 453

100
Susceptibles
Infectados
Recuperados
80

60

40

20

0
10 20 30 40 50 60 70 80 90

Figura 8.27: Evolución del número de individuos susceptibles, infectados,

y recuperados obtenidos al resolver el sistema de ecuaciones diferenciales
ordinarias 8.121 (modelo SIRS con dinámica vital), utilizando el esquema
explı́cito de diferencias finitas 8.123 y 8.124, bajo las condiciones iniciales
8.125, y asumiendo los parámetros 8.126 del modelo, para un lapso de tiempo
de 90 dı́as.

en términos de la evolución de la epidemia y a manera de escenarios, como

si se tratara de casos de la vida real.

8.13. Ejercicios
1. Un estaque se drena a través de un tubo como se observa en la Fi-
gura 8.28. Con suposiciones simplificadoras, la ecuación diferencial siguiente
describe como cambia la profundidad con el tiempo:

dh πd2 p
=− 2g(h + e)
dt 4A(h)