function pv = eval_newton(x,p,t);

   pv = p(1); n = length(x); l = 1;
   for i=2:n
     l = l*( t-x(i-1));
     pv = pv + l*p(i); 
   end