Computing the density function of complex models with randomness by using polynomial expansions and the RVT technique. Application to the SIR epidemic model | Publicación