
Figure 1.
Example (24) with solving $ \Delta v(x,y) = 1 $. Chainrecurrent set (top) approximated by the set $ \{(x,y)\, \mid\, \Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (middle) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $. $ \Delta v $ is approximately zero on the chainrecurrent set (origin) and negative everywhere else. Bottom: Constructed complete Lyapunov function $ v(x,y) $, which has a minimum at the origin

Figure 2.
Example (24) with equality and inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ (middle). Again, the orbital derivative $ \Delta v $ is correctly approximated being zero on the chainrecurrent set (origin) and negative everywhere else. Bottom: Constructed complete Lyapunov function $ v(x,y) $, which has a minimum at the origin. The point with equality constraint was $ (0.5,0) $

Figure 3.
Example (24) with inequality constraints. Chain recurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ (middle). $ \Delta v $ is approximately zero on the chainrecurrent set (origin) and negative everywhere else. Bottom: The constructed complete Lyapunov function $ v(x,y) $, which has a minimum at the origin

Figure 4.
Example (25) with solving $ \Delta v(x,y) = 1 $. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. The approximated chainrecurrent set includes the equilibria at the origin and $ (\pm 1,0) $, but is much larger, in particular around the origin. Bottom: Constructed complete Lyapunov function $ v(x,y) $, which has a saddle point at the origin and a local maximum at the unstable equilibria $ (\pm 1,0) $

Figure 5.
Example (25) with equality and inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. $ \Delta v $ is approximately zero on the chainrecurrent set, consisting of three equilibria at the origin and $ (\pm 1,0) $, and negative everywhere else. The approximation of the chainrecurrent set (the equilibria) is much better than when solving the equation $ \Delta v(x,y) = 1 $. Bottom: Constructed complete Lyapunov function $ v(x,y) $, which has a saddle point at the origin and local maxima at the unstable equilibria $ (\pm 1,0) $. Note that they have different levels, which is due to the extra point with equality constraint at $ (0.5,0) $, resulting in an unsymmetric approximation

Figure 6.
Example (25) with inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. $ \Delta v $ is approximately zero on the chainrecurrent set, consisting of three equilibria at the origin and $ (\pm 1,0) $, and negative everywhere else. The approximation of the chainrecurrent set (the equilibria) is much better than when solving the equation $ \Delta v(x,y) = 1 $. Bottom: Constructed complete Lyapunov function $ v(x,y) $, which has a saddle point at the origin and local maxima at the unstable equilibria $ (\pm 1,0) $

Figure 7.
Example (26) with solving $ \Delta v(x,y) = 1 $. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $. The approximated chainrecurrent set does not resemble the Hénon attractor very well, neither using the orbital derivative nor as the local minimum of the constructed function

Figure 8.
Example (26) with equality and inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. The characteristic shape of the Hénon attractor is clearly visible. Bottom: Constructed complete Lyapunov function $ v(x,y) $ with a local minimum at the Hénon attractor. The point with equality constraint was $ (0.5,0) $

Figure 9.
Example (26) with inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. The characteristic shape of the Hénon attractor is clearly visible. Bottom: Constructed complete Lyapunov function $ v(x,y) $ with a local minimum at the Hénon attractor

Figure 10.
Example (27) with solving $ \Delta v(x,y) = 1 $. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $. The approximated chainrecurrent set shows the Hénon repeller better than the Hénon attractor in the previous example, but still not very clearly. It is not clearly visible as local maximum of the constructed function either

Figure 11.
Example (27) with equality and inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $, showing the Hénon repeller as a local maximum. The repeller is clearly visible in all figures. The point with equality constraint was $ (0.5,0) $

Figure 12.
Example (27) with inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $, showing the Hénon repeller as a local maximum. The repeller is clearly visible in all figures

Figure 13.
Example (28) with solving $ \Delta v(x,y) = 1 $. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $. The approximated chainrecurrent set shows the attractor relatively well in the orbital derivative, but not very clearly as local minimum of the constructed function

Figure 14.
Example (28) with equality and inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $, showing the attractor as a local minimum. The attractor is clearer than in the previous method, both using the orbital derivative and as local minimum of the constructed function. The point with equality constraint was $ (0.5,0) $, where the orbital derivative is fixed at $ 1 $

Figure 15.
Example (28) with inequality constraints. Chainrecurrent set (top) approximated by the set $ \{(x,y)\,\mid\,\Delta v(x,y)\ge \gamma\} $, see Table 2, and the orbital derivative (second) $ \Delta v(x,y) $ of the constructed complete Lyapunov function $ v $ over the chainrecurrent set. The third figure shows the orbital derivative in a larger set. Bottom: Constructed complete Lyapunov function $ v(x,y) $, showing the attractor as a local minimum. The attractor is clearer than in the first method, both using the orbital derivative and as local minimum of the constructed function

Figure 16.
Example (29) with solving $ \Delta v(x,y,z) = 1 $. Top: Chainrecurrent set approximated by the set $ \{(x,y,z)\,\mid\,\Delta v(x,y,z)\ge \gamma\} $, see Table 2. The other figures show projections of this set: projections to the $ xy $ (second), $ yz $ (third) and $ xz $plane (bottom)

Figure 17.
Example (29) with equalityinequality constrains. Top: Chainrecurrent set approximated by the set $ \{(x,y,z)\,\mid\,\Delta v(x,y,z)\ge \gamma\} $, see Table 2. The other figures show projections of this set: projections to the $ xy $ (second), $ yz $ (third) and $ xz $plane (bottom). The figures are not as good as with the previous method. The point with equality constraint is $ (0.4,0.4,0) $

Figure 18.
Example (29) with inequality constrains. Top: Chainrecurrent set approximated by the set $ \{(x,y,z)\,\mid\,\Delta v(x,y,z)\ge \gamma\} $, see Table 2. The other figures show projections of this set: projections to the $ xy $ (second), $ yz $ (third) and $ xz $plane (bottom)