lagrange.m

% 拉格朗日插值
% 注释的地方需要修改
clear

x=[0.5,0.6,0.7,0.8,0.9,1.0];                      % 插值节点，逗号分隔
y=[-0.693147,-0.510826,-0.356675,-0.233144,-0.105360,0];          % 插值节点对应的函数值，逗号分隔
x0=0.75;                               % 求值节点

n=length(x);
y0=0;
for k=1:n 
    T=1;   
    for i=1:n
        if i~=k
            T=T*(x0-x(i))/(x(k)-x(i));
        end 
    end
    y0=y0+T*y(k); 
end 

vpa(y0,6)                                % 结果有效数字位数