Center Manifold Example

This example demonstrates how to compute center manifolds in the Circular Restricted Three-Body Problem (CR3BP).

Computing center manifolds in the CR3BP
 1"""Example script: computing the centre manifold Hamiltonian for the Earth-Moon hiten.system.
 2
 3Run with
 4    python examples/center_manifold.py
 5"""
 6
 7import os
 8import sys
 9
10sys.path.append(os.path.join(os.path.dirname(__file__), "..", "src"))
11
12from hiten import System
13
14
15def main() -> None:
16    """Compute and display the centre-manifold Hamiltonian."""
17    system = System.from_bodies("sun", "earth")
18    l_point = system.get_libration_point(1)
19    cm = l_point.get_center_manifold(degree=5)
20    cm.compute()
21
22    print(cm.coefficients())
23
24if __name__ == "__main__":
25    main()