LIB "primdecint.lib";
ring R=integer,(a,b,c,d),dp;
ideal I1=9,a,b;
ideal I2=3,c;
ideal I3=11,2a,7b;
ideal I4=13a2,17b4;
ideal I5=9c5,6d5;
ideal I6=17,a15,b15,c15,d15;
ideal I=intersectZ(I1,I2);
I=intersectZ(I,I3);
I=intersectZ(I,I4);
I=intersectZ(I,I5);
I=intersectZ(I,I6);
radicalZ(I);
==> _[1]=102bd
==> _[2]=78ad
==> _[3]=51bc
==> _[4]=39ac
==> _[5]=6abd
==> _[6]=3abc
ideal J=intersectZ(ideal(17,a),ideal(17,a2,b));
radicalZ(J);
==> _[1]=17
==> _[2]=a