Java code for a spatial (regular Moore neighbourhood) Gillespie simulator and approximations to this in the form of four (discrete time) cellular automata algorithms