cadojo

using Plots
using Unitful
using DifferentialEquations
using GeneralAstrodynamics
using LazySets

# Find a periodic orbit near Earth (orbit, period)
Oₑ, Pₑ = halo(SunEarth; Az=100_000u"km", L=2)

# Find a periodic orbit near Jupiter (orbit, period)