Aitken Neville

function [] = aitken_neville(y)
% paulogp
% Metodo de Aitken Neville
% Necessita de um ficheiro chamado "aitken_neville_tabela.m" com os valores
% Formato dos valores: x (tab) y
% no exemplo para a tabela 2 - y = 1.11
% resultado: 0.92240...

    % le os pontos de uma tabela de valores
    % m = csvread('aitken_neville_tabela.m');
    m = importdata('aitken_neville_tabela.m', '\t');

    % importante neste caso
    format long

    % vector
    x = m(:, 1);
    xi = x';
    f = m(:, 2);
    fi = f';

    % inicializacao do vector
    d = [];

    n = length(xi)-1;

    for i = 1:n+1
        d(i) = y-xi(i);
    end

    for k = 1:n
        for j = 1:(n-k+1)
            fi(j) = ((d(j)*fi(j+1))-(d(j+k)*fi(j)))/(d(j)-d(j+k));
        end
    end

    % saida do vector
    disp(fi);
end


% ficheiro: aitken_neville_tabela.m
1.00    1.017452
1.05    0.971622
1.10    0.930208
1.15    0.892646
1.20    0.858464
1.25    0.827269
1.30    0.798724