In this work, we apply the reproducing kernel method to coupled system of second and fourth order boundary value problems. We construct a novel algorithm to acquire the numerical results of the nonlinear boundary-value problems. We also use the Legendre polynomials. Additionally, we discuss the convergence analysis and error estimates. We demonstrate the numerical simulations to prove the efficiency of the presented method.
