Indefinite integrals of products of special functions

A method is given for deriving indefinite integrals involving squares and other products of functions which are solutions of second-order linear differential equations. Several variations of the method are presented, which applies directly to functions which obey homogeneous differential equations. However, functions which obey inhomogeneous equations can be incorporated into the products and examples are given of integrals involving products of Bessel functions combined with Lommel, Anger and Weber functions. Many new integrals are derived for a selection of special functions, including Bessel functions, associated Legendre functions, and elliptic integrals. A number of integrals of products of Gauss hypergeometric functions are also presented, which seem to be the first integrals of this type. All results presented have been numerically checked with Mathematica.