Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Population Growth

  1. Dec 8, 2008 #1
    I am trying to model the growth of a population which replicates at a rate of 160% every four hours. Also 50 000 members die every hour. so t denoting time in hours, P denoting the population, k being 0.4 I wrote:

    dP/dt=kp-50 000 *t

    Is my approach to solving the differential faulty?

    dP/dt = kP- 50 000*t

    P=kPt - (50 000*t^2)/2

    I don't know very much about differential equations so I might be doing really really wrong...

    Thanks in advance
  2. jcsd
  3. Dec 8, 2008 #2


    User Avatar
    Homework Helper

    Why do you mutiply it by t? If you do so, you are saying that for every hour which elapses, the population does not decrease constantly but incrementally, ie. more more die per hour as time elapses.
  4. Dec 8, 2008 #3
    Hmm... You're right...

    So would this be a correct equation for the problem I stated?

    Pnew = P0ekt - 50000*t
  5. Dec 8, 2008 #4


    User Avatar
    Homework Helper

    I think you forgot to divide 50,000 by k
  6. Dec 8, 2008 #5
    But regardless of the growth rate 50000 will die per hour.

    P0ekt is the number of individuals if there were no deaths but there are 50 000 deaths per hour so I added the expression -50000*t ..

    Why should I divide by k?
  7. Dec 8, 2008 #6


    User Avatar
    Homework Helper

    You shouldn't have to modify the final expression for P(t) after solving the DE, only to find the constant of integration C. Dividing by k is not an extra step you must perform after re-arranging to find an expression for P(t); it's what you must do to get P(t).
  8. Dec 8, 2008 #7
    I tried to find the equation without solving the D.E... Is my solution still wrong then?
  9. Dec 8, 2008 #8


    User Avatar
    Science Advisor

    I presume that you now realize that your differential equation is dP/dt=kP-50 000.

    If it were dP/dt= kP, could you solve that?

    Suppose P were equal to some constant, A. What would A have to be to satisfy that equation.

    The solution to the entire equation is the sum of those.
  10. Jan 11, 2009 #9
    I would like to write here some examples of my models. I would be more happy if you could check my methods. Very easy way is to have a statistics. If you have a percentage of birth and deaths, you can write:
    extent of births [tex]\mu\geq0[/tex]
    extent of deaths [tex]\sigma\geq0[/tex]
    I can write an increment of P:
    [tex]\text{d}P(t)=\mu P(t)\text{d}t-\sigma P(t)\text{d}t[/tex]
    So it is a differential equation and its solution is function:

    There is also another possibility as wrote cango91. But I use incremental change. I presume there is an extent of deaths [tex]D[/tex] (number of died people per t) and there is an extent of births [tex]\mu\geq0[/tex]. I can write:
    [tex]\text{d}P(t)=\mu P(t)\text{d}t-D\text{d}t[/tex]
    I try to solve this:
    [tex]\begin{array}{rcl}\text{d}P(t)&=&\mu P(t)\text{d}t-D\text{d}t\\\text{d}&=&(\mu P(t)-D)\text{d}t\\\dfrac{\text{d}P(t)}{\mu P(t)-D}&=&\text{d}t\end{array}[/tex]
    The integral:
    [tex]\int\dfrac{\text{d}P(t)}{\mu P(t)-D}=\left|\begin{array}{rcl}\mu P(t)-D&=&\varphi\\\text{d}P(t)&=&\dfrac{1}{\mu}\,\text{d}\varphi\end{array}\right|=\frac{1}{\mu}\int\frac{\text{d}\varphi}{\varphi}=\frac{1}{\mu}\ln|\mu P(t)-D|[/tex]
    I use this into my equation:
    [tex]\begin{array}{rcl}\dfrac{1}{\mu}\ln|\mu P(t)-D|&=&t+\ln|C|\\\mu P(t)-D&=&C^{\mu}\text{e}^{\mu t}\\P(t)&=&\dfrac{1}{\mu}C^{\mu}\text{e}^{\mu t}+\dfrac{D}{\mu}\end{array}[/tex]
    The constant C for t = 0 is:
    [tex]C=\sqrt[\mu]{\mu P(0)-D}[/tex]
    So the final function is:
    [tex]P(t)=P(0)\text{e}^{\mu t}-\frac{D}{\mu}(\text{e}^{\mu t}-1)\,;\;\mu\geq0,D\geq0, t\geq0[/tex]
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook