On the numerical solution of Brownian motion processes