@property
def disc(P):
    """
        Discriminant of a real polynomial
        """
    ans = P.resultant(P.prime()) / P[-1]
    if P.isinteger():
        ans = int(ans.round())
    if P.deg % 4 in [0, 1]:
        return ans
    else:
        return -ans
