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