更新**
我正在努力返回预期值 E[x^2]。我已经创建了正常函数来对平均值 mu 和 sigma2 的正常值进行采样。我想使用这个函数生成正常值，然后将fun应用于每个值。我是否需要在循环中定义期望值，然后应用fun？我可以在为每个 x 找到 x^2 后取平均值吗？这是我一直在做的事情：
def fun(x):
y=x**2
return y
def EV(fun,n_samples,mu,sigma2):
for i in range(1,n_samples):
... u= np.random.uniform(0,1,1) ... output.append(float(y)) ... 返回输出