function [p] = eval_poly(a, x)
% a = [a0,a1,...,an] are coefficients of the degree n polynomial
% a0 + a1*x^1 + ... + an*x^n
% x is the value to evaluate polynomial at
% p is the value of the polynomial at x

n = length(a) - 1;
p = a(n+1);
for k = n : -1 : 1
   p = x*p + a(k);
end