The method of hydrodynamic potentials is applied to sojving three-dimensional stationary viscous incompressible fluid flow problem. The solution is defined as a sum of the double layer potential with the unknown density $\varphi$ and the volume potential having nonlinear terms. The value of $\varphi$ is determined from a linear integral second kind equation with a weak singularity. The equation under consideration admits a unique solution for any right-hand side substituted from the previous iteration. The presence of Reynolds number $R$ as a multiplier in the volume potential provides convergence of iterations for sufficiently small values of $R$. The flow around a three-axial ellipsoid was considered as an example.