Question: Is it possible to solve numerically a 2D parabolic partial differential equation?

Let A a (linear) partial differential operator and f=f(x, y, t) a function of space (x, y) and time (t).
Let D_XY=[a, b] x [c, d], D_T=[0..e] and  d_XY the boundary of D_XY.
Is it possible, with Maple 2015, to solve numerically the pde Af=0 in the open domain  D_XY * D_T
with Dirichlet conditions on d_XY and initial condition f(D_XY, 0)?
It seems that the option 'numeric' of pdsolve doesn't work with more than 1 "non time" independent variable:
Error, (in pdsolve/numeric/process_PDEs) can only numerically solve PDE with two independent variables, got {t, x, y}

Nevertheles I seem to remember seeing this on Mapleprime (?)
If Maple 2015 can't do it, can Maple 2020?

Thanks in advance

• Remark: Af=0 is the heat equation with rho=Cp=lambda=1.
  I guess I could code the alternate directions method to transform  Af=0 into a sequence A_xf=0,  A_yf=0,  A_xf=0 ... of 1D diffusion equations but I'm a little bit lazzy and I'm waiting for your feedback before doing this.