In this paper we use uniform cubic spline polynomials to derive some new consistency relations. These relations are then used to develop a numerical method for computing smooth approxi-mations to the solution and its first second as well as third derivatives for a second order boundary value problem. The proesent method out-performs other collocations finite-difference and splines methods of the same order. numerical illustratiosn are provided to demonstrate the practical use of our method.