def test_rho_as_angle(): # arctan(100 km, 1 au) * 206264.806 = 0.13787950659645942 rho = rho_as_angle(100 * u.km, {'delta': 1 * u.au}) assert np.isclose(rho.to(u.arcsec).value, 0.13787950659645942)