T2 - CalcNum

x = [-2,-1,0,1,2,3,4];
y = [-14,0.5,3.1,0,-3,0,16];
n=length(x);
#inicializando a tabela
T = zeros(n,n);
#primeira coluna
T(:,1)=y;
#segunda coluna
for j = 2:n
  T(j,2)=(T(j,1)-T(j-1,1))/(x(j)-x(j-1));
  for k = 3:n
  #terceira coluna
    T(k,3)=(T(k,2)-T(k-1,2))/(x(k)-x(k-2));
    for m = 4:n
      #quarta coluna
      T(m,4)=(T(m,3)-T(m-1,3))/(x(m)-x(m-3));
      for q = 5:n
        #quinta coluna
        T(q,5)=(T(q,4)-T(q-1,4))/(x(q)-x(q-4));
        for r = 6:n
          #sexta coluna
          T(r,6)=(T(r,5)-T(r-1,5))/(x(r)-x(r-5));
          for s = n
            #setima coluna
            T(s,7)=(T(s,6)-T(s-1,6))/(x(s)-x(1));
          endfor
        endfor
      endfor
    endfor
  endfor
endfor
T

p = T(1,1) + T(2,2)*(x-x(1)) + T(3,3)*(x-x(1))*(x-x(2)) + T(4,4)*(x-x(1))*(x-x(2))*(x-x(3)) + T(5,5)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4)) + T(6,6)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4))*(x-x(5)) + T(7,7)*(x-x(1))*(x-x(2))*(x-x(3))*(x-x(4))*(x-x(5))*(x-x(6))
p