function val=poly_val(poly,X)
%poly=[2 3 4] is 2+3*x+4*x^2
N=length(poly)-1;
[r,c]=size(X);
if r==1 
   n=c;
   X=X';
else
   n=r;
end;   
val=zeros(n,1);
xi=ones(n,1);
 for h=0:N
   val=val+poly(h+1)*xi;
   xi=xi.*X;
 end;