The second order BDF ADI difference scheme applied to the present equation is new. An ADI scheme based on second order convolution quadrature applied to the present equation is new. The numerical experiments are given under homogeneous and inhomogeneous Dirichlet boundary conditions.