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.