Abstract
In this study, a new collocation method based on Bernstein polynomials defined on the interval [a, b] is introduced for approximate solutions of initial and boundary value problems involving higher order linear differential equations with variable coefficients. Error analysis of the method is demonstrated. Some numerical solutions are given to illustrate the accuracy, efficiency and implementation of the method, and the results of the proposed method are also compared with the other methods in several examples.