  Numerical modeling of a mathematical model of the cosmological evolution of
an asymmetric scalar doublet with kinetic interaction between the components
was carried out. A wide range of values of fundamental parameters and initial
conditions of the model are considered. Various types of behavior have been
identified: models with an infinite inflationary past and future - with and
without a rebound point, models with a finite past and infinite future, with an
infinite past and finite future (Big Rip), as well as models with a finite past
and future. Based on numerical analysis, the behavior of models near the
initial singularity and the Big Rip is studied; it is shown that in both cases
the barotropic coefficient tends to unity, which corresponds to an extremely
rigid state of matter near singularities. A numerical example of the
cosmological generation of the classical component of a scalar doublet by its
phantom component is given. An assessment was made of the creation of the
velocity of fermion pairs by a scalar field near the rebound points and it was
shown that a scalar field at the cold stage of the Universe can ensure the
creation of the required number of massive scalarly charged fermions. Keywords:
cosmological model, phantom and classical scalar fields, quality analysis,
asymptotic behavior, numerical modelling, scalar field generation, types of
behavior.
