文档介绍：: .
y=;
z1=sum(y(1:200))*h
z2=sum(y(2:201)): .
y=;
z1=sum(y(1:200))*h
z2=sum(y(2:201))*hz=cumsum(y);
z11=z(200)*h;z12=(z(201)-z(1))*h;z3=trapz(x,y)
：y=inline('');z1=quad(y,0,1,100)z2=quadl(y,0,1,100)
z3=quad8(y,0,1,100,trace(10))
%simpleson公式
%复化simpleson公式%simpleson8阶公式法
学号
班级
统计1001
姓名
指导教师
易昆南
实验题目
用多种方法计算数值积分
评分
1、设计（实****目的：
1•了解MATLAB在实际问题中的应用2・通过实践加深对这门语言中M文件的了解
3・熟悉简单程序结构，如循环结构（for循环、while循环）选择结构（if-else-if）、分支语句（switch-case-otherwise））
2、实验内容：
（1）.分别用左、右矩形法，梯形法，复化辛普森公式计算y=x'2在［0,1］上的定积分；
（2）•用蒙特卡罗随机投点法计算y=1/（1+x'2）在［0,1］上的定积分，并求出pi的近似值；
（3）.用蒙特卡罗均值估计法计算y=x"2在［0,1］上的定积分。
3. 详细设计：
、右矩形法和梯形法:h=1/200;x=0:h:1;
%左矩形法%右矩形法
%等同z1
%等同z2三•蒙特卡罗随机投点法：n=100000;
k=0;
fori=1:n
x=rand;
y=rand;
ify<=1/(1+x入2);
k=k+1;
end
endz=k/npi=4*k/n
三•蒙特卡罗随机投点法：n=100000;
k=0;
fori=1:n
x=rand;
y=rand;
ify<=1/(1+x入2);
k=k+1;
end
endz=k/npi=4*k/n
%产生（0,1）区间的随机数
%对y=1/（1+x入2）面积投点
%由积分pi/4=k/n而来，前者是概率，后者是频率
%梯形法，等同于z3=trapz（y）*：n=10000;x=rand(1,n);y=;
z=sum(y)/n%均