  Consideration of the latest experimental data on the searches for extended
sector of Higgs bosons produced at the LHC at a center-of-mass energy of 13
TeV, allows for computer modeling of the properties of supersymmetric particles
within 2HDM model. The experimental restrictions on model parameters accounted
in FeynHiggs code that is implemented in SusHi program, gave us the possibility
to calculate the cross sections and branching fractions for three mechanisms of
production and decay of Higgs bosons: 1) pp$\rightarrow H\rightarrow\tau\tau$,
2) pp$\rightarrow A\rightarrow Zh\rightarrow llbb$, 3) pp$\rightarrow
H\rightarrow hh\rightarrow bb\tau\tau$ at a center-of-mass energy of 14 TeV.
The considered computer modelling make it possible to draw conclusions about
the need to take into account the b-associated production process of Higgs
bosons for fermionic decay channel at large values of tan$\beta$. Differential
cross sections with respect to the Higgs transverse momentum $p_t$ and
pseudorapidity $\eta$ are calculated and the peculiarities of the kinematics of
the Higgs boson decay products are recognized.
