//y^2-13*y+30 clear; clc; close; y=poly(0,'y'); p=y^2-13*y+30; factors(p)