Equilibrium shapes of the self-gravitating gas cloud, founded by the serious of numeric experiments, based on three-dimensional nonstationary gas-dynamic model, are obtained in this work. It is shown that integral boundary condition is better for the Poison equation of gravity potential solving. Software system for solving of the Poison equation of gravity potential with different boundary condition on the external border and for solving of gas-dynamic equations, based on Roe scheme is developed. Founded equilibrium shapes can be used as initial data for three-dimensional numeric modeling of gas-dynamic processes at late stages of stars evolution.

Mathematical modelling in actual problems of science and technics