gcdCoefficients -- gcd with coefficients

Description

gcdCoefficients(a,b) -- returns {d,r,s} so that a*r + b*s is the greatest common divisor d of a and b.

Works for integers or elements of polynomial rings.