tablade diferencia divididas

%tablade diferencia divididas
x=[-1;0;1;2];
f=[-2,1,3,13];
n=length(x)-1;
for i=1:n
t(i,1)=(f(i+1)-f(i));
end
for j=2:n
for i=j:n
t(i,j)=((t(i,j-1)-t(i-1,j-1))/(x(i+1)-x(i-j+1)));
end
end
t