这是此处的代码示例:
函数的参数赋值中的相等是什么意思?喜欢N=20000
这里?N
那和简单的论点有什么区别?导入随机数,数学
def gibbs(N=20000,thin=500):
x=0
y=0
samples = []
for i in range(N):
for j in range(thin):
x=random.gammavariate(3,1.0/(y*y+4))
y=random.gauss(1.0/(x+1),1.0/math.sqrt(x+1))
samples.append((x,y))
return samples
smp = gibbs()