Solving systems of nonlinear equations is a relatively complicated problem for which a number of different approaches have been proposed. In this paper, we employ the Homotopy Analysis Method (HAM) to derive a family of iterative methods for solving systems of nonlinear algebraic equations. Our approach yields second and third order iterative methods which are more efficient than their classical counterparts such as Newton's, Chebychev's and Halley's methods. Keywords Homotopy analysis method