Computational physics

A Monte Carlo simulation generates spin configurations of a magnetic lattice according to its thermal probability distribution. How can it estimate the mean magnetization at that temperature?