《《数值分析》matlab实验报告-梯形、辛普森求积公式(共3页).doc》由会员分享,可在线阅读,更多相关《《数值分析》matlab实验报告-梯形、辛普森求积公式(共3页).doc(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
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、 分析与思考专心-专注-专业