1、实验报告 -计算科学实验室实验名称： 圆周率的计算姓名 魏萱 学号 11112228 实验地点 T5-207实验类型 综合 设计 实验要求 选修 学时量 6所用知识 数学建模 数学软件 运筹学第一部分 割圆术割圆方式 周长 周长 周长 周长 周长分割程度 6 12 24 48 192所用工具 MATLAB MATLAB MATLAB MATLAB MATLAB圆周率值 3.000000 3.105829 3.132629 3.139350 3.141452matlab 代码function pai=geyuan(count);format long gR=1;d=1/2;h=sqrt(3)/2
2、;n=6;s=0;for i=1:count;s=d*h;edge=n;l=2*d*n;area=s*n;n=n*2;d=sqrt(d2+(R-h)2)/2;h=sqrt(R2-d2);endedgepai=l/(2*R);area;）第二部分 数值积分所用函数 21/x21/x2/x21/x21/x积分方法 矩形公式 矩形方式 矩形方式 矩形方式 矩形方式分割程度 100 1000 10000 100000 100000所用工具 VC+ VC+ VC+ VC+ VC+圆周率值 3.131576 3.140592 3.141494 3.141523 3.141348第三部分 蒙特卡洛方法实验报
3、告 -计算科学实验室所用函数 循环次数 100000 200000 300000 400000 500000满足的次数 79037 158411 237748 315925 395497所用工具 VC+ VC+ VC+ VC+ VC+圆周率值 3.141680 3.168220 3.169973 3.159250 3.163976第四部分 级数方法所用级数计算项数 10000 精确位数 3所用工具 MATLAB圆周率值 3.14149265第五部分 蒙特卡洛方法的应用计算积分 （填写所计算 的积分）积分结果 （填写该积分的准确结果）循环次数 （填写总的循 环次数）满足的次数 （填写满足条 件次数）所用工具 （填使用的软 件）积分的估计值（用蒙特卡洛方法估计的值）（ ， ， ，1,0sinyxdy1,02sinyxdy1,02sinyxd等供选用参考）1,02lyx第六部分 心得体会（请写出您对这个实验过程的意见，建议及想法）第一部分 matlab 代码function y=calpi(n)syms a;for i=1:na=sqrt(2-sqrt(4-a2);enda=subs(a,a,1);实验报告 -计算科学实验室y=3*2n*vpa(a,n+5);
