This work investigates a new approach to find closed form solution to linear initial value problems (IVP). Classical Bernoulli polynomials have been used to derive a finite set of orthonormal polynomials and a finite ...