allocatemem(20*10^6);
do(P,Q)=my(v=variable()); substvec(polresultant(P,Q), v, vector(#v,i,i));

P = Pol([k,c,d,e,f,g]);
Q = P';
do(P,Q)

P = Pol([1,b,c,d,e]);
Q = Pol([1,g,h,k,j]);
do(P,Q)

P = x^7 + x^3 + b*x^2 + c*x + d;
Q = x^7 + x^3 + f*x^2 + g*x + h;
do(P,Q)

P = x^20 + a*x^5 + b;
Q = x^20 + c*x^5 + d;
do(P,Q)

P = (x+a)^8;
Q = (x+z)^8;
do(P,Q)

P = x^50 + 5*a*x^40 + 4*a*x^30 + a*x^20 + 2*a*x^10 + 3*a;
Q = x^45 + 2*b*x^35 + b*x^25 + 4*b*x^15 + 5*b*x^5 + 6*b;
do(P,Q)

P = polcyclo(31); P = subst(P, x, a*x);
Q = polcyclo(29); Q = subst(Q, x, a*x);
do(P,Q)

P = 1 + Polrev( vector(101, j, b^(101-j)) );
Q = 1 + Polrev( vector(101, j, b^(j-1)) );
do(P,Q)

P = 1 + Polrev( vector(301,i,1) );
Q = 1 + Polrev( vector(201,i,i-1) );
do(P,Q)

P = 1 + Polrev( vector(301,j, (j-1)^5) );
Q = 1 + Polrev( vector(301,j, (j-1)^4) );
do(P,Q)

n = 110;
polrandom(d, fudge) = x^d + Polrev(vector(d+1, i, round(2^32*sin(i+fudge))));
P = polrandom(n,0);
Q = polrandom(n,1/2);
do(P,Q)

polresultant(-27673*x^2+322883101*x-1876122109136,x^4+4503599627370517)
P=-27673*x^2+322883101*x-1876122109136;
for(n=2,3,for(l=1,3,print(polresultant(P*x^l+1,n*x^4+4503599627370517))))

polresultant(z^1746-1, polcyclo(1973, z))

for(i=1,15,if(rnfequation(polcyclo(17,'y),x+y^i,1)[2]^17!=1,print("error",i)))

\\ #1233
polresultantext(x^2+1, x^2-1)
\\ #1234
polresultant(x^4+x^2+1, x^4-y, x, 0)

p1=x2*(x3-x4);p2=x2*(x3-2*x4);polresultant(p1,p2,x1)
polresultant(x,x,y)

A = x*y; B = (x+y)^2;
polresultant(A, B)
polresultantext(A, B, y)

\\ #1509
polresultantext(y^3+y,3*y^2+1,y)
poldisc(x^3 + Mod(2,3)*x^2)
poldisc(x^3 + Mod(2,3)*x^2+1)

norm(Mod(Pol(Mod(0,2)), 3*x))

polresultantext(x+1, x^2)
\\#1614
polresultant(1+x*z^2,1+y*z^4,z)
