- #1

- 6

- 0

[itex]\frac{dx}{dt} = 1 - 2x + x^2y [/itex]

[itex]\frac{dy}{dt} = x-x^2y[/itex]

It asks me to find the fixed point(s), and determine their stability, also to draw the phase plane.

So to find the fixed points, I set both equations equal to zero:

[itex]1 - 2x + x^2y = 0 \Rightarrow x^2y = 2x-1[/itex]

[itex]x-x^2y = 0 \Rightarrow x^2y = x[/itex]

Setting them equal to each other:

[itex]2x-1 = x \Rightarrow x=1[/itex] and from this we can deduce [itex]y=1[/itex]. So the fixed point is (1,1).

Now this is where I'm not really sure on what to do and get confused by my lecture notes. I believe I have to use some Jacobian to work out the stability yet I have no clue why i'm doing this (I know eigenvalues/eigenvectors are involved somehow). Setting [tex]f(x,y)= 1 - 2x + x^2y[/tex] and [tex]g(x,y)= x-x^2y[/tex] I worked out the Jacobian at (1,1) to be:

\begin{equation}

\left|

\begin{array}{cc}

0 & 1 \\

-1 & -1 \end{array}

\right|

\end{equation}

(Sorry about the lack of detail about this part but I'm new to latex and struggled with doing matrices, the jacobian is f

_{x}and f

_{y}on the first row and g

_{x}and g

_{y}on the second row.

Following lecture notes, I work out the 'trace' (T) of this matrix which is simply (a

_{11}+a

_{22}) which is -1 in this case, and also work out the determinant which is 1.

Since [itex]T^2 - 4D = -3 < 0[/itex], we have complex eigenvalues and I believe this means the phase plane will be a spiral, and since T < 0 this means the general trajactory will be moving towards the fixed point (stable). But do I need to work out the eigenvectors to see exactly what this spiral will look like? I'm not sure if I've worked out enough to complete the question. Also any help on actually explaining why I'm working out the Jacobian stuff would be great.

Thanks