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

《《数值分析》matlab实验报告-梯形、辛普森求积公式(共3页).doc》由会员分享，可在线阅读，更多相关《《数值分析》matlab实验报告-梯形、辛普森求积公式(共3页).doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上东莞理工学院数值分析实验报告实验名称：牛顿插值法系 别：计算机学院专 业：2013级信息与计算科学班 级：1班姓 名：学 号：实验日期：1、 实验内容用不同数值方法计算积分。2、 算法说明梯形求积公式算法：将积分区间划分为n等份，步长分点为。积分。辛普森求积公式算法：其中h为步长。3、 Matlab软件程序清单梯形求积公式TiXing_quad(a,b,h)：function t = TiXing_quad(a,b,h)%a为积分下界，b为积分上界，h为步长。format longx = a:h:b;y = sqrt(x).*log(x);y(1) = 0;t = 0

2、;for k=1:(b-a)/h, t=t+y(k)+y(k+1);endt=t*h/2;辛普森求积公式Sinpson_quad(a,b,h)：function s=Sinpson_quad(a,b,h)format longx=a:h:b;y=sqrt(x).*log(x);z=sqrt(x+h/2).*log(x+h/2);y(1)=0;s=0;for k=1:(b-a)/h, s=s+y(k)+y(k+1)+4*z(k);ends=s*h./6;4、 运行结果真值I=-4/9=-0.4444次数I（梯形公式）I（辛普森公式）50-0.7332 -0.1150100-0.2695 -0.8956200-0.4891 -0.21025、 分析与思考专心-专注-专业