python实现统计模拟【第3章】【利用随机数求积分】

糖逗

#利用随机数求积分
#对pi进行估计
import numpy as np
import math

class comput:
    def __init__(self, n = 7000000):
        self.n = n
    
    def random_number2(self, seed = 1, m = pow(2, 31)-1, c = 19,a = pow(7, 5)):   
        res = []
        xi_1 = seed
        for i in range(self.n):
            xi = (a * xi_1 + c) % m
            res.append(xi / m)
            xi_1 = xi
        return np.array(res)
    
    def comput_pi(self):
        x = self.random_number2(seed = 23)
        y = self.random_number2(seed = 61)
        res_temp = pow((2*x-1), 2) + pow((2*y-1), 2)
        res = np.sum(res_temp <= 1) / self.n
        return res*4

res = comput().comput_pi()
print(res)

参考：
1.Sheldon M.Ross《统计模拟》