使用Python的主持人对Boltzmann分布进行采样

Question

我是python的初学者，正在学习使用MCMC采样方法，使用python的emcee包。 作为一名初学者练习，我想对Maxwell-Boltzmann分布进行抽样。 我有一个示例代码，它采样高斯，通过下面的函数定义;

    def lnprob(x, mu, icov):
diff = x-mu
return -np.dot(diff,np.dot(icov,diff))/2.0

我如何用同样的方式表达一个boltzmann分布？