In this paper, we have employed the variational iteration method to solve the boundary value problems. Numerical results reveal that it is a very effective method compared with the results obtained by using the Adomian decomposition method in Wazwaz, A. M. (2000).