高考问答csgo调fov:Romberg算法的m程序
来源：百度文库 编辑：高考问答 时间：2024/05/05 08:22:37
f为m文件定义的函数，从a到b积分，n为次数参数，求R(n,n)
最后结果以列表形式输出，同时将最后两行的结果存储在r中。
%Romberg Integration
function r = romberg(f,a,b,n)
h = b-a;
r(1,1) = h*(feval(f,a)+feval(f,b))/2
for i =2:n
sum = 0;
for k =1:2^(i-2)
sum = sum+feval(f,a+(k-0.5)*h);
end;
r(2,1)=(r(1,1)+h*sum)/2;
for j = 2:i
r(2,j)=r(2,j-1)+(r(2,j-1)-r(1,j-1))/(4^(j-1)-1);
end;
r(2,:)
h = h/2;
r(1,:)=r(2,:);
end;