A Lagrangian method for deriving new indefinite integrals of special functions

A new method is presented for obtaining indefinite integrals of common special functions. The approach is based on a Lagrangian formulation of the general homogeneous linear ordinary differential equation of second order. A general integral is derived which involves an arbitrary function and therefore yields an infinite number of indefinite integrals for any special function which obeys such a differential equation. Techniques are presented to obtain the more interesting integrals generated by such an approach, and many integrals, both previously known and completely new can be derived using the method. Sample results are given for Bessel functions, Airy functions, Legendre functions and complete elliptic integrals. Much more extensive results for specific special functions will be presented separately. Integrals can be derived which combine common special functions as separate factors.