Solving Diff. Equation of Epidemic Model: Help Needed

Click For Summary
SUMMARY

This discussion focuses on solving differential equations related to epidemic models, specifically the compartmental models in epidemiology. The user references the Wikipedia article on deterministic models and attempts to derive the relationship between susceptible (S), infected (I), and recovered (R) individuals. The correct integration leads to the equation S(t) = S(0)exp(-Ro(R - R(0))), while the user initially derives r = 1/Ro * ln(s) + c. The conversation highlights the importance of understanding the chain rule in integration and the relationship between the variables in the model.

PREREQUISITES
  • Understanding of differential equations and their applications in epidemiology
  • Familiarity with compartmental models in epidemiology
  • Knowledge of the chain rule in calculus
  • Basic concepts of exponential functions and logarithms
NEXT STEPS
  • Study the derivation of the SIR model in epidemiology
  • Learn about the application of the chain rule in solving differential equations
  • Explore the implications of the basic reproduction number (Ro) in epidemic modeling
  • Investigate numerical methods for solving differential equations in epidemiological contexts
USEFUL FOR

Students of mathematics, epidemiologists, and researchers interested in mathematical modeling of infectious diseases will benefit from this discussion.

madiha.sahar
Messages
2
Reaction score
0
I am trying to understand this article at Wikipedia on epidemic deterministic models
http://en.wikipedia.org/wiki/Compartmental_models_in_epidemiology

it is said that if
ds/dt = -b*s(t)*i(t) is divided by dr/dt= ki(t) and is integrated by using chain rule, the answer is S(t) = S(0)exp (-Ro(Rinf - R(0))

but if i solve it, answer i get is r = 1/Ro * lns + c

i have no clue how can i get this answer.

Can anyone help me understand this please.
 
Physics news on Phys.org


Your solution may well be the inverse of the wiki solution. Ignoring the constant, if s = S(0)exp(-Ro(R - R(0)) then
ln s = ln S0 - Ro r, where r = R - R0.
ln s - ln S0 = -Ro r.
r = -ln(s/S0)/Ro plus a constant. If you normalize to S0 = 1 then r = -ln(s)/Ro + constant. I haven't solved the D.E. but I realized there is a similarity between wiki's and your solution.
 
Last edited:


I got it :) thanks a bunch!
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K