In the present paper we develop an algorithm allowing to calculate line-of-sight velocity dispersions in an axisymmetric galaxy outside of the galactic plane. When constructing a self-consistent model, we take into account the galactic surface brightness distribution, stellar rotation curve and velocity dispersions. This algorithm is applied to a Sa galaxy NGC 4594 = M 104, for which there exist velocity dispersion measurements outside of the galactic major axis. The mass distribution model is constructed in two stages. In the first stage we construct a luminosity distribution model, where only galactic surface brightness distribution is taken into account. Thereafter, in the second stage we develop on the basis of the Jeans equations a detailed mass distribution model and calculate line-of-sight velocity dispersions and the stellar rotation curve. Here a dark matter halo is added to visible components. Calculated dispersions are compared with observations along different slit positions perpendicular and parallel to the galactic major axis. In the best-fitting model velocity dispersion ellipsoids are radially elongated. Outside the galactic plane velocity dispersion behaviour is more sensitive to the dark matter density distribution and allows to estimate dark halo parameters.