x = -4:1:4
n = length(x);
syms t
P = 1;
for i = 1:n
    P = P * (t-x(i));
end
ezplot(P,[-4 4])