如何计算π和e的值#

从小学就知道 $\pi=3.14159265358979\cdots$

但是一直不知道如何计算这个值。之前听说过祖冲之使用的是割圆法，但是我不会(T_T，流下了没有技术的泪水)，也就没办法自己亲自体验一下这个计算过程，好在后来我发现了别的方法可以计算出圆周率的数值

• 数学法

后来偶尔看到数学上有一种方法

$\pi=4\cdot(1-\cfrac{1}{3}+\cfrac{1}{5}-\cfrac{1}{7}+\cfrac{1}{9}\ldots+\cfrac{1}{n})$

使用这种方法只要n足够大, 精度就能提高, 所以我就想用程序来尝试计算一下pi的值, 毕竟计算机比人做这种事情合适的多

• 蒙特卡洛模拟

后来又想到可以采用蒙特卡洛采样的方法也可以计算出来 π

计算方法如下, 我们随机生成 m 个 二维坐标<x, y>, 并且满足 $0\leq x < 1, 0 \leq y < 1$。然后使用欧拉距离统计在圆内的坐标的数目 n ， 再除以生成的所有坐标的数目 $\cfrac{n}{m}$。由于我们的采样范围只覆盖了$\cfrac{1}{4}$的圆， 所以结果需要乘上4 就是圆周率的结果了，我们还可以去除随机出来的重复数字或者通过人为限定随机精度来加快计算

如何计算自然常数#

自然常数也一直是我不理解的一个数字

之前听说过一个类比,讲的挺好的,但当时并没有去仔细考虑过计算自然常数的方法数学常数e的含义

1
import math
2

3
class Calculator:
4

5
    num = 10000000
6
    a = 1
7
    res = 0
8

9
    def getPi(self):
10
        for n in range(1, self.num):
11
            m = 1 / (1 + 2 * n)
12
            if n % 2 == 0:
13
                self.a = self.a + m
14
            else:
15
                self.a = self.a - m
16
        print(self.a * 4)
17

18
    def getE(self):
19
        for n in range(0, self.num):
20
            self.res = self.res + 1 / math.factorial(n)
21
        print(self.res)
22

23

24
c = Calculator()
25
c.getPi()
26
c.getE()