p = [0.1 0 1 1];
r = roots(p)