maxp_sa function

A simulated annealing algorithm to solve the max-p-region problem