How Do You Solve the Verhulst Equation for Logistic Population Growth?

  • Thread starter Thread starter gpax42
  • Start date Start date
Join the discussion
Registration is free. Start your own thread to ask a follow-up.
2 replies · 7K views
gpax42
Messages
25
Reaction score
0
Hi all, it's been a over a year since I took my differential equations and linear algebra course and I'm currently enrolled in a class that assigned this problem as a sort of refresher on analytically solving differential equations. I can't seem to remember the proper approach to going about all this and I'm no longer in possession of my textbook from the previous course. Any help would be greatly appreciated :redface:

Homework Statement

Using methods for solving differential equations exactly, solve the Verhulst equation for logistic population growth

[tex]\frac{dP}{dt}[/tex]= rP(1 - [tex]P/K[/tex])

where r is the growth rate and K is the carrying capacity.


The attempt at a solution

I began by altering the equation so as to create a more solvable form...

Dividing through by K on both sides gives me

[tex]\frac{d}{dt}[/tex][tex]\frac{P}{K}[/tex] = r[tex]\frac{P}{K}[/tex](1 - [tex]P/K[/tex])

I then set x = [tex]P/K[/tex]

[tex]\frac{dx}{dt}[/tex] = rx(1-x)

at this point I know I need to integrate both sides

dx = rx(1-x)dt
following integration...
x(t) = [rx(1-x)]t + C

is this even remotely correct?

Thanks again for any help you can lend me

-gpax42
 
Physics news on Phys.org
gpax42 said:
Hi all, it's been a over a year since I took my differential equations and linear algebra course and I'm currently enrolled in a class that assigned this problem as a sort of refresher on analytically solving differential equations. I can't seem to remember the proper approach to going about all this and I'm no longer in possession of my textbook from the previous course. Any help would be greatly appreciated :redface:

Homework Statement

Using methods for solving differential equations exactly, solve the Verhulst equation for logistic population growth

[tex]\frac{dP}{dt}[/tex]= rP(1 - [tex]P/K[/tex])

where r is the growth rate and K is the carrying capacity.


The attempt at a solution

I began by altering the equation so as to create a more solvable form...

Dividing through by K on both sides gives me

[tex]\frac{d}{dt}[/tex][tex]\frac{P}{K}[/tex] = r[tex]\frac{P}{K}[/tex](1 - [tex]P/K[/tex])

I then set x = [tex]P/K[/tex]

[tex]\frac{dx}{dt}[/tex] = rx(1-x)

at this point I know I need to integrate both sides

dx = rx(1-x)dt
following integration...
x(t) = [rx(1-x)]t + C

is this even remotely correct?

-gpax42

No, it isn't. First, essentially changing the dependent variable from P to x accomplishes nothing except renaming it. Your x is an unknown function of t and its integral isn't xt.

You want to review separation of variables. Separate the P and t:

[tex]\frac {dP}{P(1-P)}= \frac r K dt[/tex]

Then integrate both sides. You will want to use partial fractions on the left side.
 
ahhh I see... when I separated variables, however, I found...

[tex]\int[/tex]dP/[P(1-P/k)] = [tex]\int[/tex]rdt

I then manipulated the left side to read

K/[P(K-P)] ... after partial fractions... = 1/P + 1/(K-P)

int[1/P + 1/(K-P) )dP] = int[rdt]

following through with integration leads to

ln(abs(P/(K-P))) = rt + C

e^(rt + C) = abs[P/(K-P)]

e^(-rt-C) = abs[(K-P)/P]

+/- e^(-rt) * e^(-C) = K/P - 1

if you let M = +/- e^(-C)
then
M*e^(-rt) = K/P - 1

solving further...


P(t) = K / (1+Me^(-rt))

at t = 0 ; P = P0

P0 = K / (1+M)

M =( K - P0 ) / P0