This paper is devoted to the study of the Cubic B-splines to find the numerical solution of linear and non-linear 8th order BVPs that arises in the study of astrophysics, magnetic fields, astronomy, beam theory, cylindrical shells, hydrodynamics and hydro-magnetic stability, engineering, applied physics, fluid dynamics, and applied mathematics. The recommended method transforms the boundary problem to a system of linear equations. The algorithm we are going to develop in this paper is not only simply the approximation solution of the 8th order BVPs using Cubic-B spline but it also describes the estimated derivatives of 1st order to 8th order of the analytic solution. The strategy is effectively applied to numerical examples and the outcomes are compared with the existing results. The method proposed in this paper provides better approximations to the exact solution.