in Eq. (\[approx-rho\]) - (\[approx-q\]), respectively. Deceleration parameter $q$ seems to be positive everywhere.

Results {#final}
=======

In present work, we have considered the evolution of a full causal bulk viscous flat, isotropic and homogeneous Universe with bulk viscosity parameters and equation of state taken from recent lattice QCD data and heavy-ion collisions. Three classes of solutions of the evolution equation have been obtained.

In Fig. \[Figg1\], $H(t)$ and $a(t)$ are depicted in dependence on the comoving time $t$. We compare $H(t)$, given by Eq. (\[eq:mysolut1\]), and $a(t)$, given by Eq. (\[approx-a\]), with the counterpart parameters obtained in the case when the background matter is assumed to be an ideal and non-viscous fluid, described by the equations of state of the non-interacting ideal gas, $$\begin{aligned}
H(t) &=& \frac{1}{2t} \label{htideal1}, \\
a(t) &=& \sqrt{t}. \label{atideal1}\end{aligned}$$

In the left panel of Fig. \[Figg1\], $H(t)=\dot a/a$ has an exponential decay, whereas in the non-viscous case, $H(t)$ is decreasing according to Eq. (\[htideal1\]). The latter is much slower than the former, reflecting the nature of the exponential and linear dependencies. The other difference between the two cases is obvious at small $t$. We notice a divergence, or singularity, associated with the ideal non-viscous fluid, 