We extend for the first time the applicability of the Optimal Homotopy Asymptotic Method (OHAM) to find approximate solution of a system of two-point boundary-value problems (BVPs). The OHAM provides us with a very simple way to control and adjust the convergence of the series solution using the auxiliary constants which are optimally determined. Comparisons made show the effectiveness and reliability of the method