求高手看一下matlab编的龙贝格求积公式function [I,step]=Roberg(f,a,b,eps)M=1;h=b-a;err=1J=0;R=zeros(4,4);R(1,1)=h*(feval('f',a)+feval('f',b))/2delta)&(J
网友回答
试试我为你编的这个代码
function I=romberg(f,a,b,e)
%%对给定函数f(x),积分区间[a,b]和误差界e,用Romberg积分法计算积分
T1=0;T2=0;S1=0;S2=0;C1=0;C2=0;R1=0;R2=0;
h=b-a;
T1=h/2*(f(a)+f(b));
ee=1;while ee>efor k=1:2
R1=R2;
for j=1:2
C1=C2;
for i=1:2
S1=S2;
s=0;x=a+h/2;
while x