龙贝格积分-matlab通用程序(共2页).docx

精选优质文档-倾情为你奉上数值分析课程的大作业，教材数值分析李乃成.梅立泉clearclcformat longf=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);for k=2:4 sum1=0; for i=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;endfor k=1:3 S(k)=T(k+1)+1/(4-1)*(T(k+1)-T(k);endfor k=1:2 C(k)=S(k+1)+1/(42-1)