The plasmaphysical code is presented that makes it possible to simulate the energy distributions of atoms escaping from plasma of the toroidal facilities. The code includes calculations of the particle spatial distributions in the charge-exchange plasma target consisting of atoms and hydrogen-like impurity ions. The calculations of the charge-exchange target are performed in two ways: using the Monte Carlo method and solving the ionization balance equation in the coronal approximation. The code can be used to interpret the neutral particle analyzers (NPAs) data. In particular, with its help, the ion temperatures determined using the NPAs can be corrected. The code can be also used to estimate the atomic fluxes from plasma onto the first wall expected at the currently designed fusion facilities.