“…It is important to emphasize that the analysis of second-order linear random differential equations has been performed in previous contributions using the random mean square calculus [5] (Ch. 4) (see [7][8][9][10][11][12] for the Airy, Hermite, Legendre, Laguerre, Chebyshev, and Bessel equations, respectively) and using other approaches like the random homotopy method [13,14], the Adomian method [15], the differential transformation method [16], the variational iteration method [17], etc. However, in all these contributions, only approximations for the two first statistical moments (mean and variance) of the solution were calculated.…”