Using the Simulated Annealing method with 400 fixed iterations, we got estimates for the model parameters as 227 kgm-3 for the density contrast and 2.47 km as the source depth of the Covington pluton (fig.7). We reached the estimates after four iterations with a minimum error of 1 X 10-4 (fig. 8a), the distribution for the predicted and observed gravity anomaly had two different distribution patterns (fig.8b).