《数值分析》matlab实验报告.docx

文档介绍：东莞理工学院
《数值分析》
实验报告
实验名称：牛顿插值法
系 别：计算机学院
专 业：2013级信息与计算科学
班 级：1班
姓名：
学 号：
实验日期：
1、实验内容
用不同数值方法计算积分J \ x
In 东莞理工学院
《数值分析》
实验报告
实验名称：牛顿插值法
系 别：计算机学院
专 业：2013级信息与计算科学
班 级：1班
姓名：
学 号：
实验日期：
1、实验内容
用不同数值方法计算积分J \ x
In xdx =-
2、算法说明
梯形求积公式算法：
将积分区间［a,b］划分为n等份，步长匕=h分点
n
为 x = a + kh, k = 1,2,..., n。
k
积分I =
辛普森求积公式算法:
鳥f（x）dx=3（y+4 犷打-9of（4）（c）
其中h为步长。
3、Matlab软件程序清单
梯形求积公式 TiXing_quad(a,b,h)：
function t = TiXing_quad(a,b,h)
%a为积分下界，b为积分上界，h为步长。
format long
x = a:h:b;
y = sqrt(x).*log(x);
y(1) = 0;
t = 0;
for k=1:(b-a)/h,
t=t+y(k)+y(k+1);
end
t=t*h/2;
辛普森求积公式Sinpson quad(a,b,h):
function s=Sinpson_quad(a,b,h) format long
x=a:h:b;
y=sqrt(x).*log(x);
z=sqr