docs: def density(self, r, rho0, gamma):
          """
          computes the density
          :param x:
          :param y:
          :param rho0:
          :param a:
          :param s:
          :return:
          """
          rho = rho0 / r**gamma
          return rho