www.gusucode.com > MGS求解方程组的MATLAB程序 > code6/lagelangrichazhi.m

    clc;clear;close;
n=input('输入 n:');
x=input('输入矩阵x');
y=input('输入矩阵y');
L=zeros(100,1);
syms t
for k=1:100
for i=1:n
    mul=1;
    for j=1:n
        if(i~=j)
            x0=0+k*0.004;
            mul=mul*(x0-x(j))/(x(i)-x(j));
        end
    end
    L(k)=L(k)+y(i)*mul;
end
end
plot(L)
hold on
L1=zeros(100,1);
for k=1:100
    x0=0+k*0.004;
    L1(k)=1/sqrt(2*pi)*int(exp(-t^2/2),t,-1000,x0);
end
plot(L1)