# Ticket #178: in_singular3

File in_singular3, 2.4 KB (added by seelisch, 14 years ago)
Line
1option(redSB);
2ring R=0,(x(1..16),y(1..16)),lp;
3poly f1 = x(1);
4poly f2 = y(1);
5poly f3 = x(2) -1;
6poly f4 = y(2);
7poly f5 = (x(1)-x(4))^2 + (y(1)-y(4))^2 - 1;
8poly f6 = (x(2)-x(4))^2 + (y(2)-y(4))^2 - 1;
9poly f7 = (x(2)-x(5))^2 + (y(2)-y(5))^2 - 1;
10poly f8 = (x(4)-x(5))^2 + (y(4)-y(5))^2 - 1;
11poly f9 = (x(2)-x(3))^2 + (y(2)-y(3))^2 - 1;
12poly f10 = (x(5)-x(3))^2 + (y(5)-y(3))^2 - 1;
13poly f11 = (x(1)-x(6))^2 + (y(1)-y(6))^2 - 1;
14poly f12 = (x(6)-x(7))^2 + (y(6)-y(7))^2 - 1;
15poly f13 = (x(7)-x(3))^2 + (y(7)-y(3))^2 - 1;
16poly f14 = (x(12)-x(3))^2 + (y(12)-y(3))^2 - 1;
17poly f15 = (x(7)-x(12))^2 + (y(7)-y(12))^2 - 1;
18poly f16 = (x(7)-x(11))^2 + (y(7)-y(11))^2 - 1;
19poly f17 = (x(12)-x(11))^2 + (y(12)-y(11))^2 - 1;
20poly f18 = (x(11)-x(13))^2 + (y(11)-y(13))^2 - 1;
21poly f19 = (x(12)-x(13))^2 + (y(12)-y(13))^2 - 1;
22poly f20 = (x(10)-x(6))^2 + (y(10)-y(6))^2 - 1;
23poly f21 = (x(10)-x(11))^2 + (y(10)-y(11))^2 - 1;
24poly f22 = (x(1)-x(8))^2 + (y(1)-y(8))^2 - 1;
25poly f23 = (x(6)-x(8))^2 + (y(6)-y(8))^2 - 1;
26poly f24 = (x(10)-x(8))^2 + (y(10)-y(8))^2 - 1;
27poly f25 = (x(9)-x(8))^2 + (y(9)-y(8))^2 - 1;
28poly f26 = (x(10)-x(9))^2 + (y(10)-y(9))^2 - 1;
29poly f27 = (x(14)-x(9))^2 + (y(14)-y(9))^2 - 1;
30poly f28 = (x(14)-x(13))^2 + (y(14)-y(13))^2 - 1;
31poly f29 = (x(15)-x(9))^2 + (y(15)-y(9))^2 - 1;
32poly f30 = (x(15)-x(14))^2 + (y(15)-y(14))^2 - 1;
33poly f31 = (x(15)-x(16))^2 + (y(15)-y(16))^2 - 1;
34poly f32 = (x(16)-x(14))^2 + (y(16)-y(14))^2 - 1;
35poly f33 = (x(16)-x(13))^2 + (y(16)-y(13))^2 - 1;
36poly f34 = (x(1)-x(9))^2 + (y(1)-y(9))^2 - 4;
37poly f35 = (x(9)-x(13))^2 + (y(9)-y(13))^2 - 4;
38poly f36 = (x(13)-x(3))^2 + (y(13)-y(3))^2 - 4;
39poly f37 = (x(1)-x(10))^2 + (y(1)-y(10))^2 - 3;
40poly f38 = (x(9)-x(6))^2 + (y(9)-y(6))^2 - 3;
41poly f39 = (x(16)-x(9))^2 + (y(16)-y(9))^2 - 3;
42poly f40 = (x(15)-x(13))^2 + (y(15)-y(13))^2 - 3;
43poly f41 = (x(11)-x(3))^2 + (y(11)-y(3))^2 - 3;
44poly f42 = (x(7)-x(13))^2 + (y(7)-y(13))^2 - 3;
45poly f43 = 2*x(4)-1;
46poly f44 = 4*y(4)-3;
47poly f45 = 2*x(5)-3;
48poly f46 = y(4)-y(5);
49poly f47 = x(3)-2;
50poly f48 = y(3);
51poly f49 = x(13)-9;
52poly f50 = (x(7)-x(6))*(y(11)-y(10))-(x(11)-x(10))*(y(7)-y(6));
53poly f51 = (x(10)-x(6))*(y(11)-y(7))-(x(11)-x(7))*(y(10)-y(6));
54ideal 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;
55list L = facstd(I);
56size(L);
57L;