This paper proposes an efficient numerical method to obtain analytical-numerical solutions for a class of system of boundary value problems. This new algorithm is based on a reproducing kernel Hilbert space method. The analytical solution is calculated in the form of series in reproducing kernel space with easily computable components. In addition, convergence analysis for this method is discussed. In this sense, some numerical examples are given to show the effectiveness and performance of the proposed method. The results reveal that the method is quite accurate, simple, straightforward, and convenient to handle a various range of differential equations