traprl.m

数值分析常用几种方法

function s=traprl(a,b)
%区间逐次分半的复化梯形计算积分
m=1;
z=0;
T=(b-a)/2*(exp(-a^2)+exp(-b^2));
while 1
    t=T;
h=(b-a)/2^m;
for i=1:2^(m-1)
    c(i)=a+(2*i-1)*h;
    d(i)=exp(-c(i)^2);
end
for j=1:2^(m-1)
    z=z+d(j);
end
    T=0.5*T+h*z;   
if abs((T-t)/T)<10e-4                    
  break;
else
    m=m+1;
    z=0;
end
end
m                              
T