function [y] = polynomial(coeff, x)
n = length(coeff);
L = length(x);
x0 = ones(1,L);
y = zeros(1,L);
for i = 1:n
    y = y + coeff(i)*x0;
    x0=x0.*x;
end
end

