fn main() {
    let orbit = deimoss::physics::orbit::DeimosOrbit::new();
    println!("Semi-major axis: {:.0} m", orbit.semi_major_axis_m);
    println!("Eccentricity: {}", orbit.eccentricity);
    println!(
        "Orbital period: {:.4} days",
        orbit.orbital_period_s() / 86_400.0
    );
    println!("Mean motion: {:.2e} rad/s", orbit.mean_motion_rad_s());
    println!("Periapsis: {:.0} m", orbit.periapsis_m());
    println!("Apoapsis: {:.0} m", orbit.apoapsis_m());
    println!(
        "Speed at semi-major: {:.1} m/s",
        orbit.orbital_speed_m_s(orbit.semi_major_axis_m)
    );

    let rotation = deimoss::physics::rotation::DeimosRotation::new();
    println!("Synchronous: {}", rotation.synchronous);
    println!(
        "Equatorial speed: {:.4} m/s",
        rotation.equatorial_speed_m_s()
    );
}
