matlab中 迭代 源代码

function T=romberg(a,b,n)
h=b-a;
f=inline('x');
T=zeros(2,n);
T(1,1)=(f(a)+f(b))*h/2;
for i=2:n
s=0.0;
for j=1:(2^(i-2));
s=s+f(a+(j-0.5)*h);
end
T(2,1)=T(1,1)/2+s*h/2;
for k=2:i
T(2,k)=((4^(k-1))*T(2,k-1)-T(1,k-1))/((4^(k-1))-1);
end
h=h/2;
for j=1:i
T(1,j)=T(2,j);
end
end
这是romberg迭代法的源代码