Aimed at the widespread significance of solving high dimensional nonlinear equations in practical engineering, for the faultiness of classical methods about initial value sensitivity, systems of nonlinear equations are transformed to function optimization, combining with the advantages of PSO, SA and LM, a new method is proposed, in which, according to the idea of SA, LM operator is embedded in PSO. Besides, adaptive damp factor is led into LM operator, which improves efficiency of convergence. The method takes full advantages of above three methods, solving the problems such as initial value sensitivity of LM and falling into local extreme value of PSO. Experimental results of comparing with other methods show that the proposed method has reliable convergence probability and high computation rate, which is an effective method of dealing with systems of nonlinear equations in practical engineering.