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

《《数值分析》matlab实验报告-梯形、辛普森求积公式.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、东莞理工学院《数值分析》实验报告实验名称：牛顿插值法系别：计算机学院专业：2013级信息与计算科学班级：1班姓名：学号：实验日期：1、实验内容用不同数值方法计算积分。2、算法说明梯形求积公式算法：将积分区间划分为n等份，步长分点为 。积分。辛普森求积公式算法：其中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

2、;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;1、运行结果真值I=-4/9=-0.4444次数I（梯形公式）I（辛普森公式）5

3、0-0.7332-0.1150100-0.2695-0.8956200-0.4891-0.21022、分析与思考