read recden; read MGCD; read PGCD; g := 12*x-11*y+23*z-13: abar := randpoly([x,y,z],dense): bbar := randpoly([x,y,z],dense): a := expand(g*abar): b := expand(g*bbar): A := convert(a,POLYNOMIAL,[x,y,z]): B := convert(b,POLYNOMIAL,[x,y,z]): MGCD(A,B);