This paper provides a full probabilistic solution of the randomized fractional linear nonhomogeneous differential equation with a random initial condition via the computation of the first probability density function of the solution stochastic process. To account for most generality in our analysis, we assume that uncertainty appears in all input parameters (diffusion coefficient, source term, and initial condition) and that a wide range of probabilistic distributions can be assigned to these parameters. Throughout our study, we will consider that the fractional order of Caputo derivative lies in ]0,1], that corresponds to the main standard case. To conduct our analysis, we take advantage of the random variable transformation technique to construct approximations of the first probability density function of the solution process from a suitable infinite series representation. We then prove these approximations do converge to the exact density assuming mild conditions on random input parameters. Our theoretical findings are illustrated through 2 numerical examples.
KEYWORDSfirst probability density function, random fractional differential equations, random variable transformation technique