Markov chain Monte Carlo method