def METHOD_NAME(self):
    x = np.linspace(0.01, 2, 10)
    y = np.zeros_like(x)
    kwargs = {"alpha_Rs": 2, "Rs": 2, "center_x": 0, "center_y": 0}
    f_nfw = self.nfw.function(x, y, **kwargs)
    f_cse = self.nfw_cse.function(x, y, e1=0, e2=0, **kwargs)
    npt.assert_almost_equal(f_cse, f_nfw, decimal=5)
    f_cse_low = self.nfw_cse_low_accuracy.function(x, y, e1=0, e2=0, **kwargs)
    npt.assert_almost_equal(f_cse_low / f_nfw, 1, decimal=3)