使用Python进行高斯分布的蒙特卡洛方法

Question

我正在尝试为高斯实现cdf近似的蒙特卡罗方法。但它不起作用。我正在尝试从x = 0近似到x = inf，但我认为1000足够好。它应该给我接近一个答案。

from scipy import random
from scipy.stats import norm
import numpy as np


def guassian_approx(mean, var):
    sd = var**0.5 
    n=10000;
    a = 0;
    b = 1000
    x = random.uniform(-1000, 1000, n)
    sm = norm(mean, sd).pdf(x)
    p = (b-a) * np.mean(sm)   
    return p