Qualitative Analysis and Numerical Simulation of Equations of the   Standard Cosmological Model: $\Lambda\not=0$
  On the basis of qualitative analysis of the system of differential equations of the standard cosmological model it is shown that in the case of zero cosmological constant this system has a stable center corresponding to zero values of potential and its derivative at infinity. Thus, the cosmological model based on single massive classical scalar field in infinite future would give a flat Universe. The carried out numerical simulation of the dynamic system corresponding to the system of Einstein - Klein - Gordon equations showed that at great times of the evolution the invariant cosmological acceleration has an oscillating character and changes from $-2$ (braking), to $+1$ (acceleration). Average value of the cosmological acceleration is negative and is equal to $-1/2$. Oscillations of the cosmological acceleration happen on the background of rapidly falling Hubble constant. In the case of nonzero value of the cosmological constant depending on its value there are possible three various qualitative behavior types of the dynamic system on 2-dimensional plane $(\Phi,\dot{\Phi})$, which correspond either to zero attractive focus or to stable attractive knot with zero values of the potential and its derivative. Herewith the system asymptotically enters the secondary inflation. Carried out numerical simulation showed that at cosmological constant $\Lambda<m^2 3\cdot10^{-8}$ the macroscopic value of the cosmological acceleration behaves itself similar to the case $\Lambda=0$, i.e. in the course of the cosmological evolution there appears a lasting stage when this value is close to $-1/2$ which corresponds to non-relativistic equation of state. In this article, the results of qualitative and numerical analysis, obtained in Yu. Ignat'ev, arXiv:1609.00745 [gr-qc], common to the case of a non-zero cosmological term. 
Yurii Ignat'ev
gr-qc