@(p,q)poly([roots(p);roots(q)])*p(1)*q(1) ans([1,1], [1,4,6,4,1])