A novel discriminant approximation of periodic differential equations

A new approximation of the discriminant of a second order periodic differential equation is presented as a recursive summation of the evaluation of its excitation function at different values of time. The new approximation is obtained, at first, by means of Walsh functions and then, by using some algebraic properties the dependence on the Walsh functions is eliminated. This new approximation is then used to calculate the boundaries of stability. We prove that by letting the summation elements number approach to infinite, the discriminant approximation can be rewritten as a summation of definite integrals. Finally we prove that the definite integrals summation is equivalent to the discriminant approximation made by A. Lyapunov which consists in an alternating series of coefficients defined by multiple definite integrals, that is, a series of the form $A=A_0?A_1+\dots +{(?1)}^n{A}_n$, where each coefficient $A_n$ is defined as an n-multiple definite integral.