In this study, the Legendre operational matrix method based on collocation point is introduced to solve high order ordinary differential equations with some nonlinear terms arising in physics and mechanics. This technique transforms the nonlinear differential equation via mixed conditions into a matrix equation with unknown Legendre coefficients. This solution of this matrix equation yields the Legendre coefficients of the solution function. Thus, the approximate solution is obtained in terms of Legendre polynomials. Some test problems together with residual error estimation are given to show the usefulness and applicability of the method and the numerical results are compared.