I have Lotka Volterra Equations where the prey equation is modified. The Equations are
x'[t]=a*x[t]*(1-x[t]/k)-b*x[t]*y[t]
y'[t]=-c*y[t]+d*x[t]*y[t]
To plot the limit cycles, i used NDSolve
to solve the coupled nonlinear differential equations for x[t] and y[t]. I then used ParametricPlot
to plot for x[t]
and y[t]
with respect to t
, but I'm not getting appropriate limit cycles.
I did a steady state analysis of the equations and got certain conditions on the value of a
, b
, c
, d
and k
. The equation has a non-trivial steady state, at $(x,y)=(c/d,(a/b)-(ac/bdk))$. The commands are:
Equations = {x'[t] == ax[t]*a*x[t]*(1 - x[t]/k) - bx[t]*y[t]b*x[t]*y[t], x[0] == 2, y[0] == 1, y'[t] == -cy[t]c*y[t] + dx[t]*y[t]d*x[t]*y[t]};
s = NDSolve[{Equations}NDSolve[Equations, {x, y}, {t, 0, time}];
ParametricPlot[Evaluate[{x[t], y[t]} /. s], {t, 0, 30}]
For the value of a = 2; b = 1; c = 2.5; d = 1.2; k = 2.8
, I'm supposed to get a unstable spiral. But instead, the trajectory approaches the steady state of $(x,y) = (2.08333 , 0.511905)$.
Any help will be greatly appreciated.