Abstract: This study has been devoted to investigate the semi-analytical solution of nonlinear differential equations with boundary value problems (BVPs). Modified variation iteration method has been utilized to solve some BVPs nonlinear differential equations. In this method, general lagrange multipliers have been introduced to construct correction functions for the problems. The multipliers can be identified optimally via the variational theory. The results have been compared with those of exact solutions and Homotopy Analysis Method (HAM). A clear conclusion can be drawn from the numerical results that proposed method provides excellent approximations to the solutions of this kind of problems in the terms of simplicity and accuracy, thus, it can be easily extended to other BVPs nonlinear differential equations and so can be found widely applicable in engineering sciences.