1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | option(prot); |
---|
5 | //------------------------------------------------------------ |
---|
6 | |
---|
7 | ring R=(0,w),(x1,y1,z1,x2,y2,z2,x3,y3,z3,l,d,xs,ys,zs,x,y,z), Dp; |
---|
8 | minpoly=w^2-3; |
---|
9 | poly s1=x1^2+y1^2+(2-z1)^2-9; |
---|
10 | poly s2=(-w-x2)^2+y2^2+(-1-z2)^2-9; |
---|
11 | poly s3=(w-x3)^2+y3^2+(-1-z3)^2-9; |
---|
12 | poly s4=(x-x1)*x+(y-y1)*y+(z-z1)*z; |
---|
13 | poly s5=(x-x2)*x+(y-y2)*y+(z-z2)*z; |
---|
14 | poly s6=(x-x3)*x+(y-y3)*y+(z-z3)*z; |
---|
15 | poly s7=(x1-x2)^2+(y1-y2)^2+(z1-z2)^2-12; |
---|
16 | poly s8=(x3-x2)^2+(y3-y2)^2+(z3-z2)^2-12; |
---|
17 | poly s9=(x1-x3)^2+(y1-y3)^2+(z1-z3)^2-12; |
---|
18 | poly s10=(x1+x2+x3)-3*x; |
---|
19 | poly s11=(y1+y2+y3)-3*y; |
---|
20 | poly s12=(z1+z2+z3)-3*z; |
---|
21 | poly s13=x^2-xs^2; |
---|
22 | poly s14=y^2-ys^2; |
---|
23 | poly s15=z^2-zs^2; |
---|
24 | poly s16=x^2+y^2+z^2-l^2; |
---|
25 | poly s17=y^2+z^2-d^2; |
---|
26 | poly s18=x-x1; |
---|
27 | poly s19=-3*d*ys-zs*6-3*d*y1; |
---|
28 | poly s20=-3*d*zs+ys*6-3*d*z1; |
---|
29 | poly s21=-2*l*xs+d*2*w+2*l*x2; |
---|
30 | poly s22=-6*ys*xs*w-6*d*l*ys+6*zs*l-6*d*l*y2; |
---|
31 | poly s23=-6*zs*xs*w-6*zs*d*l-6*ys*l-6*d*l*z2; |
---|
32 | poly s24=2*w*d+2*xs*l-2*l*x3; |
---|
33 | poly s25=6*w*ys*xs-6*ys*d*l+6*zs*l-6*d*l*y3; |
---|
34 | poly s26=6*w*zs*xs-6*zs*d*l-6*ys*l-6*d*l*z3; |
---|
35 | ideal I=s1,s2,s3,s4,s5,s6,s7,s8,s9,s10,s11,s12,s13,s14,s15,s16,s17,s18,s19,s20,s21,s22,s23,s24,s25,s26; |
---|
36 | ideal J=groebner(I); |
---|
37 | J; |
---|
38 | kill R; |
---|
39 | |
---|
40 | tst_status(1);$ |
---|