function y=poly_nest(p,x)

% Evaluate polynomal p(x) = p(1) x^{n-1} + p(2) x^{n-2} + ... + p(n)

n=length(p);
y= p(1);
for i=n-1:-1:1
   y = x*y + p(n-i+1);
end