In this paper, we present two new methods for finding the numerical solutions of systems of the nonlinear equations. The basic idea depend on founding relationship between minimum of a function and the solution of systems of the nonlinear equations. The first method seeks the numerical solution with a sequence of search directions, which is depended on gradient and Hessian matrix of function, while the second method is based on a sequence of conjugate search directions. The study shows that our two methods are convergent, and they can find exact solutions for quadratic functions, so they can find high accurate solutions for over quadratic functions. The purposed two algorithms are programmed by Mathematica Version9. The approximate solutions of some test problems are given. Comparisons of our results with other methods illustrate the efficiency and highly accurate of our suggested methods.