Solutions of linear second order initial value problems by using Bernoulli polynomials

In this paper we use Bernoulli polynomials to derive a new spectral method to find the numerical solutions of second order linear initial value problems. Stability and error analysis of this method are studied. Numerical examples are presented which support theoretical results and provide favorable comparisons with other existing methods.