数值分析龙贝格定积分程序.docVIP

编写程序用Romberg方法计算定积分 ，使其误差≤10-4。 function[quad,R]=Romberg(f,a,b,eps)h=b-a;R(1,1)=h*(feval(f,a)+feval(f,b))/2;M=1;J=0;err=1;while erreps J=J+1; h=h/2; S=0; for p=1:M x=a+h*(2*p-1); S=S+feval(f,x); end R(J+1,1)=R(J,1)/2+h*S; M=2*M; for k=1:J R(J+1,k+1)=R(J+1,k)+(R(J+1,k)-R(J,k))/(4^k-1); end err=abs(R(J+1,J)-R(J+1,J+1));endquad=R(J+1,J+1);