def METHOD_NAME(self, x, y, q, n_sersic, R_sersic, k_eff, u=1, power=1):
    b_n = self._sersic.b_n(n_sersic)
    elliptical_coord = self._elliptical_coord_u(x, y, u, q) ** power
    elliptical_coord *= R_sersic**-power
    exponent = -b_n * (elliptical_coord ** (1.0 / n_sersic) - 1)
    return k_eff * np.exp(exponent)