In recent years spectral methods have been found to be a powerful tool for the numerical solution of hynamic differential equations. The main attraction of spectral method is accuracy even though it is generally difficult to implement and solve the complex problems using spectral method. We introduced diffusion equations describing the state of air pollution and solved by pseutospectral method in dimensionless form. The results were compared with both those of other numerical methods and analytical solutions. Comparing with finite difference method and finite element method, spectral method shows the highest accuracy for one dimension problem in this study. Also, the results of two dimensional diffusion problems show good agreement with analytical solutions.