code: def density(self, x, y, n_sersic, r_eff, k_eff, center_x=0, center_y=0):
          
          raise ValueError("not implemented! Use a Multi-Gaussian-component decomposition.")
docs: de-projection of the Sersic profile based on
          Prugniel & Simien (1997)
          :return: