LIB "realrad.lib";
ring r1=0,(x,y,z),lp;
//dimension 0
ideal i0=(x2+1)*(x3-2),(y3-2)*(y2+y+1),z3+2;
//dimension 1
ideal i1=(y3+3y2+y+1)*(y2+4y+4)*(x2+1),(x2+y)*(x2-y2)*(x2+2xy+y2)*(y2+y+1);
ideal i=intersect(i0,i1);
realrad(i);
==> _[1]=xy4+5xy3+7xy2+3xy+2x+y4z+5y3z+7y2z+3yz+2z
==> _[2]=-236x4y-590x4+236x2y3+354x2y2-590x2y+873y5z+573y5+873y4z2+4938y4z+31\
   01y4+4365y3z2+8976y3z+4601y3+6111y2z2+6630y2z+1719y2+2619yz2+3465yz+1146y\
   +1746z2+1146z
==> _[3]=68676x4yz-844x4y+171690x4z-2110x4-68676x2y3z+844x2y3-103014x2y2z+126\
   6x2y2+171690x2yz-2110x2y-316581y5-385257y4z-1582061y4-1754595y3z-2213957y\
   3-2216067y2z-949743y2-949743yz-633162y-633162z
==> _[4]=3x4y2+8x4y-x4+10x2y3+30x2y2+8x2y+6x2-3y4z2+3y4z+7y4-15y3z2+15y3z+22y\
   3-21y2z2+21y2z+9y2-9yz2+9yz+6y-6z2+6z
==> _[5]=y4z3+2y4+5y3z3+10y3+7y2z3+14y2+3yz3+6y+2z3+4