1、龙贝格积分-matlab通用程序__________________________________________________数值分析课程的大作业,教材《数值分析》李乃成.梅立泉clear
clc
formatlong
f=input('请输入原函数f=','s');
a=input('积分下限a=');
b=input('积分上限b=');
eps1=input('精度eps1=');
T(1)=double((b-a)/2*(limit(sym(f),findsym(sym(f)),a)+limit(sym(f),findsym(sym(f)),b)));
f
2、ork=2:4
sum1=0;
fori=1:2^(k-2)
sum1=sum1+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k-1)));
end
T(k)=1/2*T(k-1)+(b-a)/(2^(k-1))*sum1;
end
fork=1:3
____________________________________________________________________________________________________S(k)=T(k+1)+1/(4-1)*(T(
3、k+1)-T(k));
end
fork=1:2
C(k)=S(k+1)+1/(4^2-1)*(S(k+1)-S(k));
end
R(1)=C(2)+1/(4^3-1)*(C(2)-C(1));
k=3;
while1
T(1)=T(2);
T(2)=T(3);
T(3)=T(4);
sum2=0;
fori=1:2^k
sum2=sum2+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k+1)));
end
T(4)=1/2*T(4)+(b-a)/2^(k+1)*sum2;