The extension of the previous paper [Can. J. Phys. Vol. 88, (2010), 501–511] has been made. Therefore, the effect of the neutral atoms collisions with electrons and with positive ions is taken into consideration, which was ignored, for the sake of simplicity, in the earlier work. Thus, we will have multi-collision terms (electron–electron, electron–ion, electron– neutral) instead of one term, as was studied before for the sake of facilitation. These collision terms are needed to obtain the real physical situation. The new procedures will increase the ability of the research applications. This study is based on the solution of the BGK (Bhatnager–Gross–Krook) model of the nonlinear partial differential Boltzmann equations coupled with Maxwell’s partial differential equations. The initial-boundary value problem of the Rayleigh flow problem applied to the system of the plasma (positive ions + electrons+ neutral atoms), bounded by a moving plate, is solved. For this purpose, the traveling wave solution method is used to get the exact solution of the nonlinear partial differential equations system. The ratios between the different contributions of the internal energy changes are predicted via the extended Gibbs equation for both dia-magnetic and para-magnetic plasma. The results are applied to a typical model of laboratory argon plasma. 3D-Graphics illustrating the calculated variables are drawn to predict their behavior and the results are discussed.