Newtons Interpolation: Finding the value of polynomials from divided difference table

Question

So I have this function to create a divided difference table but i need to write another function to find the values of the polynomials and i am completely lost as to how to do it.

function [F] = divdiff(x,y)
n=length(x)

F=zeros(n,n)
    
    for i = 1:n
        F(i,1) = y(i)
    end
    
    for j = 2:n
        for i = j:n
            F(i,j) = (F(i,j-1) - F(i-1, j-1))./ (x(i)-x(i-j+1));
        end
    end

now i have `function v = polyvalue(a,x,t)

    pval = zeros(length(t),1);
   for ind_t = 1:length(t);
       t_val = t(ind_t);
       
       for ind_a = 1: length(a)
           xterm_val = 1,
           
           for ind_x = 1:(ind_a - 1)
               xterm_val = xterm_val .* (t_val -x(ind_x));
               pval(ind_t) = pval(ind_t) + a(ind_a) .* xterm_val
           end
       end
   end
   end`

which is not working