Abstract
Optimal homotopy asymptotic method (OHAM) is proposed to solve linear and nonlinear systems of second-order boundary value problems. OHAM yields exact solutions in just single iteration depending upon the choice of selecting some part of or complete forcing function. Otherwise, it delivers numerical solutions in excellent agreement with exact solutions. Moreover, this procedure does not entail any discretization, linearization, or small perturbations and therefore reduces the computations a lot. Some examples are presented to establish the strength and applicability of this method. The results reveal that the method is very effective, straightforward, and simple to handle systems of boundary value problems.