def METHOD_NAME(self):
    geod_circle = self.geod.circle(40, 50, 500000, n_samples=3)
    assert_almost_equal(geod_circle,
                        np.array([[40., 54.49349757],
                                  [34.23766162, 47.60355349],
                                  [45.76233838, 47.60355349]]), decimal=5)