Tau spectral methods and Adomian’s decomposition can be fruitfully combined to quickly approximate the analytical solution of any time-dependent partial differential equation with boundary conditions defined on the four sides of a rectangle. In this work, combinations of Legendre polynomials have been used to generate orthogonal two-dimensional polynomials on a rectangular domain. The time evol...