This article discusses the process of spreading harmful impurities from traffic congestion. The process of spreading harmful impurities is considered by a three-dimensional model. The simulation of the spread of harmful impurities from traffic congestion is given, and algorithms for solving the three-dime nsional diffusion problem by the recurrent operator method are also given. In previous studies, we also used the recurrent operator method. The solution of the three-dimensional linear differential diffusion equation is sought in the form of a series. Where this series is inserted into each term of the linear differential equation, then making the indices shift in accordance with the recurrent operator method and taking out the common multiplier, we obtain the recurrent equation. Based on this recurrent equation, the coefficients of the series are searched for. Turning to the procedures for finding the coefficients of a series in the Maple software environment, the calculation of finding the coefficients was performed. Next, inserting the found coefficients into a series, we get the sum, which leads to a simpler solution convenient for differentiation. Further, partial derivatives in spatial coordinates are found based on software procedures in Maple. With this approach, the general solutions are expressed in terms of arbitrary functions and are not related to the solution of another equation. The obtained simulation results are shown in tables and illustrated in figures which are also obtained in the Maple software environment.