| ring r=3,(x,y,z),dp;
poly f=(x-y)^3*(x+z)*(y-z);
sqrfree(f);
==> [1]:
==> _[1]=1
==> _[2]=-xy+xz-yz+z2
==> _[3]=-x+y
==> [2]:
==> 1,1,3
sqrfree(f,1);
==> _[1]=-xy+xz-yz+z2
==> _[2]=-x+y
sqrfree(f,2);
==> [1]:
==> _[1]=-xy+xz-yz+z2
==> _[2]=-x+y
==> [2]:
==> 1,3
sqrfree(f,3);
==> x2y-xy2-x2z-xyz-y2z-xz2+yz2
|