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