I'm a beginner at python and am learning to use MCMC sampling methods, using python's emcee package. As a beginner exercise I want to sample a Maxwell-Boltzmann Distribution. I have an example code which samples a Gaussian, defined through the function below;
def lnprob(x, mu, icov):
diff = x-mu
return -np.dot(diff,np.dot(icov,diff))/2.0
How would I express a boltzmann distribution in the same way?