Differential Equations - Population Dynamics

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
4 replies · 4K views
hawks32
Messages
3
Reaction score
0
1. The problem statement
The DE governing a fish pop. P(t) with harvesting proportional to the population is given by:
P'(t)=(b-kP)P-hP
where b>0 is birthrate, kP is deathrate, where k>0, and h is the harvesting rate. Model assumes that the death rate per individual is proportional to the pop. size. An equilibrium point for the DE is a value of P so that P'(t)=0.

Find general solution of the DE, when..
a) h>b
b) h=b
c) h<b

The Attempt at a Solution


I'm having problems figuring out how to set up parts a) and c). I'm horrible at DE, so if anyone could help point me in the right direction, it would be greatly appreciated.
 
Physics news on Phys.org
P'(t)=(b-kP)P-hP
P'(t)=bP-kP2-hP
P'(t)= (b-h)P-kP2

P'(t)= ((b-h)-kP)P

P'(t)= dP/dt

so put it in the form

f(P) dP= f(t) dt

then integrate both sides.
 
okay, i worked the integral of dp/dt = ((b-h)-kP)P out as...



1/(b-h) * ln(p/b-h-kP) + C = t

is that correct?
 
Umm you have a [tex]P^2[/tex] in there. You should try partial fractions.
 
I did use partial fractions.

1/((b-h)-kP)P dp Let a = b-h

integral of 1/(a-kP)P = integral of A/a-kP + B/P

Solved for A & B, A = k/a, B = 1/a

So integral (k/a)/(a-kP) + (1/a)/p

end up with -(1/a)ln(a-kP) + (1/a)lnP
==> 1/a ln(P/(a-kP)) + C

sub back a = b-h

1/(b-h) * ln(p/(b-h-kP)) + C

Did I do something wrong?