We consider the conservative finite difference method for calculating the three-dimensional stationary irrotational flows of an ideal gas. We solve the system of finite difference equations by the relaxation method, the convergence of iterations being accelerated by the method of least squares. We give the examples of calculations of sub-and supercritical flows around an ellipsoid of revolution.