数值分析实验报告(5)

err=10e-5; N=6;

t=zeros(N+1);

t(1,1)=(b-a)*(fa+fb)/2; j=1;

for i=j+1:N n=2^(i-2);

h=(b-a)*2^(2-i); for k=1:n

x(k)=a+(2*k-1)*h/2; end

y=4./(1+(x).^2); H=h*sum(y);

t(i,j)=(t(i-1,j)+H)/2; end

for n=j:N

for m=n:N-1

t(m+1,n+1)=((4^(n))*t(m+1,n)-t(m,n))/(4^(n)-1); end

if abs(t(n+1,n+1)-t(n,n))

disp(t(n,n));

七、个人实验总结

三者比较龙贝格的效果是最好的。

２０