In this paper, an approach based on the variational iteration method (VIM) is proposed with an auxiliary parameter to predict the multiplicity of the solutions of homogeneous nonlinear ordinary differential equations with boundary conditions. The proposed approach is capable to predict and calculate all branches of the solutions simultaneously. Four practical problems are chosen to show the efficiency and importance of the proposed method. The proposed approach successfully detects multiple solutions to