LIB "tst.lib";tst_init();
proc test5prim (int p)
{
  ring r= p, x(1..11), dp;
  poly d1= x(1);
  poly d2= x(2);
  poly f1= x(1);
  poly f2= x(2);
  poly g1= x(1);
  poly g2= x(2);
  poly d= d1*d2 - 1;
  poly f= f1*f2 + 2;
  poly g= g1*g2 - 2;
  f= f*d;
  g= g*d;
  poly testgcd= gcd (f,g);
  testgcd= testgcd/leadcoef (testgcd);
  testgcd == d;
  int i;
  for (i= 3; i <= 11; i++)
  {
    d2= d2*x(i);
    f2= f2*x(i);
    g2= g2*x(i);
    
    d= d1*d2 - 1;
    f= f1*f2 + 2;
    g= g1*g2 - 2;
    f= f*d;
    g= g*d;
    testgcd= gcd (f,g);
    testgcd= testgcd/leadcoef (testgcd);
    testgcd == d;
  }
}

test5prim (3);
tst_status();
test5prim (43051);
tst_status();
test5prim (0);
tst_status();
tst_status(1);
$