metropSB function

Metropolis-Szymura-Barton algorithm