实验名称 圆周率计算

《实验名称 圆周率计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验报告----计算科学实验室实验名称：圆周率的计算姓名魏萱学号11112228实验地点T5-207实验类型综合设计实验要求选修学时量6所用知识数学建模数学软件运筹学第一部分割圆术割圆方式周长周长周长周长周长分割程度6122448192所用工具MATLABMATLABMATLABMATLABMATLAB圆周率值3.0000003.1058293.1326293.1393503.141452matlab代码functionpai=geyuan(count);formatlonggR=1;d=1/2;h=sqrt(3)/2;n=6;s=0;fori=1:co

2、unt;s=d*h;edge=n;l=2*d*n;area=s*n;n=n*2;d=sqrt(d^2+(R-h)^2)/2;h=sqrt(R^2-d^2);endedgepai=[l/(2*R);area];）第二部分数值积分所用函数积分方法矩形公式矩形方式矩形方式矩形方式矩形方式分割程度100100010000100000100000所用工具VC++VC++VC++VC++VC++圆周率值3.1315763.1405923.1414943.1415233.141348第三部分蒙特卡洛方法实验报告----计算科学实验室所用函数循环次数1000002000

3、00300000400000500000满足的次数79037158411237748315925395497所用工具VC++VC++VC++VC++VC++圆周率值3.1416803.1682203.1699733.1592503.163976第四部分级数方法所用级数计算项数10000精确位数3所用工具MATLAB圆周率值3.14149265第五部分蒙特卡洛方法的应用计算积分（填写所计算的积分）积分结果（填写该积分的准确结果）循环次数（填写总的循环次数）满足的次数（填写满足条件次数）所用工具（填使用的软件）积分的估计值（用蒙特卡洛方法估计的值）（，，，等

4、供选用参考）第六部分心得体会（请写出您对这个实验过程的意见，建议及想法）第一部分matlab代码functiony=calpi(n)symsa;fori=1:na=sqrt(2-sqrt(4-a^2));enda=subs(a,’a’,’1’);y=3*2^n*vpa(a,n+5);实验报告----计算科学实验室