用龙贝格算法计算∫（0到1）[x/（4+x²）]dx的近似值.用龙贝格算法计算定积分0-1x/(4+x²）的近似值.在线求大神!

Question

用龙贝格算法计算∫（0到1）[x/（4+x²）]dx的近似值. 用龙贝格算法计算定积分0-1 x/(4+x²） 的近似值.在线求大神!

Answer 1

% 自定义函数
function y = f(x)
y = x/(4+x^2);
end
% 龙贝格算法
function result = romberg(f,a,b,e)
h = b-a;
k = 0;
T(1,1)=(b-a)*(f(a)+f(b))/2;
while true
k = k+1;
sum = 0;
for i=1:2:2^k
sum = sum + f(a+i*h/2^k);
end
T(k+1,1) = T(k,1)/2 + (h/2^k)*sum;
for i=1:k
T(k+1,i+1) = ((4^i)*T(k+1,i)-T(k,i))/(4^i-1);
end
if abs(T(k+1,1)-T(k,1))>