A finite difference scheme is presented to solve the Fokker-Planck equation in (2+1) variables numerically. This scheme is applied to study stochastic beam dynamics in two-dimensional phase space.