na-lab-newtondivided differences

code

clc;

x=input("enter x");

y=input("enter y");

n=length(x);

f=zeros(n,n);

f(:,1)=y;

p=input("enter point");

for i=2:n

for j=i:n

f(j,i)=(f(j,i-1)-f(j-1,i-1))/(x(j)-x(j-i+1));

end

end

for i=1:n

prod(i)=1;

for j=1:i-1

prod(i)=prod(i)*(p-x(j));

end

end

sum=0;

for i=1:n

sum=sum+(f(i,i)*prod(i));

end

disp(sum);