In this paper, we will apply cubic B-splines on a uniform mesh to explore the numerical solutions and numerical derivatives of a class of nonlinear second-order boundary value problems with two dependent variables. Our new method is based on the cubic spline interpolation. The analytical solutions and any-order derivatives can be well approximated with 4th order accuracy. Furthermore, our new m...