We discuss the geophysical inversion methodology for volcanic muography based on the Simulated Annealing algorithm, using a semi-empirical model of the muon flux reaching the volcano, its surrounding topography and a framework for the energy loss of muons in rock. We determined the minimum muon energy-as function of the arrival direction-needed to cross the volcanic building, the emerging integrated flux of muons and the density profile inside a model of Cerro Machín volcano (Tolima, Colombia) within a maximum error of 1% concerning the theoretical model.