A Monte Carlo method for large reaction–diffusion systems