ring r = (0,a,b), (s1,c1,s2,c2), lp; ideal i = a - c1*c2+s1*s2-c1, b-s1*c2-s2*c1-s1, c1^2+s1^2-1, c2^2+s2^2-1; std(i);