《数值分析》matlab实验报告-梯形、辛普森求积公式

东莞理工学院《数值分析》实验报告实验名称：牛顿插值法系别：计算机学院专业：2013级信息与计算科学班级：1班姓名：学号：实验日期：1、实验内容用不同数值方法计算积分104ln9xxdx。2、算法说明梯形求积公式算法：将积分区间[,]ab划分为n等份，步长bahn分点为,1,2,...,kxakhkn。积分10()2nkkkhIxx。辛普森求积公式算法：05(4)012()(4)()390xxhhfxdxyyyfc其中h为步长。3、Matlab软件程序清单梯形求积公式TiXing_quad(a,b,h)：functiont=TiXing_quad(a,b,h)%a为积分下界，b为积分上界，h为步长。formatlongx=a:h:b;y=sqrt(x).*log(x);y(1)=0;t=0;fork=1:(b-a)/h,t=t+y(k)+y(k+1);endt=t*h/2;辛普森求积公式Sinpson_quad(a,b,h)：functions=Sinpson_quad(a,b,h)formatlongx=a:h:b;y=sqrt(x).*log(x);z=sqrt(x+h/2).*log(x+h/2);y(1)=0;s=0;fork=1:(b-a)/h,s=s+y(k)+y(k+1)+4*z(k);ends=s*h./6;4、运行结果真值I=-4/9=-0.444444444444444次数I（梯形公式）I（辛普森公式）50-0.441090226387332-0.443793798301150100-0.443117905322695-0.444194510818956200-0.443925359444891-0.4443490454521025、分析与思考