Lagrange_interpol(라그랑주 보간법)

function [inter,error]=lagrange_interpol(x, y, new_x)

% x,y: 주어진 좌표점

% new x_ : 보간이 필요한 x값

n=length(x);

sum=0;

for i=1:n

    temp=y(i); %y(i)값 temp로 저장

for j=1:n %y(i)*Li(x)계산

    if(i~=j)

        temp=temp*(new_x-x(j))/(x(i)-x(j));

    end

end

sum=sum+temp; %Pn(x)계산

end

inter=sum; %Pn(x)값 출력

error=abs(tan(1.15)-inter)