Abstract
An efficient and accurate algorithm for Monte Carlo simulation of a solution to the diffusion equation with drift and branching terms is derived by expanding the Green function operator in time.