Numerical Methods For Reaction Diffusion Systems In High Spatial Dimensions