首先可以了解一下Monte Carlo 方法,这里是链接https://en.wikipedia.org/wiki/Monte_Carlo_method。

根据圆的面积公式S=πR^2,当R=1时,S=π。
我们在一个边长为1的正方形中划定一个半径为1的1/4圆区域。我们向这个正方形区域内随机投点,根据几何概型(高中数学知识,我在此不再赘述),点落到1/4圆区域内的概率为P=π/4. 这个概率的四倍便可以视为π的近似值。
根据这个思路,我们要想用程序模拟出π的近似值,就可以转换为模拟这样一个几何概型求它的概率的问题。以下是用c语言实现这个过程的代码:

/*这个精度有点看概率,一开始我宏定义N为2000,得到了3.16,
改成10000得到了3.19,
改成20000得到了3.17,
改成200000,得到了3.14,
之后又用200000测试了4次,得到了三次3.14,一次3.13*/
改成10000得到了3.19,
改成20000得到了3.17,
改成200000,得到了3.14,
之后又用200000测试了4次,得到了三次3.14,一次3.13*/