Add _eval_simplify helper for AlgebraicNumber

fixes #4401
Sorry...I'm not sollowing the logic. RootOf will only return algebraic numbers in the case of quadratics, binomials or linear functions. But can't you have algebraic numbers from higher polynomials?