Ticket #178: in_singular3

option(redSB);
ring R=0,(x(1..16),y(1..16)),lp;
poly f1 = x(1);
poly f2 = y(1);
poly f3 = x(2) -1;
poly f4 = y(2);
poly f5 = (x(1)-x(4))^2 + (y(1)-y(4))^2 - 1;
poly f6 = (x(2)-x(4))^2 + (y(2)-y(4))^2 - 1;
poly f7 = (x(2)-x(5))^2 + (y(2)-y(5))^2 - 1;
poly f8 = (x(4)-x(5))^2 + (y(4)-y(5))^2 - 1;
poly f9 = (x(2)-x(3))^2 + (y(2)-y(3))^2 - 1;
poly f10 = (x(5)-x(3))^2 + (y(5)-y(3))^2 - 1;
poly f11 = (x(1)-x(6))^2 + (y(1)-y(6))^2 - 1;
poly f12 = (x(6)-x(7))^2 + (y(6)-y(7))^2 - 1;
poly f13 = (x(7)-x(3))^2 + (y(7)-y(3))^2 - 1;
poly f14 = (x(12)-x(3))^2 + (y(12)-y(3))^2 - 1;
poly f15 = (x(7)-x(12))^2 + (y(7)-y(12))^2 - 1;
poly f16 = (x(7)-x(11))^2 + (y(7)-y(11))^2 - 1;
poly f17 = (x(12)-x(11))^2 + (y(12)-y(11))^2 - 1;
poly f18 = (x(11)-x(13))^2 + (y(11)-y(13))^2 - 1;
poly f19 = (x(12)-x(13))^2 + (y(12)-y(13))^2 - 1;
poly f20 = (x(10)-x(6))^2 + (y(10)-y(6))^2 - 1;
poly f21 = (x(10)-x(11))^2 + (y(10)-y(11))^2 - 1;
poly f22 = (x(1)-x(8))^2 + (y(1)-y(8))^2 - 1;
poly f23 = (x(6)-x(8))^2 + (y(6)-y(8))^2 - 1;
poly f24 = (x(10)-x(8))^2 + (y(10)-y(8))^2 - 1;
poly f25 = (x(9)-x(8))^2 + (y(9)-y(8))^2 - 1;
poly f26 = (x(10)-x(9))^2 + (y(10)-y(9))^2 - 1;
poly f27 = (x(14)-x(9))^2 + (y(14)-y(9))^2 - 1;
poly f28 = (x(14)-x(13))^2 + (y(14)-y(13))^2 - 1;
poly f29 = (x(15)-x(9))^2 + (y(15)-y(9))^2 - 1;
poly f30 = (x(15)-x(14))^2 + (y(15)-y(14))^2 - 1;
poly f31 = (x(15)-x(16))^2 + (y(15)-y(16))^2 - 1;
poly f32 = (x(16)-x(14))^2 + (y(16)-y(14))^2 - 1;
poly f33 = (x(16)-x(13))^2 + (y(16)-y(13))^2 - 1;
poly f34 = (x(1)-x(9))^2 + (y(1)-y(9))^2 - 4;
poly f35 = (x(9)-x(13))^2 + (y(9)-y(13))^2 - 4;
poly f36 = (x(13)-x(3))^2 + (y(13)-y(3))^2 - 4;
poly f37 = (x(1)-x(10))^2 + (y(1)-y(10))^2 - 3;
poly f38 = (x(9)-x(6))^2 + (y(9)-y(6))^2 - 3;
poly f39 = (x(16)-x(9))^2 + (y(16)-y(9))^2 - 3;
poly f40 = (x(15)-x(13))^2 + (y(15)-y(13))^2 - 3;
poly f41 = (x(11)-x(3))^2 + (y(11)-y(3))^2 - 3;
poly f42 = (x(7)-x(13))^2 + (y(7)-y(13))^2 - 3;
poly f43 = 2*x(4)-1;
poly f44 = 4*y(4)-3;
poly f45 = 2*x(5)-3;
poly f46 = y(4)-y(5);
poly f47 = x(3)-2;
poly f48 = y(3);
poly f49 = x(13)-9;
poly f50 = (x(7)-x(6))*(y(11)-y(10))-(x(11)-x(10))*(y(7)-y(6));
poly f51 = (x(10)-x(6))*(y(11)-y(7))-(x(11)-x(7))*(y(10)-y(6));
ideal I=f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12,f13,f14,f15,f16,f17,f18,f19,f20,f21,f22,f23,f24,f25,f26,f27,f28,f29,f30,f31,f32,f33,f34,f35,f36,f37,f38,f39,f40,f41,f42,f43,f45,f46,f47,f48;
list L = facstd(I);
size(L);
L;