计算机模拟在数学实验报告,数学实验报告2.docx

文档介绍：

数学实验报告2

篇一:数学实验报告 2

数学实验报告

实验序号: 日期:XX年 11 月 30 日

篇二:数学实验报告2

实验目的:

练****的求解方法。

实验内容:

1、用反正切函数的幂级数展开式结合有关公式求?,若要精确到以40位、50位数字,

试比较简单公式和Machin公式所用的项数。

2、用数值积分计算?,分别给出用梯形法和Simpson法精确到10位数字、用Simpson

法精确到15位数字时所用的项数n及的?近似值

3、用计算机模拟Buffon实验,给出n=1,000、10,000、1,000,000时的模拟结果。

实验要求:

撰写实验报告

写出试验过程中所使用的Mathematica程序或语句和计算结果

第一题:

Mathematica程序如下:

1. Arctan幂级数展开法:

In[1]:=k=100000;

S1=N[4*Sum[(-1)^(n-1)/(2n-1),{n,1,k}],40]

Out[2]= 3.141582653589793488462643352029502893728

In[3]:=k=100000;

S1=N[4*Sum[(-1)^(n-1)/(2n-1),{n,1,k}],50]

Out[4]= 3.1415826535897934884626433520295028937284193939649

2. 简单公式:

In[5]:=k=62;

S2=??[4?Sum[(?1)^(???1)?(1 2)^(2???1) (2???1)+(?1)^(???

1)?(1 3)^(2???1) (2???1) ,{??,1,??}],40]

Out[6]= 3.141592653589793238462643383279502884197

In[7]:=k=79;

S2=??[4?Sum[(?1)^(???1)?(1 2)^(2???1) (2???1)+(?1)^(???

1)?(1 3)^(2???1) (2???1) ,{??,1,??}],50]

Out[8]= 3.1415926535897932384626433832795028841971693993751

3. Machin公式

In[9]:=k=27;

S3=N[4*Sum[(-1)^(n-1)*(1/2)^(2n-1)/(2n-1)+(-1)^(n-1)*(1/3)^(2n-1)/(2n-1),{n,1,k}],40]

Out[10]= 3.141592653589793238462643383279502884197

In[11]:=k=35;

S3=N[4*Sum[(-1)^(n-1)*(1/2)^(2n-1)/(2n-1)+(-1)^(n-1)* (1/3)^(2n-1)/(2n-1),{n,1,k}],50]

Out[12]= 3.1415926535897932384626433832795028841971693993751

第二题:

1、梯形法:

In[13]:=y[x_]:=4/(1+x^2);

n=3000;

S3=N[1/(2*n)*(Su

内容来自淘豆网www.taodocs.com转载请标明出处.