function [p] = my_poly_eval(c,x,t)

%MY_POLY_EVAL evaluates the polynomial at the points (t_0,...,t_m)
%given the coefficients c=(co,...,c_n) of the newton
%interpolation and the nodes (x_0,...,x_n)

m=max(size(c));
p=c(m)*ones(size(t));
for k=m-1:-1:1
    p=c(k)+(t-x(k)).*p;
end


end

