A self-consistent gasdynamic model of the solar wind interaction with the local interstellar medium (LISM), which took into account the mutual influence of the plasma component (electrons and protons) of the LISM and the LISM H atoms that penetrate into the heliosphere was constructed by Baranov et al. (1981) in the approximation of axial symmetry. This model, however, had a number of defects. In particular, the motion of the H atoms was described by hydrodynamical equations, although the mean free path of the H atoms and the characteristic length of the problem were comparable. An iterative method, that used a Monte Carlo simulation of H atom motion in the field of the plasma component hydrodynamic parameters, was suggested by Baranov et al. (1991) and only the first step of the iteration was realized (non-self-consistent problem solution). In this paper the results of the self-consistent problem solution for a single set of the undisturbed solar wind and LISM parameters are presented. The structure of the upwind as well as wake regions of the flow is calculated. The geometrical pattern of the flow (bow shock, heliopause, termination shock, Mach disc, etc.), the bulk velocity and the number densities of H atoms and plasma component are obtained and analyzed as a function of the distance from the Sun for different values of the polar angle. The effects of resonance charge exchange of the LISM H atoms as well as energetic H atoms ''born'' in the solar wind are taken into account. It is interesting to note that the effect of H atoms penetrating the solar wind results in the disappearance of the complicated flow structure as well as the supersonic regions between the heliopause and termination shock in the downwind region. In future we are going to compare our theoretical results with the results of Voyager 1/2, Pioneer 10/11, Ulysses spacecraft, and other experiments. ¿ American Geophysical Union 1993