[Matlab] LiczRzedy

function [ R ] = LiczRzedy( X, Y )
%Y = [2,-1,-37,-58,464,-3382];
%X=[-4,-3,-1,0,2,4];
n = length(X);
R = zeros(n);
R(:,1) = Y;
%(-18) - (-3) = -15
%-1 - (-4) = 3
for i = 1:n-1
    for j = 1:n-i
        R(j,i+1)=(R(j+1,i)-R(j,i))/(X(i+j)-X(j));
    end
end