def METHOD_NAME(
    self, x, y, n_sersic, R_sersic, k_eff, phi_G, q, center_x=0, center_y=0
):
    alpha_x = (
        x
        * q
        * quad(
            self._integrand_J, 0, 1, args=(x, y, n_sersic, q, R_sersic, k_eff, 0)
        )[0]
    )
    alpha_y = (
        y
        * q
        * quad(
            self._integrand_J, 0, 1, args=(x, y, n_sersic, q, R_sersic, k_eff, 1)
        )[0]
    )
    return alpha_x, alpha_y