rootsur.lib
LIB "rootsur.lib"; ring r = 0,x,dp; poly p = (x+2)*(x-1)*(x-5); allrealst(p); ==> 1 p = p*(x2+1); allrealst(p); ==> 0