Summary: Axisymmetric viscous fluid flows past spheroids are considered. The time-independent complete Navier-Stokes equations written in a spherical coordinate system are used for describing the flow. The problem is solved by the stabilization method on the basis of a variable direction scheme. The input coordinate system is transformed in order to construct a regular computational grid. As a result of the numerical investigation, the stream patterns of flow past elongated and oblate spheroids are obtained for various values of the determining parameters. Numerical values of the dimensions of the circulation zone and the drag coefficient are given for various values of the spheroid semi-axis ratio in the domain of moderate Reynolds numbers.