ring R=0,(s,t,x(0..4)),dp; ideal I=x(0)-t4,x(1)-t3s,x(2)-t2s2,x(3)-ts3,x(4)-s4; ideal J=eliminate(I,st); J; //->J[1]=x(3)^2-x(2)*x(4) //->J[2]=x(2)*x(3)-x(1)*x(4) //->J[3]=x(1)*x(3)-x(0)*x(4) //->J[4]=x(2)^2-x(0)*x(4) //->J[5]=x(1)*x(2)-x(0)*x(3) //->J[6]=x(1)^2-x(0)*x(2) ring S=0,(x(0..4)),dp; ideal J=imap(R,J); J=std(J); hilb(J); // 1 t^0 // -6 t^2 // 8 t^3 // -3 t^4 // 1 t^0 // 3 t^1 // dimension = 1 // degree = 4 LIB"poly.lib"; hilbPoly(J); //->1,4 dim(J); //->2 degree(J); //->4