Python练习之圆周率的计算

法一：拉马努金公式

import numpy
a=input("k=")
k=int(a)
print ("圆周率的值为")
i=0
S=0
for i in range (k):
    a=numpy.math.factorial(4*i)
    b=numpy.math.factorial(i)
    S+=a*(26390*i+1103)/(pow(b,4)*pow(396,4*i))
A=1/(2*(2**(1/2))/(99**2)*S)
print (A)

法二：蒙特卡洛方法

import random
N=input("投掷次数为：")
n=int(N)
i=0
k=0
for i in range(n):
    x=random.uniform(0,1)
    y=random.uniform(0,1)
    if x**2+y**2<1:#在四分之一圆内的点
        k=k+1
print ("pi=",4*k/n)