The Smooth Composite Chebyshev Finite Difference method is generalized for higher order initial and boundary value problems. Round-off and truncation error analyses and convergence analysis of the method are also extended to higher order. The proposed method is applied to obtain the highly precise numerical solutions of boundary or initial value problems of the Bratu and higher order Lane Emden types. To visualize the competency of the presented method, the obtained results are compared with nine different methods, namely, Bezier curve method, Adomian decomposition method, Operational matrix collocation method, Direct collocation method, Haar Wavelet Collocation, Bernstein Collocation Method, Improved decomposition method, Quartic B-Spline method and New Cubic B-spline method. The comparisons show that the presented method is highly accurate than the other numerical methods and also gets rid of the singularity of the given problems.(c) 2022 International Association for Mathematics and Computers in Simulation (IMACS). Published by Elsevier B.V. All rights reserved.