//SINGULAR Example4.4.4

//========================= we need==============================
option(redSB);            //to obtain a reduced normal form 
ring R=0,(x,y),lp;            
ideal a1=x;               //preparation of the example 
ideal a2=y2+2y+1,x-1;
ideal a3=y2-2y+1,x-1;
ideal I=intersect(a1,a2,a3);
poly h=y4-2y2+1;       //the lcm of the leading coefficients

//===============================================================

dim(std(I+ideal(h)));