Numerical solution of a two-dimensional IHCP based on Duhamel's principle