A new perspective on single and multi-variate differential equations

We elaborate upon a new method of solving linear differential equations, of arbitrary order, which is applicable to a wide class of single and multi-variate equations. Our procedure separates the operator part of the equation under study in to a part containing a function of the Euler operator and constants, and another one retaining the rest. The solution of the equation is then obtained from the monomials (or the monomial symmetric functions, for the multi-variate case), which are the eigenfunctions of the Euler operator. Novel exponential forms of the solutions of the differential equations enable one to analyze the underlying symmetries of the equations and explore the algebraic structures of the solution spaces in a straightforward manner. The procedure allows one to derive various properties of the orthogonal polynomials and functions in a unified manner. After showing how the generating functions and Rodriguez formulae emerge naturally in this method, we briefly outline the generalization of the present approach to the multi-variate case.