编程MATLAB程序实现复化梯形和辛普森数值积分.docxVIP

-- -- 数值分析实验报告—— [1] 掌握复化梯形和辛普森数值积分法的基本原理和方法； 验 [2] 编程 MATLAB 程序实现复化梯形和辛普森数值积分 的 编程序实现复化梯形数值积分求积公式 function y=f(x) y=sqrt(x).*log(x); 实 验 function T_n=F_H_T(a,b,n) 内 h=(b-a)/n; 容 for k=0:n 与 x(k+1)=a+k*h; 步 if x(k+1)==0 骤 x(k+1)=10^(-10); end end T_1=h/2*(f(x(1))+f(x(n+1))); for i=2:n F(i)=h*f(x(i)); end T_2=sum(F); T_n=T_1+T_2; 实 验 内  运行结果： T_n=F_H_T(0,1,20) T_n = -0.4336 2.编程序实现复化辛普森数值积分求积公式 与 步 function y=f(x) y=sqrt(x).*log(x); function S_n=S_P_S(a,b,n) h=(b-a)/n; for k=0:n x(k+1)=a+k*h; x_k(k+1)=x(k+1)+1/2*h; if (x(k+1)==0)|(x_k(k+1)==0) x(k+1)=10^(-10); x_k(k+1)=10^(-10); end S_1=h/6