The homotopy analysis method (HAM) is applied to solve the nonlinear Schrödinger (NLS)
equations. In this paper, we will reduce the NLS equation to a system of two nonlinear
equations contain the real and imaginary parts of the solution. The method provides the
solution in the form of a rapidly convergent series with easily computable components using
symbolic computation software such as Mathematica. The scheme shows importance of
choice of convergence-control parameter ħ to guarantee the convergence of the solutions of
nonlinear differential equations. This scheme is tested on two cases study, the cubic nonlinear
Schrödinger (CNLS) equation and a system of coupled nonlinear Schrödinger equations. The
results demonstrate reliability and efficiency of the algorithm developed.