In the past three decades numerous multi-objective evolution- ary algorithms have been developed and tested as trustable and efficient solution methods to solve multi-objective models (Deb, 2001). However, these algorithms are best known to their cap- ability of obtaining good or near optimal solutions and attainment of the exact optimal solution(s) is never guaranteed. In order to solve the multi-objective model (16), we consider goal programming method as a subroutine of the solution approach. The major drawback of the standard goal programming method is that the method can obtain only one non-dominated solution which is highly dependent to the decision maker's choice of the goals and the weights of deviation from the predefined goals. To rectify this dependability and in order to obtain the true Pareto- optimal front, the following hybrid Monte Carlo simulation model is proposed in which randomly generated objective goals and deviation weights are used in the goal programming submodel in each simulation replication.