?? traprl.m
字號(hào):
function s=traprl(a,b)
%區(qū)間逐次分半的復(fù)化梯形計(jì)算積分
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
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -