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