//x^2+13*x+36 clear; clc; close; x=poly(0,'x'); p=x^2+13*x+36; factors(p)