The immune system in vertebrates is composed of individual cells called lymphocytes which work together to combat antigens such as bacteria and viruses. Upon detecting foreign molecules these immune cells secrete soluble factors that attract other immune cells to the site of the infection. The soluble factors are governed by a system of reaction-diffusion equations with sources that are centered on the cells. The motion of the cells is inherently stochastic, but biased toward the gradient of the soluble factors. I will discuss a numerical method for solving the reaction-diffusion stochastic system based on a first order splitting scheme. This method makes use of known first order schemes for solving the diffusion, the reaction and the stochastic differential equations separately.