@pytest.mark.parametrize(('primary_mass', 'secondary_mass', 'semimajor_axis', 'expected'), [p(100000000000.0, 100, 10, 76.9102, id='arbitrary period')])
def test_orbital_period(primary_mass, secondary_mass, semimajor_axis, expected):
    assert (OrbitalPeriod(primary_mass, secondary_mass, semimajor_axis, G).evalf() == pytest.approx(expected, 0.001))
@pytest.mark.parametrize(('primary_mass', 'secondary_mass', 'semimajor_axis', 'expected'), [p(100000000000.0, 100, 10, 76.9102, id='arbitrary period')])
def test_orbital_period(primary_mass, secondary_mass, semimajor_axis, expected):
    assert (OrbitalPeriod(primary_mass, secondary_mass, semimajor_axis, G).evalf() == pytest.approx(expected, 0.001))