  By choosing appropriate deformed Maxwellian ion and electron distribution
functions depending on the two particle constants of motion, i.e. the energy
and toroidal angular momentum, we reduce the Vlasov axisymmetric equilibrium
problem for quasineutral plasmas to a transcendental Grad-Shafranov-like
equation. This equation is then solved numerically under the Dirichlet boundary
condition for an analytically prescribed boundary possessing a lower X-point to
construct tokamak equilibria with toroidal sheared ion flow and anisotropic
pressure. Depending on the deformation of the distribution functions these
steady states can have toroidal current densities either peaked on the magnetic
axis or hollow. These two kinds of equilibria may be regarded as a bifurcation
in connection with symmetry properties of the distribution functions on the
magnetic axis.
