?? romberg.m
字號:
function X=romberg(f1,a,b,n)
l=1;
if n>=4
f=inline(f1);
h=b-a;
m=1;
for i=1:n
s1=0;
for k=1:m-1
s1=f(a+h*k)+s1;
end
T(i)=h/2*(f(a)+f(b))+h*s1;
h=h/2;
m=m*2;
if i>=2
S(i-1)=4/3*T(i-1)-1/3*T(i);
end
if i>=3
C(i-2)=-1/15*S(i-2)+16/15*S(i-1);
end
if i>=4
R(i-3)=-1/63*C(i-3)+64/63*C(i-2);
fprintf('%d%13.8f\n',l,R(i-3));
l=l+1;
end
end
X=R(n-3);
else
disp('n必須為大于3的整數');
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -