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

Build logistic differential equation for given data

  1. Feb 17, 2008 #1
    http://img258.imageshack.us/img258/3518/01ck5.png [Broken]

    I would like to build a logistic differential equation for US oil production from the given data above.
    dP/dt = k(1-P/N)P

    I tried many different models but nothing came close to the given data.

    Any help and suggestion?
    Last edited by a moderator: May 3, 2017
  2. jcsd
  3. Feb 17, 2008 #2
    Why should it be logistic?
  4. Feb 18, 2008 #3
    This is a project from the textbook, there are many requirements (answering about 12 questions) and one of them is the equation must be logistic.
  5. Feb 18, 2008 #4
    I think you will have to model the resulting function with the data using a least squares fitting. This equation is:

    [tex]P(t)=\frac{N\cdot P_0 \cdot e^{kt}}{N+P_0\cdot \left(e^{kt}-1\right)}[/tex]

    It's very late over here, I will have a look at it tomorrow, it might be that it is non-linear and then I have to adapt my program to do this.
  6. Feb 18, 2008 #5
    If the above equation is correct, linear least squares while suffice. Just cross multiply:

    [tex]P(t)(N+P_0\cdot \left(e^{kt}-1\right))=N\cdot P_0 \cdot e^{kt}[/tex]

    Then rearrange:

    [tex]P(t)P_0\cdot \left(e^{kt}-1\right)=N (P_0 \cdot e^{kt} - P(t))[/tex]

    [tex]P_0[/tex] is the initial production level and use least squares to solve for N. I'm still not convinced oil production should be logistic. The only reason it should be logistic is because the production is driven by the demand which is driven by the population.
    Last edited: Feb 18, 2008
  7. Feb 19, 2008 #6
    You need some code to do this least squares fitting in case one wants to calculate N, P0 and k at the same time. When this is done you will end up with N=14.9395; P0=0.23349 and k=0.102758 for the US and N=116.435; P0=0.04219 and k=0.12401 for the world production. This in case the time is considered the center of the intervals as is usually done in statistics. I used my own code for this, however I assume that this must be possible in some other mathematical package available. I can't help you in this though. One final remark is the following, I'm not giving any advice on whether this is a good curve for describing the data, this is something you need to do yourself. Other curves or equations might be more suitable. If you look at the one for the US, it is strange data, it's going down. Interpret this very carefully.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook