Morse function is used as follows. We have set of nonlinear first order ODE’s. We’d like to determine if the nonlinear system is stable or not at its critical (or equilibrium) points. This is done by first finding the critical points and then linearizing the system around each one of these points, since it is much easier to check if the system is stable or not when it is linear by finding the eigenvalues of the Jacobian matrix.
But how do we know that behavior of the linear system around the critical points is the same as the nonlinear system? In other words, if the linearized system is stable at a critical point, does this mean also the nonlinear is stable? By finding Morse function, called \(F\left ( x_{1},x_{2},\cdots \right ) \) below, and then by finding the determinant of its Hessian evaluated at each critical point, i.e. \(\det \left ( \nabla ^{2}F\right ) \), then if the determinant is not zero, then the critical point is called non-degenerate, which means the nonlinear system behaves the same as the linear system around the critical point. So if the linear system says the critical point is stable, we know the nonlinear system is also stable. But if \(\det \left ( \nabla ^{2}F\right ) =0\) then the point is called degenerate. In this case, the nonlinear system is not guaranteed to behave the same was as the linearized system.
Morse function is basically the first integral at the critical point, and it always come out having a quadratic form (in the state variables of the system). If all terms in the quadratic form have positive sign, then the system is stable at the critical point. If one term in the quadratic form have negative sign, then the system is unstable there. But this all makes sense only if the point is non-degenerate, so we should always check the point is non-degenerate before trying to find Morse function.
To find Morse function, we need to first find what is called the first integral, called \(\tilde {F}\left ( x_{1},x_{2},\cdots \right ) \) below, of the set of given ODE’s. The first integral \(\tilde {F}\left ( x_{1},x_{2},\cdots \right ) \) normally represents the energy of the system. For example, the Hamiltonian is a first integral. But in general, it is any function \(\tilde {F}\left ( x_{1},x_{2},\cdots \right ) \) which is constant along each solution path. To check that a given \(\tilde {F}\left ( x_{1},x_{2},\cdots \right ) \) is a first integral of the system, we can do this test\[ \frac {\partial \tilde {F}}{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \tilde {F}}{\partial _{x_{2}}}\dot {x}_{2}+\cdots +\frac {\partial \tilde {F}}{\partial _{x_{n}}}\dot {x}_{n}=0 \] If the above test fail to give zero then \(\tilde {F}\left ( x_{1},x_{2},\cdots \right ) \) is not a first integral of the set of ODE’s, and most likely we made mistake in finding it in first place. Also, not all set of ODE’s have a first integral and it can be hard sometimes to find the first integral analytically since it involves solving nonlinear set of ODE’s if the system is nonlinear.
\[ \ddot {x}+\sin x=0 \] step 1 convert to state space. Let \(x_{1}=x,x_{2}=\dot {x}\). Taking derivative w.r.t. gives \(\dot {x}_{1}=x_{2},\dot {x}_{2}=\ddot {x}=-\sin x_{1}\). Hence\[\begin {pmatrix} \dot {x}_{1}\\ \dot {x}_{2}\end {pmatrix} =\begin {pmatrix} x_{2}\\ -\sin x_{1}\end {pmatrix} =\begin {pmatrix} f_{1}\\ f_{2}\end {pmatrix} \] step 2 Find first integral\[ \frac {dx_{2}}{dx_{1}}=\frac {f_{2}}{f_{1}}=\frac {-\sin x_{1}}{x_{2}}\] Integrating\begin {align*} \int x_{2}dx_{2} & =-\int \sin x_{1}dx_{1}\\ \frac {1}{2}x_{2}^{2}-\cos x_{1} & =E \end {align*}
Where \(E\) is the constant of integration. Let \(\bar {F}\left ( x_{1},x_{2}\right ) \equiv \frac {1}{2}x_{2}^{2}-\cos x_{1}\).
step 2.1 Verify that the above is indeed first integral by checking that \(\frac {\partial \tilde {F}}{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \tilde {F}}{\partial _{x_{2}}}\dot {x}_{2}+\cdots +\frac {\partial \tilde {F}}{\partial _{x_{n}}}\dot {x}_{n}=0\). Hence \[ \frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2}=\sin x_{1}\dot {x}_{1}+x_{2}\dot {x}_{2}\] But \(\dot {x}_{2}=\ddot {x},x_{1}=x,\dot {x}_{1}=x_{2}\). Hence the above becomes\begin {align*} \frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2} & =\sin xx_{2}+x_{2}\ddot {x}\\ & =x_{2}\left ( \sin x+\ddot {x}\right ) \end {align*}
But \(\sin x+\ddot {x}=0\) since this is the ODE given. Therefore \(\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2}=0\) as expected.
step 3 Find critical points by solving\begin {align*} \begin {pmatrix} f_{1}\\ f_{2}\end {pmatrix} & =\begin {pmatrix} 0\\ 0 \end {pmatrix} \\\begin {pmatrix} x_{2}\\ -\sin x_{1}\end {pmatrix} & =\begin {pmatrix} 0\\ 0 \end {pmatrix} \end {align*}
Hence \(x_{2}=0\) and \(x_{1}=\left \{ 0,\pi ,-\pi \right \} \). Therefore the critical points are \(\left ( 0,0\right ) ,\left ( \pi ,0\right ) ,\left ( -\pi ,0\right ) \)
step 4 Find Hessian of \(\bar {F}\left ( x_{1},x_{2}\right ) =\frac {1}{2}x_{2}^{2}-\cos x_{1}\) at each critical point and determine if the determinant of the Hessian at each critical point is non-degenerate. \begin {align*} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) & =\begin {pmatrix} \frac {\partial ^{2}F}{\partial x_{1}\partial x_{1}} & \frac {\partial ^{2}F}{\partial x_{1}\partial x_{2}}\\ \frac {\partial ^{2}F}{\partial x_{2}\partial x_{1}} & \frac {\partial ^{2}F}{\partial x_{2}\partial x_{2}}\end {pmatrix} \\ & =\begin {pmatrix} \frac {\partial }{\partial x_{1}}\left ( \sin x_{1}\right ) & \frac {\partial F}{\partial x_{1}}x_{2}\\ \frac {\partial F}{\partial x_{2}}\left ( \sin x_{1}\right ) & \frac {\partial F}{\partial x_{2}}x_{2}\end {pmatrix} \\ & =\begin {pmatrix} \cos x_{1} & 0\\ 0 & 1 \end {pmatrix} \end {align*}
Hence \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=0,x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} 1 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(1\neq 0\). Hence non-degenerate.
\(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=\pi ,x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} \cos \pi & 0\\ 0 & 1 \end {pmatrix} =\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(-1\neq 0\). Hence non-degenerate.
\(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=-\pi ,x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} \cos \left ( -\pi \right ) & 0\\ 0 & 1 \end {pmatrix} =\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(-1\neq 0\). Hence non-degenerate.
Since \(\det \left ( \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \right ) \neq 0\) at each critical point, then Morse function exist.
step 5 Find Morse function \(F\left ( x_{1},x_{2}\right ) \) associated with \(\bar {F}\left ( x_{1},x_{2}\right ) \) at each critical point.
For \(\left ( 0,0\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( 0,0\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes\begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\\ & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=0\\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=0\\x_{2}=0}}}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\cos x_{1}\right ) _{\substack {x_{1}=0\\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} \cos x_{1} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=0\\x_{2}=0}}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =-1+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} 1 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\left ( x_{1}^{2}+x_{2}^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =x_{1}^{2}+x_{2}^{2}\\ & =E \end {align*}
Where \(1\) and \(\frac {1}{2}\) are combined into the constant \(E\). Since there is no negative term in the quadratic form above, then the index of the critical point is \(0\) which means this is stable critical point.
For \(\left ( \pi ,0\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( \pi ,0\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes\begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=\pi \\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1}-\pi & x_{2}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=\pi \\x_{2}=0}}}\begin {pmatrix} x_{1}-\pi \\ x_{2}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\cos x_{1}\right ) _{\substack {x_{1}=\pi \\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1}-\pi & x_{2}\end {pmatrix}\begin {pmatrix} \cos x_{1} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=\pi \\x_{2}=0}}\begin {pmatrix} x_{1}-\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\begin {pmatrix} x_{1}-\pi & x_{2}\end {pmatrix}\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}-\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\begin {pmatrix} -x_{1}+\pi & x_{2}\end {pmatrix}\begin {pmatrix} x_{1}-\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\left ( \left ( -x_{1}+\pi \right ) \left ( x_{1}-\pi \right ) +x_{2}^{2}\right ) \\ & =1+\frac {1}{2}\left ( -\left ( x_{1}-\pi \right ) ^{2}+x_{2}^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =-\left ( x_{1}-\pi \right ) ^{2}+x_{2}^{2}\\ & =E \end {align*}
Since there is a negative term in the quadratic form above, then the index of this critical point is \(-1\) which means this is unstable critical point. (index must be zero for stable critical point)
For \(\left ( -\pi ,0\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( -\pi ,0\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes\begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=-\pi \\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1}+\pi & x_{2}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=-\pi \\x_{2}=0}}}\begin {pmatrix} x_{1}+\pi \\ x_{2}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\cos x_{1}\right ) _{\substack {x_{1}=-\pi \\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1}+\pi & x_{2}\end {pmatrix}\begin {pmatrix} \cos x_{1} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=-\pi \\x_{2}=0}}\begin {pmatrix} x_{1}+\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\begin {pmatrix} x_{1}+\pi & x_{2}\end {pmatrix}\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}+\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\begin {pmatrix} -\left ( x_{1}+\pi \right ) & x_{2}\end {pmatrix}\begin {pmatrix} x_{1}+\pi \\ x_{2}\end {pmatrix} \\ & =1+\frac {1}{2}\left ( -\left ( x_{1}+\pi \right ) ^{2}+x_{2}^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =-\left ( x_{1}+\pi \right ) ^{2}+x_{2}^{2}\\ & =E \end {align*}
Since there is a negative term in the quadratic form above, then the index of this critical point is \(-1\) which means this is unstable critical point. (index must be zero for stable critical point)
\begin {align*} \dot {x} & =y\\ \dot {y} & =x-2x^{3} \end {align*}
step 1 \[\begin {pmatrix} \dot {x}_{1}\\ \dot {x}_{2}\end {pmatrix} =\begin {pmatrix} x_{2}\\ x_{1}-2x_{1}^{3}\end {pmatrix} =\begin {pmatrix} f_{1}\\ f_{2}\end {pmatrix} \] step 2 Find first integral\[ \frac {dx_{2}}{dx_{1}}=\frac {f_{2}}{f_{1}}=\frac {x_{1}-2x_{1}^{3}}{x_{2}}\] Integrating\begin {align*} \int x_{2}dx_{2} & =\int x_{1}-2x_{1}^{3}dx_{1}\\ \frac {1}{2}x_{2}^{2} & =\frac {1}{2}x_{1}^{2}-\frac {2}{4}x_{1}^{4}+E\\ \frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4} & =E \end {align*}
Where \(E\) is the constant of integration. Let \(\bar {F}\left ( x_{1},x_{2}\right ) \equiv \frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4}\). The above is also the equation of orbit.
step 2.1 Verify that the above is indeed first integral by checking that \(\frac {\partial \tilde {F}}{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \tilde {F}}{\partial _{x_{2}}}\dot {x}_{2}+\cdots +\frac {\partial \tilde {F}}{\partial _{x_{n}}}\dot {x}_{n}=0\). Hence \[ \frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2}=\left ( -x_{1}+2x_{1}^{3}\right ) \dot {x}_{1}+x_{2}\dot {x}_{2}\] But \(x_{1}=x,\dot {x}_{1}=y,x_{2}=y,\dot {x}_{2}=\dot {y}\) Hence the above becomes\begin {align*} \frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2} & =\left ( -x_{1}+2x_{1}^{3}\right ) y+y\dot {y}\\ & =y\left ( \left ( -x_{1}+2x_{1}^{3}\right ) +\dot {y}\right ) \end {align*}
But \(\dot {y}=x-2x^{3}\) as given. The above becomes \begin {align*} \frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{1}}}\dot {x}_{1}+\frac {\partial \bar {F}\left ( x_{1},x_{2}\right ) }{\partial _{x_{2}}}\dot {x}_{2} & =y\left ( -x_{1}+2x_{1}^{3}+x-2x^{3}\right ) \\ & =0 \end {align*}
as expected for first integral.
step 3 Find critical points by solving\begin {align*} \begin {pmatrix} f_{1}\\ f_{2}\end {pmatrix} & =\begin {pmatrix} 0\\ 0 \end {pmatrix} \\\begin {pmatrix} x_{2}\\ x_{1}-2x_{1}^{3}\end {pmatrix} & =\begin {pmatrix} 0\\ 0 \end {pmatrix} \end {align*}
Therefore the critical points are \(\left ( x_{1},x_{2}\right ) =\left \{ \left ( 0,0\right ) ,\left ( \frac {1}{\sqrt {2}},0\right ) ,\left ( -\frac {1}{\sqrt {2}},0\right ) \right \} \)
step 4 Find Hessian of \(\bar {F}\left ( x_{1},x_{2}\right ) =\frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4}\) at each critical point and determine if the determinant of the Hessian at each critical point is non-degenerate. \begin {align*} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) & =\begin {pmatrix} \frac {\partial ^{2}F}{\partial x_{1}\partial x_{1}} & \frac {\partial ^{2}F}{\partial x_{1}\partial x_{2}}\\ \frac {\partial ^{2}F}{\partial x_{2}\partial x_{1}} & \frac {\partial ^{2}F}{\partial x_{2}\partial x_{2}}\end {pmatrix} \\ & =\begin {pmatrix} \frac {\partial }{\partial x_{1}}\left ( -x_{1}+2x_{1}^{3}\right ) & \frac {\partial F}{\partial x_{1}}x_{2}\\ \frac {\partial F}{\partial x_{2}}\left ( -x_{1}+2x_{1}^{3}\right ) & \frac {\partial F}{\partial x_{2}}x_{2}\end {pmatrix} \\ & =\begin {pmatrix} -1+6x_{1}^{2} & 0\\ 0 & 1 \end {pmatrix} \end {align*}
Hence \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=0,x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(-1\neq 0\). Hence non-degenerate.
\(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=\frac {1}{\sqrt {2}},x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} -1+3 & 0\\ 0 & 1 \end {pmatrix} =\begin {pmatrix} 2 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(2\neq 0\). Hence non-degenerate.
\(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \) evaluated at \(x_{1}=-\frac {1}{\sqrt {2}},x_{2}=0\) is \(\nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) =\begin {pmatrix} -1+3 & 0\\ 0 & 1 \end {pmatrix} =\begin {pmatrix} 2 & 0\\ 0 & 1 \end {pmatrix} \). The determinant of this is \(2\neq 0\). Hence non-degenerate.
Since \(\det \left ( \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) \right ) \neq 0\) at each critical point, then Morse function exist.
step 5 Find Morse function \(F\left ( x_{1},x_{2}\right ) \) associated with \(\bar {F}\left ( x_{1},x_{2}\right ) \) at each critical point.
For \(\left ( 0,0\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( 0,0\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes \begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=0\\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=0\\x_{2}=0}}}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4}\right ) _{\substack {x_{1}=0\\x_{2}=0}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} -1+6x_{1}^{2} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=0\\x_{2}=0}}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =0+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} -1 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =\frac {1}{2}\begin {pmatrix} -x_{1} & x_{2}\end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\end {pmatrix} \\ & =\frac {1}{2}\left ( -x_{1}^{2}+x_{2}^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =-x_{1}^{2}+x_{2}^{2}\\ & =E \end {align*}
Since there is a negative term in the quadratic form above, then the index of this critical point is \(-1\) which means this is unstable critical point. It is a saddle. (index must be zero for stable critical point)
For \(\left ( 0,\frac {1}{\sqrt {2}}\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( 0,\frac {1}{\sqrt {2}}\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes\begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=0\\x_{2}=\frac {1}{\sqrt {2}}}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=0\\x_{2}=\frac {1}{\sqrt {2}}}}}\begin {pmatrix} x_{1}\\ x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4}\right ) _{\substack {x_{1}=0\\x_{2}=\frac {1}{\sqrt {2}}}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} -1+6x_{1}^{2} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=0\\x_{2}=\frac {1}{\sqrt {2}}}}\begin {pmatrix} x_{1}\\ x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} 2 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\begin {pmatrix} 2x_{1} & x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}-\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\left ( 2x_{1}^{2}+\left ( x_{2}-\frac {1}{\sqrt {2}}\right ) ^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =2x_{1}^{2}+\left ( x_{2}-\frac {1}{\sqrt {2}}\right ) ^{2}\\ & =E \end {align*}
Since there is no negative term in the quadratic form above, then the index of the critical point is \(0\) which means this is stable critical point. (Center node)
For \(\left ( 0,\frac {-1}{\sqrt {2}}\right ) \): Taylor expansion of \(F\left ( x_{1},x_{2}\right ) \) around \(a=\left ( 0,\frac {-1}{\sqrt {2}}\right ) \) gives\[ F\left ( x_{1},x_{2}\right ) =\bar {F}\left ( a\right ) +\left ( x-a\right ) \nabla F\left ( a\right ) +\frac {1}{2}\left ( x-a\right ) \nabla ^{2}F\left ( a\right ) \left ( x-a\right ) ^{T}\] But \(\nabla F\left ( a\right ) =0\) since \(a\) is critical point. The above becomes \begin {align*} F\left ( x_{1},x_{2}\right ) & =\bar {F}\left ( x_{1},x_{2}\right ) _{\substack {x_{1}=0\\x_{2}=\frac {-1}{\sqrt {2}}}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix} \nabla ^{2}\bar {F}\left ( x_{1},x_{2}\right ) _{_{\substack {x_{1}=0\\x_{2}=\frac {-1}{\sqrt {2}}}}}\begin {pmatrix} x_{1}\\ x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\left ( \frac {1}{2}x_{2}^{2}-\frac {1}{2}x_{1}^{2}+\frac {1}{2}x_{1}^{4}\right ) _{\substack {x_{1}=0\\x_{2}=\frac {-1}{\sqrt {2}}}}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} -1+6x_{1}^{2} & 0\\ 0 & 1 \end {pmatrix} _{\substack {x_{1}=0\\x_{2}=\frac {-1}{\sqrt {2}}}}\begin {pmatrix} x_{1}\\ x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\begin {pmatrix} x_{1} & x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} 2 & 0\\ 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\begin {pmatrix} 2x_{1} & x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}+\frac {1}{\sqrt {2}}\end {pmatrix} \\ & =\frac {1}{4}+\frac {1}{2}\left ( 2x_{1}^{2}+\left ( x_{2}+\frac {1}{\sqrt {2}}\right ) ^{2}\right ) \end {align*}
Since \(F\left ( x_{1},x_{2}\right ) \) is constant, say \(E\), then the above can be written as\begin {align*} F\left ( x_{1},x_{2}\right ) & =2x_{1}^{2}+\left ( x_{2}+\frac {1}{\sqrt {2}}\right ) ^{2}\\ & =E \end {align*}
Since there is no negative term in the quadratic form above, then the index of the critical point is \(0\) which means this is stable critical point. (center node)
\begin {align} \dot {x}_{1} & =x_{1}-x_{1}x_{2}-x_{2}^{3}+x_{3}\left ( x_{1}^{2}+x_{2}^{2}-1-x_{1}+x_{1}x_{2}+x_{2}^{3}\right ) \tag {A}\\ \dot {x}_{2} & =x_{1}-x_{3}\left ( x_{1}-x_{2}+2x_{1}x_{2}\right ) \nonumber \\ \dot {x}_{3} & =\left ( x_{3}-1\right ) \left ( x_{3}+2x_{3}x_{2}^{2}+x_{3}^{3}\right ) \nonumber \end {align}
step 1 \[\begin {pmatrix} \dot {x}_{1}\\ \dot {x}_{2}\\ \dot {x}_{3}\end {pmatrix} =\begin {pmatrix} x_{1}-x_{1}x_{2}-x_{2}^{3}+x_{3}\left ( x_{1}^{2}+x_{2}^{2}-1-x_{1}+x_{1}x_{2}+x_{2}^{3}\right ) \\ x_{1}-x_{3}\left ( x_{1}-x_{2}+2x_{1}x_{2}\right ) \\ \left ( x_{3}-1\right ) \left ( x_{3}+2x_{3}x_{2}^{2}+x_{3}^{3}\right ) \end {pmatrix} =\begin {pmatrix} f_{1}\\ f_{2}\\ f_{3}\end {pmatrix} \] step 2 Find first integral
TO DO.
References