Diferenças Divididas: script no Scinotes do Scilab

/* Diferenças divididas

Polinômio Interpolador de Newton
Script por Marcos Pinnto */
function [D]=Divdif(A)
    //F é uma matriz nx2 com os pontos (x,f(x))
    F=[];
    X=[];
    //M=[];
    for i=1:size(A,1)
        X(i,1) = A(i,1);
        F(i,1) = A(i,2);
    end
    for k=2:size(A,1)
        for i=1:(size(A,1)-k+1)
            F(i,k) = ( F(i+1,k-1)-F(i,k-1) )/( X(i+k-1,1)-X(i,1) );    
        end
    end
    D = F;
   // disp(D);
endfunction