Bernoulli polynomial and the numerical solution of high-order boundary value problems

In this work we present a fast and accurate numerical approach for the higher-order boundary value problems via Bernoulli collocation method. Properties of Bernoulli polynomial along with their operational matrices are presented which is used to reduce the problems to systems of either linear or nonlinear algebraic equations. Error analysis is included. Numerical examples illustrate the pertinent characteristic of the method and its applications to a wide variety of model problems. The results are compared to other methods.