na-lab-lagrange

code

clc;

close all;

x=input("Enter values of xi");

y=input("Enter values of fi");

n=length(x);

p=input("Enter point");

for i=1:n

l(i)=1.0;

for j=1:n

if(j~=i)

l(i)=((p-x(j))/(x(i)-x(j)))*l(i);

end

end

end

sum=0.0;

for i=1:n

sum=sum+l(i)*y(i);

end

disp(sum);