In this article, we propose a high order method for solving steady and unsteady two-dimensional laminar boundary-layer equations. This method is convergent of sixth-order of accuracy. It is shown that this method is unconditionally stable. The unsteady separated stagnation point flow, the Falkner–Skan equation and Blasius equation are considered as special cases of these equations. Numerical experiments are given to illustrate our method and its convergence.