In this work, a new solution is developed for the problem of contaminant transport in a single fracture-matrix system, where the first-order reaction rate constants are different in both fracture and matrix. It takes a form of convolution with three functions as a basis to consider different transport mechanisms separately. The statistical nature of the three functions, as well as the interpretation of the solution as a marginal probability distribution in the case of no first-order reactions, allows us to develop a simple Time Domain Random Wa...