In the paper, we consider the simultaneous approximation of a collection of analytic functions by a collection of shifts of the Riemann zeta-function , where
is the Gram function and
are different positive numbers. It is obtained that the set of such shifts has a positive lower density.