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