Fractional wave equations with attenuation have been proposed by Caputo [5], Szabo [28], Chen and Holm [7], and Kelly et al. [11]. These equations capture the power-law attenuation with frequency observed in many experimental settings when sound waves travel through inhomogeneous media. In particular, these models are useful for medical ultrasound. This paper develops stochastic solutions and weak solutions to the power law wave equation of Kelly et al. [11].