The big number of digits asked for leads to that. You may see what does the following: clc %eps format(16) x = poly(0,'x'); f = (x-1)^2, factors(f) format(19) x = poly(0,'x'); f = (x-1)^2, factors(f) Regards Samuel