gcd - gcd calculation
computes the gcd of components of p and a unimodular matrix (with polynomial inverse) U, with minimal degree such that
p*U=[0 ... 0 pgcd]
s=poly(0,'s'); p=[s,s*(s+1)^2,2*s^2+s^3]; [pgcd,u]=gcd(p); p*u