Abstract:Hamiltonian dynamical systems can be studied from a variety of viewpoints. Our intention in this paper is to show some examples of usage of two Maxima packages for symbolical and numerical analysis (\texttt{pdynamics} and \texttt{poincare}, respectively), along with the set of scripts \KeTCindy\ for obtaining the \LaTeX\ code corresponding to graphical representations of Poincaré sections, including animation movies.
What problem does this paper attempt to address?
This paper aims to solve several key problems encountered in the research of Hamiltonian dynamical systems, specifically including:
1. **Determining closed stable orbits**: In the theory of dynamical systems, an important goal is to find the possible closed and stable orbits. These orbits are very important for understanding the long - term behavior of the system.
2. **Calculating adiabatic invariants**: Another main goal is to calculate adiabatic invariants, which are crucial for describing the behavior of the system under slowly changing conditions.
3. **Dealing with perturbed Hamiltonian systems**: The paper pays special attention to the case where the Hamiltonian function \(H\) is a perturbation of an integrable Hamiltonian function \(H_0\), that is, \(H = H_0+\sum_{j = 1}^n\varepsilon_jH_j\). A common method for studying such perturbed systems is to write them in the so - called "normal form", that is, to represent them as a formal series:
\[
H=\sum_{j = 0}^{\infty}\varepsilon^jN_j
\]
where \(N_0 = H_0\), and each \(N_j\) commutes with the unperturbed Hamiltonian function \(H_0\), that is, satisfies the Poisson bracket relation:
\[
\{H_0,N_j\}=0
\]
4. **Constructing Poincaré sections**: In order to study the existence and stability of closed orbits, Poincaré sections provide a direct and intuitive method. However, the analytical expressions of Poincaré sections are usually difficult to obtain, so numerical methods are required for calculation.
5. **Improving the quality of graphic output**: Traditional computer algebra systems (CAS) have deficiencies in graphic output, especially when dealing with phase - space portraits. For this reason, the paper introduces the KETCindy tool, which can generate high - quality graphics and can be inserted into LaTeX documents, and can even generate animations.
By combining Maxima's symbolic computing and numerical analysis toolkits (such as pdynamics and poincare), as well as KETCindy's graphic processing capabilities, the paper shows how to efficiently study the symbolic, numerical, and graphic characteristics of Hamiltonian dynamical systems. This not only improves the computational efficiency but also significantly improves the quality of graphic output, making the research results more intuitive and easier to understand.