DividedDifferences.m

clc
clear
x0=input('x0=');
y0=input('y0=');
x=input('x=');
n=length(x0);
for k=1:n
F(k,1)=y0(k);
end
for i=2:n
for j=2:i
F(i,j)=(F(i,j-1)-F(i-1,j-1))/(x0(i)-x0(i-j+1));
end
end
p=F(n,1);
for a=2:n
Q=1;
for b=n:-1:n-a+2
Q=Q*F(n,a)*(x-x0(b))
end
p=p+Q;
end
disp(p);