计算π的值

圆周率π是一个无理数，没有任何一个精确公式能够计算π值，π的计算只能采用近似算法。国际公认采用蒙特卡洛方法计算。蒙特卡洛(Monte Carlo)方法，又称随机抽样或统计试验方法。当所求解问题是某种事件出现的概率，或某随机变量期望值时，可以通过某种“试验”的方法求解。简单说，蒙特卡洛是利用随机试验求解问题的方法。

首先构造一个单位正方形 和 1/4圆。随机向单位正方形和圆结构抛洒大量点，对于每个点，可能在圆内或者圆外，当随机抛点数量达到一定程度，圆内点将构成圆的面积，全部抛点将构成矩形面积。圆内点数除以圆外点数就是面积之比，即π/4。随机点数量越大，得到的π值越精确。

由于DARTS点数量较少，π的值不是很精确。通过增加DARTS数量继续试验，同时，运行时间也逐渐增加。

代码及执行结果

以上是Python语言编写的程序，运行较慢。采用Fortran语言编写程序，会快很多，以下是抛洒不同的点，程序运行时间比较。

蒙特卡洛方法提供了一个利用计算机中随机数和随机试验解决现实中无法通过公式求解问题的思路。它广泛应用在金融工程学，宏观经济学，计算物理学（如粒子输运计算、量子热力学计算、空气动力学计算）等领域。