romberg.m

来自「various matlab programs to slove various」· M 代码 · 共 5 行

function [A]=romberg(a,b,n,fun);
for i=1:(n+1), A(i,1)=trapezc(a,b,2^(i-1),fun); end;
for j=2:(n+1), for i=j:(n+1),  
A(i,j)=(4^(j-1)*A(i,j-1)-A(i-1,j-1))/(4^(j-1)-1); end; end;