Numerical simulations of stochastic circulatory models