Simulating a partial differential equation — reaction-diffusion systems and Turing patterns