A method to derive general standard and null Lagrangians for second-order differential equations whose solutions are special function of mathematical physics is presented. The general null Lagrangians are used to find the corresponding general gauge functions. All derived Lagrangians are new and in special cases they reduce to those published in the literature. The obtained results are applied to the Bessel, Hermite and Legendre equations, which have many applications in physics, applied mathematics and engineering.