 ring r=3,(x,y,z),dp;
poly f=(xy)^3*(x+z)*(yz);
sqrfree(f);
==> [1]:
==> _[1]=1
==> _[2]=xy+xzyz+z2
==> _[3]=x+y
==> [2]:
==> 1,1,3
sqrfree(f,1);
==> _[1]=xy+xzyz+z2
==> _[2]=x+y
sqrfree(f,2);
==> [1]:
==> _[1]=xy+xzyz+z2
==> _[2]=x+y
==> [2]:
==> 1,3
sqrfree(f,3);
==> x2yxy2x2zxyzy2zxz2+yz2
