sampling

This is a Metropolis-Hastings algorithm that could be used to sample simple posterior surfaces. It can be tuned to run with a search phase, employing a simulated annealing scheme for more efficient search of the parameter space.