source:git/Tst/Short/factorizeQa_s.tst@276c3f

spielwiese
Last change on this file since 276c3f was faf58f, checked in by Martin Lee <martinlee84@…>, 12 years ago
• Property mode set to `100644`
File size: 3.5 KB
Line
1LIB "tst.lib";
2tst_init();
3
4// from gorzel
5ring S=(0,j),(u(1..2),x(1..2),z(1..2)),dp;
6minpoly=j2+1;
7poly P=u(1)^4*z(1)^2-2*u(1)^3*x(1)*z(1)^2
8+2*u(1)*x(1)^3*z(1)^2-x(1)^4*z(1)^2
9+2*u(1)*x(1)*x(2)^2*z(1)^2
10-2*x(1)^2*x(2)^2*z(1)^2-x(2)^4*z(1)^2
11-4*u(1)^3*x(2)*z(1)*z(2)
12+4*u(1)*x(1)^2*x(2)*z(1)*z(2)
13+4*u(1)*x(2)^3*z(1)*z(2)
14+u(1)^4*z(2)^2+2*u(1)^3*x(1)*z(2)^2
15-2*u(1)*x(1)^3*z(2)^2-x(1)^4*z(2)^2
16-2*u(1)*x(1)*x(2)^2*z(2)^2
17-2*x(1)^2*x(2)^2*z(2)^2-x(2)^4*z(2)^2;
18
19factorize(P);
20
21// ticket #330
22ring rc =(0,c),(x,y),dp;
23minpoly = c2+1/4c-1/16;
24poly fx = x2+(-1/4c)*x+(-1/256c+1/1024);
25factorize(fx);
26poly fxy = homog(fx,y);
27fxy;
28factorize(fxy);
29
30//ticket #236
31ring r = (0,a),(x,y),dp;
32minpoly = a2+1;
33poly f = (x-a)*(y-a);
34factorize(f);
35
36//ticket #89
37ring ri =(0,i),(a,b,c),dp;
38minpoly = i2+1;
39poly f = (a+b)^2+c^2;
40factorize(f);
41facstd(f);
42poly g = (a+c)^2+b^2;
43factorize(g);
44facstd(g);
45
46
47kill r;
48ring r= (0,a),(x,y,z,s),dp;
49minpoly= a2+1;
50poly f= 8730/37*a*x^5*y^14*z^5*s^5+17460/37*a*x^8*y^11*z^5*s^5+8730/37*a*x^11*y^8*z^5*s^5+4860/37*x^4*y^13*z^4*s^5+9720/37*x^7*y^10*z^4*s^5+4860/37*x^10*y^7*z^4*s^5-4559/48*a*x^4*y^13*z^5*s^4-4559/24*a*x^7*y^10*z^5*s^4-4559/48*a*x^10*y^7*z^5*s^4-423/8*x^3*y^12*z^4*s^4-423/4*x^6*y^9*z^4*s^4-423/8*x^9*y^6*z^4*s^4-9000/37*a*x^7*y^11*z^2*s^4-18000/37*a*x^10*y^8*z^2*s^4-9000/37*a*x^13*y^5*z^2*s^4-3492/49*a*x^5*y^15*z^5*s^3-3783/28*a*x^6*y^14*z^5*s^3-7081/37*x^6*y^14*z^5*s^3-6984/49*a*x^8*y^12*z^5*s^3-3783/14*a*x^9*y^11*z^5*s^3-14162/37*x^9*y^11*z^5*s^3-3492/49*a*x^11*y^9*z^5*s^3-3783/28*a*x^12*y^8*z^5*s^3-7081/37*x^12*y^8*z^5*s^3+2716/23*a*x^5*y^14*z^4*s^3-1944/49*x^4*y^14*z^4*s^3+3942/37*a*x^5*y^13*z^4*s^3-1053/14*x^5*y^13*z^4*s^3+5432/23*a*x^8*y^11*z^4*s^3-3888/49*x^7*y^11*z^4*s^3+7884/37*a*x^8*y^10*z^4*s^3-1053/7*x^8*y^10*z^4*s^3+2716/23*a*x^11*y^8*z^4*s^3-1944/49*x^10*y^8*z^4*s^3+3942/37*a*x^11*y^7*z^4*s^3-1053/14*x^11*y^7*z^4*s^3+1512/23*x^4*y^13*z^3*s^3+3024/23*x^7*y^10*z^3*s^3+1512/23*x^10*y^7*z^3*s^3+1175/12*a*x^6*y^10*z^2*s^3-2070/37*a*x^9*y^9*z^2*s^3+1175/6*a*x^9*y^7*z^2*s^3-4140/37*a*x^12*y^6*z^2*s^3+1175/12*a*x^12*y^4*z^2*s^3-2070/37*a*x^15*y^3*z^2*s^3-5940/37*a*x^5*y^10*z*s^3-11880/37*a*x^8*y^7*z*s^3-5940/37*a*x^11*y^4*z*s^3-6120/37*x^7*y^11*z^4*s^2-12240/37*x^10*y^8*z^4*s^2-6120/37*x^13*y^5*z^4*s^2+3600/49*a*x^7*y^12*z^2*s^2+975/7*a*x^8*y^11*z^2*s^2+7300/37*x^8*y^11*z^2*s^2+7200/49*a*x^10*y^9*z^2*s^2+1950/7*a*x^11*y^8*z^2*s^2+14600/37*x^11*y^8*z^2*s^2+1081/48*a*x^8*y^8*z^2*s^2+3600/49*a*x^13*y^6*z^2*s^2+975/7*a*x^14*y^5*z^2*s^2+7300/37*x^14*y^5*z^2*s^2+1081/24*a*x^11*y^5*z^2*s^2+1081/48*a*x^14*y^2*z^2*s^2-2800/23*a*x^7*y^11*z*s^2+517/8*a*x^4*y^9*z*s^2-5600/23*a*x^10*y^8*z*s^2+517/4*a*x^7*y^6*z*s^2-2800/23*a*x^13*y^5*z*s^2+517/8*a*x^10*y^3*z*s^2+799/12*x^6*y^10*z^4*s+799/6*x^9*y^7*z^4*s+799/12*x^12*y^4*z^4*s+828/49*a*x^9*y^10*z^2*s+897/28*a*x^10*y^9*z^2*s+1679/37*x^10*y^9*z^2*s+1656/49*a*x^12*y^7*z^2*s+897/14*a*x^13*y^6*z^2*s+3358/37*x^13*y^6*z^2*s+828/49*a*x^15*y^4*z^2*s+897/28*a*x^16*y^3*z^2*s+1679/37*x^16*y^3*z^2*s+2376/49*a*x^5*y^11*z*s+1287/14*a*x^6*y^10*z*s+4818/37*x^6*y^10*z*s-28*a*x^9*y^9*z*s+4752/49*a*x^8*y^8*z*s+1287/7*a*x^9*y^7*z*s+9636/37*x^9*y^7*z*s-56*a*x^12*y^6*z*s+2376/49*a*x^11*y^5*z*s+1287/14*a*x^12*y^4*z*s+4818/37*x^12*y^4*z*s-28*a*x^15*y^3*z*s-1848/23*a*x^5*y^10*s-3696/23*a*x^8*y^7*s-1848/23*a*x^11*y^4*s+2448/49*x^7*y^12*z^4-4964/37*a*x^8*y^11*z^4+663/7*x^8*y^11*z^4+4896/49*x^10*y^9*z^4-9928/37*a*x^11*y^8*z^4+1326/7*x^11*y^8*z^4+2448/49*x^13*y^6*z^4-4964/37*a*x^14*y^5*z^4+663/7*x^14*y^5*z^4-1904/23*x^7*y^11*z^3-3808/23*x^10*y^8*z^3-1904/23*x^13*y^5*z^3;
51f= f*f;
52factorize (f);
53
54tst_status(1); \$
Note: See TracBrowser for help on using the repository browser.