A steady two-dimensional laminar boundary layer flow of an incompressible viscous fluid over a semi-infinite surface is considered to investigate the accuracy of the homotopy analysis method. The governing coupled nonlinear system of differential equations is solved by means of the HAM approach. Explicit analytical series solutions are obtained and compared with numerical solutions. Good agreement is observed between the numerical results and HAM analytical solutions.