[Date Index]
[Thread Index]
[Author Index]
Equilibrium points for 3 nonlinear ODE
*To*: mathgroup at smc.vnet.net
*Subject*: [mg70135] Equilibrium points for 3 nonlinear ODE
*From*: Virgil Stokes <vs at it.uu.se>
*Date*: Thu, 5 Oct 2006 03:32:28 -0400 (EDT)
I am trying to find the equilibrium points for this system of 3
nonlinear ODEs:
vars = {x[t], y[t], z[t]};
bi = r1*x[t]*(1 - x[t]/K1) - Î·*x[t]*y[
t] - a1*x[t]*z[t]/(1 + Î±*x[t] + Î²*y[t] + Î³*z[t]);
ra = r2*y[
t]*(1 - y[t]/K2) - a2*y[t]*z[t]/(1 + Î±*x[
t] + Î²*y[t] + Î³*z[t]) - Î¼r*y[t];
ca = -d*z[t] + ( b1*x[t]*z[t] +
b2*y[t]*z[t])/(1 + Î±*x[t] + Î²*y[t] + Î³*z[t]) - Î¼c*z[t];
eqns = {x'[t] == bi, y'[t] == ra, z'[t] == ca}
I then give,
equi = Solve[{bi == 0, ra == 0, ca == 0}, vars]
which gives:
\!\({x[t] -> 0, y[t] -> 0, z[t] -> 0}, {x[t] -> 0, z[t] -> 0, y[t] ->
\(K2 \((
r2 - Î¼r)\)\)\/r2}, \)
with another very, very, ... very long result. These first two solutions
are indeed correct; but, I do not get the solution
{x[t]->K1,y[t]->0,z[t]->0}, which on examination of the first equation
(bi) is clearly correct.
And I believe that these are the only solutions for equilibrium. How can
I "tweak" Mathematica to give me all three correct solutions?
--Thanks,
V. Stokes
Prev by Date:
**Re: Easy way to 're-order' the axes in Plot3D**
Next by Date:
**Re: Draw approximated normal distribution based on mean, median, percentiles?**
Previous by thread:
**Re: RE:Beginner--[Please Help]How to get coefficient list issued from a Series Command**
Next by thread:
**FrameTicks**
| |