How i can solve limit this ODE question .

4 vues (au cours des 30 derniers jours)
Jumanah Al omran
Jumanah Al omran le 1 Mai 2022
Commenté : Sam Chak le 8 Mai 2022
The population P(x) of a species satisfies the logistic differential equation.the initial population P(0)= 3000 and t is the in years. How I can find the lim x -> infinity P(t) ? this is the equation
  4 commentaires
Jan
Jan le 7 Mai 2022
Modifié(e) : Jan le 7 Mai 2022
@Jumanah Al omran: You are funny. The unit does not matter. For ->Inf you get the same result for years or seconds, and even for bananas. :-)
You can determine the equation for x(t) easily. Then it is easy to set x to Inf are determine the result.
A deeper look on the equation for dx/dt helps also. x is growing, -x^2 is shrinking even faster. What dows a curve with a rapidly shrinking derivative do for large values?
Sam Chak
Sam Chak le 7 Mai 2022
Modifié(e) : Sam Chak le 7 Mai 2022
Based on your description, P(x) is the population as a function of x (the years...).
And dx/dt?
By the way, I plotted this to help you visualize the rate dx/dt. You see, x grows initially but the growth starts declining after x > 500K. It is still growing but getting slower and slower until dx/dt comes to a halt when it reaches 0.
x = linspace(0, 1e6, 100001);
dx = x/100 - x.^2/1e8;
plot(x, dx)

Connectez-vous pour commenter.

Réponses (3)

Bjorn Gustavsson
Bjorn Gustavsson le 7 Mai 2022
This is a Riccati equation. Convert it to a standard linear second-order differential equation and solve it that way, check what transformation to use from: Riccati-equation
Solve the ODE by separation of variables.
At what level will the source and loss-terms on the RHS be equal?
HTH

Torsten
Torsten le 7 Mai 2022
I think you mean
lim (t->oo) P(t)
instead of
lim (x->oo) P(t)
don't you ?
Seems the limit is not 0, but 1e6:

Sam Chak
Sam Chak le 8 Mai 2022
Looks like the population saturates at 1 million.
t = linspace(0, 1500);
P = (3000000*exp(t/100))./(997 + 3*exp(t/100));
plot(t, P)
  2 commentaires
Bjorn Gustavsson
Bjorn Gustavsson le 8 Mai 2022
Yes. This is obvoius if we rewrite the ODE as:
dxdt = x/100.*(1-x/1e6);
First factor will be positive for all positive x, second factor goes to zero as x->1e6.
Sam Chak
Sam Chak le 8 Mai 2022
Thanks @Bjorn Gustavsson. I forgot to factor the quadratic function.

Connectez-vous pour commenter.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by