1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | LIB "elim.lib"; |
---|
5 | |
---|
6 | option(redSB); |
---|
7 | option(redTail); |
---|
8 | |
---|
9 | ring R = (0),( k, l, m, n, o, p, q, r, a, b, c, d, e, f, g, h, i, j), lp; |
---|
10 | |
---|
11 | ideal I = // your ideal |
---|
12 | -1 - e*i*m + f*j*m + c*k*m - d*l*m + f*i*n + e*j*n - d*k*n - c*l*n + e*g*o - f*h*o - a*k*o + b*l*o - |
---|
13 | f*g*p - e*h*p + b*k*p + a*l*p - c*g*q + d*h*q + a*i*q - b*j*q + d*g*r + c*h*r - b*i*r - a*j*r, |
---|
14 | -(f*i*m) - e*j*m + d*k*m + c*l*m - e*i*n + f*j*n + c*k*n - d*l*n + f*g*o + e*h*o - b*k*o - a*l*o + |
---|
15 | e*g*p - f*h*p - a*k*p + b*l*p - d*g*q - c*h*q + b*i*q + a*j*q - c*g*r + d*h*r + a*i*r - b*j*r, |
---|
16 | -1 + a^2 + b^2 + c^2 + d^2 + e^2 + f^2, a*g + b*h + c*i + d*j + e*k + f*l, |
---|
17 | a*m + b*n + c*o + d*p + e*q + f*r, a*g + b*h + c*i + d*j + e*k + f*l, |
---|
18 | -1 + g^2 + h^2 + i^2 + j^2 + k^2 + l^2, g*m + h*n + i*o + j*p + k*q + l*r, |
---|
19 | a*m + b*n + c*o + d*p + e*q + f*r, g*m + h*n + i*o + j*p + k*q + l*r, |
---|
20 | -1 + m^2 + n^2 + o^2 + p^2 + q^2 + r^2, b*g - a*h + d*i - c*j + f*k - e*l, |
---|
21 | b*m - a*n + d*o - c*p + f*q - e*r, -(b*g) + a*h - d*i + c*j - f*k + e*l, |
---|
22 | h*m - g*n + j*o - i*p + l*q - k*r, -(b*m) + a*n - d*o + c*p - f*q + e*r, |
---|
23 | -(h*m) + g*n - j*o + i*p - l*q + k*r; |
---|
24 | |
---|
25 | I = groebner(I); // GB(I) |
---|
26 | |
---|
27 | elim(I, 1, nvars(R)); // eliminate all variables from I |
---|
28 | |
---|
29 | tst_status(1);$ |
---|