Stochastic simulation of the spatio-temporal dynamics of reaction-diffusion systems: the case for the bicoid gradient