In this study, a new collocation method based on Bernstein polynomials defined on the interval [a, b] is introducedfor approximate solutions of initial and boundary value problems involving higher order linear differential equations with variable coefficients. Error analysis of the…