In this paper, both Least Squares minimization (LSM) and Variational Monte-Carlo (VMC) techniques have been implemented to determine the co-efficients of Semi-Empirical Mass Formula (SEMF). First, the experimental binding energies (BEs) are determined for all the available nuclei from Atomic Mass Evaluation (AME2016) data. Then, LSM technique is implemented in Gnumeric worksheet to obtain the SEMF co-efficients by considering only the first three co-efficients which are deduced from Liquid Drop Model (LDM). The mean squared error (MSE) value, between obtained BEs from the optimized co-efficients and the experimental BEs, has been determined. Then, to emphasize the relevance of empirical terms, they have been introduced successively one after other and the procedure is repeated. A reduction in MSE-value has been observed after each iteration. This same procedure has also been employed using Monte-Carlo approach to obtain SEMF co-efficients by minimizing MSE-value as in variational principle. A comparative analysis has shown that the optimized parameters using VMC have resulted in smaller MSE-value than that of LSM technique.