function [y,R]=lagranzi(X,Y,x,M)          %構(gòu)造拉格朗日函數(shù)
%X為節(jié)點(diǎn)的橫坐標(biāo)向量;
%Y為節(jié)點(diǎn)的縱坐標(biāo)向量;
%x為最終所要求的節(jié)點(diǎn);
n=length(X);m=length(x);                 %n為節(jié)點(diǎn)X的長(zhǎng)度,m為插值點(diǎn)x的長(zhǎng)度
for i=1:m                                %i取從1到m的值,即對(duì)i進(jìn)行循環(huán)
    z=x(i);s=0.0;                        %令z等于每一個(gè)插值點(diǎn),s取初值為0.0
    for k=1:n                            %k取從1到n的值,即對(duì)k進(jìn)行循環(huán)
        p=1.0;q1=1.0;c1=1.0;             %p取初值為1.0,q1取值為1.0,c1取初值為1.0
        for j=1:n                        %對(duì)j進(jìn)行循環(huán),從1到n進(jìn)行循環(huán)
            if j~=k                      %當(dāng)j不等于k時(shí)
                p=p*(z-X(j))/(X(k)-X(j));%得出拉格朗日的式子的系數(shù)
            end
            q1=abs(q1*(z-X(j)));c1=c1*j;
        end
        s=p*Y(k)+s;                      %得出拉格朗日的插值函數(shù)
    end
    y(i)=s;
end
R=M*q1/c1;                               %R為插值余項(xiàng)