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)