A new approach to find an approximate solution of linear initial value problems with high degree of accuracy

Abstract

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 operational matrix to simplify derivatives in IVP. These orthonormal polynomials together with the operational matrix of relevant order provides a robust approximation to the solution of a linear initial value problem by converting the IVP into a set of algebraic equations. Depending upon the nature of a problem, a polynomial of degree n or numerical approximation can be obtained. The technique has been demonstrated through four examples. In each example, obtained solution has been compared with available exact or numerical solution. High degree of accuracy has been observed in numerical values of solutions for considered problems.