An innovative harmonic numbers operational matrix method for solving initial value problems.

In this paper a novel operational matrix of derivatives of certain basis of Legendre polynomials is established. We show that this matrix is expressed in terms of the harmonic numbers. Moreover, it is utilized along with the collocation method for handling initial value problems of any order. The convergence and the error analysis of the proposed expansion are carefully investigated. Numerical examples are exhibited to confirm the reliability and the high efficiency of the proposed method. [ABSTRACT FROM AUTHOR]