// interpola.sci

// devolve um vetor com os coeficientes
// do polinomio interpolador da tabela T.
// o grau do polinomio interpolador eh
// regido pelo tamanho da tabela.

function c = interpola (T)
  // importando as funcoes necessarias
  getf("gausspv.sci");
  
  [nl, nc] = size(T);
 
  // montando o sistema linear a ser resolvido
  A = zeros (nc, nc);
  
  for i = 1:nc do
    for j = 1:nc do
      A(j,i) = T(1,j)^(i - 1);
    end
  end
  
  b = T(2,:)';
  
  // resolvendo o sistema linear Ac = b;
  [p, B, c] = gausspv(A, b);
endfunction

