The purpose of this paper is to use a kind of analytical method called new homotopy perturbation method (NHPM), to solve nonlinear system of second order boundary value problems. The NHPM yields solutions in convergent series form with easily computable terms, and in some cases, yields exact solutions in one iteration. This method can be applied directly to the second order boundary systems needless of converting to the first order initial systems. To illustrate the application of the method, numerical results are derived using the calculated components of the NHPM. Comparisons with exact solutions show the efficiency and accuracy of NHPM in solving nonlinear system of second order boundary value problems.