In this paper we extend the study on the linear advection equation with independent stochastic velocity and initial condition performed in Dorini and Cunha (2011). By using both existing and novel results on the stochastic chain rule, we solve the random linear advection equation in the mean square sense. We provide a new expression for the probability density function of the solution stochastic process, which can be computed as accurately as wanted via Monte Carlo simulations, and which does not require the specific probability distribution of the integral of the velocity. This allows us to solve the non-Gaussian velocity case, which was not treated in the aforementioned contribution. Several numerical results illustrate the computations of the probability density function by using our approach. On the other hand, we derive a theoretical partial differential equation for the probability density function of the solution stochastic process. Finally, a shorter and easier derivation of the joint probability density function of the response process at two spatial points is obtained by applying conditional expectations appropriately.