primdec.lib
LIB "primdec.lib"; ring R=(5,t,s),(x,y,z),dp; poly f=(x^25-t*x^5+t)*(x^3+s); sep(f,1); ==> [1]: ==> x8+(s5)*x5+(-t)*x4+(t)*x3+(-ts5)*x+(ts5) ==> [2]: ==> 1