I've been hearing about using the sign of the real part of eigenvalues of the multivariate derivative of the change vector field of an ODE at an equilibrium point to determine whether that equilibrium is stable or unstable. Unfortunately, the source is not very good at things like getting the statement of theorems correct, so I was wondering:


Is there a statement and proof of such a theorem I could access online?


I'm particularly interested in if anything can be said about the non-diagonalizable case where the eigenvalues have zero real part.



Linearization at equilibrium points
I got curious about "the non-diagonalizable case" with purely imaginary eigenvalues
and came up with the system
$$ y' = A y $$ where
$$ 
A \; = \; \left( \begin{array}{rrrr}
 0 & -1 & 1 & 0 \\\
 1 & 0 & 0 & 1 \\\
 0 & 0 & 0 & -1 \\\
 0 & 0 & 1 & 0 
\end{array} 
 \right) .
 $$
Taking the matrix of eigenvectors and "generalized eigenvectors"
(see <http://en.wikipedia.org/wiki/Generalized_eigenvector> )
$$ 
B \; = \; \left( \begin{array}{rrrr}
 1 & i & 0 & 0 \\\
 i & 1 & 0 & 0 \\\
 0 & 0 & 1 & i \\\
 0 & 0 & i & 1 
\end{array} 
 \right) .
 $$
we get a permuted Jordan form in the shape I prefer for this case,
$$ 
B^{-1} A B \; = \; \left( \begin{array}{rrrr}
 -i & 0 & 1 & 0 \\\
 0 & i & 0 & 1 \\\
 0 & 0 & -i & 0 \\\
 0 & 0 & 0 & i 
\end{array} 
 \right) .
 $$
This is something I made up years ago and forgot, with an original real matrix and a repeated characteristic value $\alpha$ with non-diagonal Jordan block, we get the same block for the complex conjugate $\bar{\alpha}.$ For both the original real matrix $A$ and 
$B^{-1} A B $ we can rearrange everything into convenient 2 by 2 blocks, in particular the off-diagonal 1's become little 2 by 2 identity matrices. Try it for a six by six example where there is a 3 by 3 block for $\alpha$ with two off-diagonal ones, then a 3 by 3 block for $\bar{\alpha}$ with two off-diagonal ones. Permute the diagonal elements to
$ \alpha, \; \bar{\alpha}, \; \alpha, \; \bar{\alpha}, \; \alpha, \; \bar{\alpha} $ and see what happens. By the way, these "forms" probably have standard names.


After a bunch of work I found that the "fundamental matrix" for the system is
$$ 
e^{A t} \; = \; \left( \begin{array}{rrrr}
 \cos t & - \sin t & t \cos t & - t \sin t \\\
 \sin t & \cos t & t \sin t & t \cos t \\\
 0 & 0 & \cos t & - \sin t \\\
 0 & 0 & \sin t & \cos t 
\end{array} 
 \right) .
 $$
which means that any solution of the system is $ y = e^{A t} y\_0.$
So, if the third and fourth entries in $y\_0$ are $0,$ the orbit is a circle. If not, the orbit leaves the origin as $t$ increases. The trick from my first answer, making a not quite linear system, can force the periodic orbits to switch to attracting.


