We present nearly self-consistent stellar-halo models of Fornax dwarf spheroidal galaxy associated with the Milky Way galaxy. Such satellites are dominated by dark matter and have almost no gas in the system. Therefore, they are excellent objects for N-body modelling that takes into account visible and dark matter halo components. We perform seven dynamical models in total based on two N-body source code sets. Three of them are made by self-consistent realizations containing a spheroidal body modeled by truncated King or Plummer models for stellar components and the lowered Evans dark halo model. These components are given by the distribution functions which are the functions of energy $E$ and angular momentum $L_z$ only. The others are made by the AGAMA code based on the calculation of orbits. To guess the parameters which make the N-body models close to the visible object we use hydrodynamic stellar-dark model of the Fornax galaxy taking into account the velocity anisotropy parameter. The AGAMA models show the best agreement of the resulting velocity dispersion profiles with the observed data.