Based on Legendre–Gauss–Lobatto zeros and tensor product formulation, we provide a Legendre pseudo–spectral approach for solving one–dimensional parabolic advection–diffusion equations with constant parameters subject to stated beginning and boundary conditions in this chapter. To estimate the unknown function, we first employ differentiation matrices and their derivatives with respect to...