数学软件作业圆周率的计算.docVIP

2007数学应用软件设计内容 ——怎样计算圆周率π的值 实验报告 姓名：。。 学号：；；；； 实验目的： 利用所学的数学应用软件，解决下列问题： 1. 用反正切函数的幂级数展开式结合有关公式求π，若要精确到以40位、50位数字，试比较简单公式和Machin公式所用的项数； 2. 用数值积分计算π，分别用梯形法和Simpson法精确到10位数字，用Simpson法精确到15位数字； 3. 用Monte Carlo 法计算π，除了加大随机数，在随机数一定时可重复算若干次后求平均值，看能否求得5位精确数字？ 4. 设计方案用计算机模拟Buffon实验； 5.利用学习过的知识(或查阅资料)，提出其他计算π的方法。 实验的基本理论与方法： 1.利用反正切函数的幂级数表达式： 利用积分公式两边积分得： 在上式中令x=1得： 用Mathematica计算In[1] k=1000; S1=N[4*Sum[(-1)^(n-1)/(2n-1),{n,1,k}],18] [Out2] 3.14059265383979293 In[3] k=10000; [Out4] 3.14149265359004324 In[5] k=15000; [Out6] 3.14152598692320065 In[7]