Differential equation for exponential growth

Click For Summary
SUMMARY

The discussion centers on modeling the exponential growth of viruses against a constant elimination rate of 50,000 viruses per hour. The initial population of viruses is set at 10^6, with a growth rate calculated as k = ln(1.6)/240, reflecting a 160% increase over four hours. The differential equation is established as dP/dt = P*k - 833.33, where 833.33 represents the elimination rate in viruses per second. Participants emphasize the need to integrate this equation correctly to determine the time when the virus population reaches 10^12.

PREREQUISITES
  • Understanding of differential equations
  • Familiarity with exponential growth models
  • Knowledge of integration techniques
  • Basic concepts of population dynamics
NEXT STEPS
  • Study the integration of differential equations in population models
  • Learn about the implications of constant elimination rates in exponential growth
  • Explore the application of the logistic growth model in biological systems
  • Investigate the use of numerical methods for solving differential equations
USEFUL FOR

Students in mathematics or biology, researchers in virology, and anyone interested in mathematical modeling of population dynamics will benefit from this discussion.

nrslmz
Messages
14
Reaction score
0

Homework Statement


Viruses are reproducing exponentially, while the body eliminates the viruses.
The elimination rate is constant, 50000 per hour. I decided to take down on the minute level, so it would be 50000/60.
Pinitial is 10^6
k is ln(1,6)/240, since the growth rate is 160% in 4 hours.
Then the exponential and differential equation would be:
P(t) = 10^6 * (ln(1,6) * t / 240)




Homework Equations


Then the exponential and differential equation would be:
P(t) = 10^6 * (ln(1,6) * t / 240)



The Attempt at a Solution


dP/dt = P*k ..., right? How can I add up the elimination rate to the differential equation.
And I should integrate this to find the specific moment when the population reaches 10^12? I am somehow not feeling good about this solution. I would appreciate if you help.
 
Physics news on Phys.org
Look at the units of dP/dt. P is number of virus, t is time, in seconds. So dP/dt has units of "viruses per second". The right hand side of the equation dP/dt= ... must also have units of "viruses per second" and tells how the number of viruses changes. There are two reasons why that would change:
1) The viruses reproduce at a rate proportional to their number.

2) Some viruses are killed- "50000 per hour" so -5000 viruses per hour has precisely the correct units. Just subtract 5000 from the right side.
 
Sorry, I posted wrong again. The function is P(t) = 10^6 * e^(ln(1,6) * t / 240)
You mean 50000/60 = 833,33? and
Should I write:
dP/dt = P*k -833,33 and than integrate? Bu that would be the same as
P(t) = 10^6 * e^(ln(1,6) * t / 240) - 833,33t
I don't think that would work. Or should I take the limit of
P(t) = 10^6 * e^(ln(1,6) * (t/n) / 240) - 833,33(t/n) when n goes to infinity?
 
nrslmz said:
Sorry, I posted wrong again. The function is P(t) = 10^6 * e^(ln(1,6) * t / 240)
You mean 50000/60 = 833,33? and
Should I write:
dP/dt = P*k -833,33 and than integrate?
You can use 50k directly. As Halls said there is no need to do it per minute.

Bu that would be the same as
P(t) = 10^6 * e^(ln(1,6) * t / 240) - 833,33t
I don't think that would work. Or should I take the limit of
P(t) = 10^6 * e^(ln(1,6) * (t/n) / 240) - 833,33(t/n) when n goes to infinity?
What is n supposed to be here? And what did you work k out to be?
 
Wouldn't the operation be
((( * e^(kt))-50000)* e^(kt))-50000)-50000...

and e^(kt) must be 1.6 since it grows into %160 of its initial population in 4 hours. Therefore t is 4 or 240 depending on the time interval.

I solved e^(k * 240) = 1,6 to find k.
 
nrslmz said:
Sorry, I posted wrong again. The function is P(t) = 10^6 * e^(ln(1,6) * t / 240)
You mean 50000/60 = 833,33? and
Should I write:
dP/dt = P*k -833,33 and than integrate? Bu that would be the same as
P(t) = 10^6 * e^(ln(1,6) * t / 240) - 833,33t
No, it would not. You are not integrating the equation correctly.

I don't think that would work. Or should I take the limit of
P(t) = 10^6 * e^(ln(1,6) * (t/n) / 240) - 833,33(t/n) when n goes to infinity?
 

Similar threads

Replies
4
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 7 ·
Replies
7
Views
4K
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
4K
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 10 ·
Replies
10
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K