gcd(p,diff(p,x)) gives x-sqrt(2) with gcd = subres or red and algebraic = true On Fri, Jan 11, 2013 at 12:40 PM, Raymond Toy <toy.raymond at gmail.com> wrote: > p:expand((x-sqrt(2))^2*(x^2+2*x+2)); > gcd(p,diff(p,x)) >