option("redSB"); LIB "primdec.lib"; ring rpe=(0,ah,bh),(ch3,sh3,ch2,sh2,ch1,sh1,ch0,sh0),(dp(4),dp(4)); poly q0=ch0^2+sh0^2-1; poly q1=ch1^2+sh1^2-1; poly q2=ch2^2+sh2^2-1; poly q3=ch3^2+sh3^2-1; poly p1=ah*ch0+bh*ch1+ah*ch2+bh*ch3; poly p2=ah*sh0+bh*sh1+ah*sh2+bh*sh3; ideal i1=p1,p2,q0,q1,q2,q3; list l1=minAssGTZ(i1); ideal j0=l1[1]; ideal j1=l1[2]; print("eka");