  By the example of a mathematical model of a biochemical process, the
structural instability of dynamical systems is studied by calculating the full
spectrum of Lyapunov indices with the use of the generalized Benettin
algorithm. For the reliability of the results obtained, the higher Lyapunov
index determined with the orthogonalization of perturbation vectors by the
Gram-Schmidt method is compared with that determined with the overdetermination
of only the norm of a perturbation vector. Specific features of these methods
and the comparison of their efficiencies for a multidimensional phase space are
presented. A scenario of the formation of strange attractors at a change of the
dissipation parameter is studied. The main regularities and the mechanism of
formation of a deterministic chaos due to the appearance of a fold or a funnel,
which leads to the uncertainty of the evolution of a biosystem, are determined.
