Population Dynamics how to solve a particular ODE

jerro
Messages
7
Reaction score
0

Homework Statement



I have a population problem where:

\frac{dy}{dt} = ay - by^{2}-\frac{c*y^{3}}{d+y^{3}}

I need to find an expression for y(t). I'm not looking for the answer, just some advice/ helpful hints.

Thank you.

Homework Equations





The Attempt at a Solution



I know that the integrating factor method is out of the question, as is separation of variables. Bernoulli's equation will also not work. All of these require forms that are different than what is written above.
 
Physics news on Phys.org
jerro said:

Homework Statement



I have a population problem where:

\frac{dy}{dt} = ay - by^{2}-\frac{c*y^{3}}{d+y^{3}}

I need to find an expression for y(t). I'm not looking for the answer, just some advice/ helpful hints.

Thank you.

Homework Equations





The Attempt at a Solution



I know that the integrating factor method is out of the question, as is separation of variables. Bernoulli's equation will also not work. All of these require forms that are different than what is written above.

Getting t in terms of y is not too hard, but getting y in terms of t is horrible. If you write
F(y) \equiv \int\frac{dy}{f(y)} = \int dt,\\<br /> f(y) = a y - b y^2 - \frac{c y^3}{d + y^3},<br />
the y-integral is doable. Using 'r' instead of 'd' (because 'd' is a reserved symbol), Maple gets the y-integral as

F(y) = 1/2*a*x^2-1/3*b*x^3-c*x+1/3*c*r^(1/3)*ln(x+r^(1/3))
-1/6*c*r^(1/3)*ln(x^2-x*r^(1/3)+r^(2/3))
+1/3*c*r^(1/3)*3^(1/2)*arctan(1/3*3^(1/2)*(2/r^(1/3)*x-1))

where Maple writes u/v*w to mean (u/v)*w. So, F(y) = t+C is a "solution", but inverting F to get y(t) is probably only doable numerically.
 
Ah, I see. Separation of variables is completely doable! I don't know what I was thinking. Thank you so much.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top