@staticmethod
def density(r, rho0, e1=0, e2=0):
    """
        computes the density
        :param r: radius in angles
        :param rho0: density at angle=1
        :return: density at r
        """
    rho = rho0 / r ** 2
    return rho
