In this study, an alternative numerical method having regard to the Bernstein operator is generated for approximate solutions of linear differential equations in the most general form under the initial and boundary conditions. Some applications are also revealed to show how the procedure can be performed for the problems.