matlab 龙格现象,利用MATLAB分析数值积分中的龙格（Runge）现象（1）

实验目的： 观察Lagrange插值及数值积分中的龙格(Runge)现象。了解数值不稳定现象。

实验题目：(1)对于函数f(x)=1/(1+x^2),-4<=x<=4进行Lagrange插值。取不同结点数n，在区间[-4,4]上取等距间隔的结点为插值点，把f(x)和插值多项式的曲线画在同一张图上进行比较。(2)利用复化矩形或梯形公式计算f(x)=1/(1+x^2),-4<=x<=4在定义区间上的数值积分值，说明Runge现象。

(先完成题目(1)，题目(2)下次续上。)

在同一目录下编制一下的函数：f.m,langrange.m,runge1.m。

下面贴出每个文件的内容清单。

f.m:

function f= f( x )

f=1./(1+x.^2);

end

langrange.m：

function langrange= langrange( x,n )

langrange=0;

xx=linspace(-4,4,n+1);

for i=1:n+1

lix=1;

for

j=1:n+1

if j~=i

lix=lix.*((x-xx(j))./(xx(i)-xx(j)));

end

end

langrange=f(xx(i)).*lix+langrange;

end

end

runge1.m：

function r