1、计算方 法上 机报 告 27 附录 4 龙贝 格积分 法的 matlab 程序 clc;clear; s = input( 请输入被 积函 数 表达式 ：nf(x) = ,s); a = input(n 请输入积 分区 间下限 ：n); b = input(n 请输入积 分区 间上限 ：n); e = input(n 请输入积 分误 差：n); N = -log10(e); f = inline(s); i = 1; error = 1; syms x f_a = limit(s,x,a); f_b = limit(s,x,b); T(1) = (b-a)/2*(double(f_a)+dou
2、ble(f_b); S(1) = 0; C(1:2) = 0; R(1:3) = 0; while errore | i=3 C(i) = S(i)+1/15*(S(i)-S(i-1); if i=4 R(i) = C(i)+1/63*(C(i)-C(i-1); error = abs(R(i)-R(i-1); end end end for k=1:i T_int = fix(T(k); T_dec = round(T(k)-T_int)*10N); S_int = fix(S(k); S_dec = round(S(k)-S_int)*10N); C_int = fix(C(k); C_d
3、ec = round(C(k)-C_int)*10N); R_int = fix(R(k); R_dec = round(R(k)-R_int)*10N); if k=1 fprintf(n%0.f.%0.fn,T_int,T_dec); elseif k=2 fprintf(%0.f.%0.f %0.f.%0.fn,T_int,T_dec,S_int,S_dec); elseif k=3 fprintf(%0.f.%0.f %0.f.%0.f %0.f.%0.fn,T_int,T_dec,S_int,S_dec,C_int,C_dec); else fprintf(%0.f.%0.f %0.f.%0.f %0.f.%0.f %0.f.%0.fn,T_int,T_dec,S_int,S_dec,C_i nt,C_dec,R_int,R_dec); end end fprintf(n 所求积 分值 为 %0.f.%0.f, 误 差为%.e 。n,R_int,R_dec,e);
