mathematica数学实验报告 实验二.docxVIP

数 学 实 验 报 告 实验二 学院：数学与统计学院 班级：信息与计算科学(1)班 姓名： 郝玉霞 学号：201171020107 实验二 一、实验名称：π的计算 二、实验目的：首先在 Mathematica 环境中用多种方法计算圆周率? 的值，通过实验来体会各种方法的区别，比较各种方法的优劣，接着尝试自己提出新的方法来计算圆周率? 的值。 三、实验环境：学校机房，Mathematica 软件。 四、实验的基本理论和方法 1、用Mathematica绘图函数Plot绘制圆周率? ； 2、计算圆周率? 的数值积分法、泰勒级数法、蒙特卡罗法，并且利用特定的公式来计算圆周率? 。 五、实验的内容和步骤及实验的结果和结果分析步骤一、数值积分法计算? 因为单位圆的半径为 1，它的面积等于? ，所以只要计算出单位圆的面积， 就算出了? 。在坐标轴上画出以圆点为圆心，以 1 为半径的单位圆，则这个单位圆在第一象限的部分是一个扇形，而且面积是单位圆的 1/4，于是，我们只要算出此扇形的面积，便可以计算出? 。 当 n=5000 时； 语句： n=5000;y[x_]:=4/(1+x*x); s1=(Sum[y[k/n],{k,1,n-1}]+(y[0]+y[1])/2)/n; s2=(y[0]+y[1]+2*Sum[y[k/n],{k,1,n-1}]+4*Sum[y[(k-1/2)/n],{k,