We illustrate in the present paper how steady-state and unsteady-steady three-dimensional transport phenomena problems can be solved using the Chebyshev orthogonal collocation technique. All problems are elucidated using the ubiquitous software: MATLAB 1 . The treated case studies include: (1) unsteadystate heat conduction in a parallelepiped body, (2) quenching of a brick, (3) transient diffusion of a colorant in a liquid, and (4) steady-state natural convection in an enclosure. Whenever possible, the results obtained with the spectral method are bench-marked against the analytical solution in order to validate the proposed numerical technique. ß 2016 Wiley Periodicals, Inc. Comput Appl Eng Educ 24:866-875, 2016; View this article online at wileyonlinelibrary.com/journal/cae; 1 , POLYMATH, MATHCAD, JAVA, ASPEN-HYSYS or EXCEL is common in chemical engineering pedagogical literature [1 -12]. A literature survey over last three decades shows that only two papers report on the application of the Chebyshev orthogonal collocation technique, implemented in Mathematica © , for educational purposes [12,13], restricted however to one-and two-dimensional problems.Our main goal in this paper is to illustrate the utilization of this numerical method, commonly applied in research [14][15][16][17], in graduate-level classrooms in order to solve miscellaneous problems Correspondence to H. Binous (binoushousam@yahoo.com).