Does the following method suffice for proving that a system of three first-order, non-linear ODEs always exhibits damped oscillations (as opposed to undamped), using the system's characteristic polynomial? If not, why?
Method
Calculate the system's non-trivial equilibrium (assuming there is one).
Calculate the characteristic polynomial corresponding to that equilibrium. This polynomial is third order and therefore has three eigenvalues
Show that the characteristic polynomial's discriminate is always negative (this shows that the equilibrium exhibits oscillations, right?)
Show that the real part of the leading eigenvalue is always negative (this shows that said oscillations decay in amplitude toward zero, right?)
I am not exactly sure about the method you are describing, but here is how I would do it (for a system $\dot{x} = f(x)$):
All of this follows from the general solution of matrix differential equations, which describe the linear approximation of the dynamics around $x_0$. Of course, the dynamics may exhibit damped oscillations far from the fixed point that you cannot detect with this method.
Some comments on your method as you describe it:
If you have three second-order differential equations, you get six first-order ones and thus your characteristic polynomial would have degree six – not three.
Since the polynomial has degree six, it does not have a clear discriminant.