def METHOD_NAME(self, x, y, n_sersic, R_sersic, k_eff, e1, e2, center_x=0, center_y=0):
    raise Exception("not yet implemented")
    # phi_G, q = param_util.ellipticity2phi_q(e1, e2)
    #
    # if isinstance(x, float) and isinstance(y, float):
    #
    #     x_, y_ = self._coord_rotate(x, y, phi_G, center_x, center_y)
    #     integral = quad(self._integrand_I, 0, 1, args=(x_, y_, q, n_sersic, R_sersic, k_eff, center_x, center_y))[0]
    #
    # else:
    #
    #     assert isinstance(x, np.ndarray) or isinstance(x, list)
    #     assert isinstance(y, np.ndarray) or isinstance(y, list)
    #     x = np.array(x)
    #     y = np.array(y)
    #     shape0 = x.shape
    #     assert shape0 == y.shape
    #
    #     if isinstance(phi_G, float) or isinstance(phi_G, int):
    #         phiG = np.ones_like(x) * float(phi_G)
    #         q = np.ones_like(x) * float(q)
    #     integral = []
    #     for i, (x_i, y_i, phi_i, q_i) in \
    #             enumerate(zip(x.ravel(), y.ravel(), phiG.ravel(), q.ravel())):
    #
    #         integral.append(quad(self._integrand_I, 0, 1, args=(x_, y_, q, n_sersic,
    #                                                             R_sersic, k_eff, center_x, center_y))[0])
    #
    #
    # return 0.5 * q * integral