using PlanetOrbits

# This file defines and calculates the orbit for the experiment

exper1 = orbit(
    M = 2,
    a = 15,
    i = 0,
    Ω = 0,
    ω = 0,
    plx=50
)