def METHOD_NAME():
    params = ExoplanetOrbitDatabase.query_planet('HD 209458 b ')
    simbad_coords = SkyCoord(ra='22h03m10.77207s', dec='+18d53m03.5430s')
    sep = params['sky_coord'].separation(simbad_coords)
    print(sep, type(sep))
    assert abs(sep) < 5 * u.arcsec