We study scalar perturbations of a recently found 3+1-dimensional FLRW quantum space-time solution in Yang-Mills matrix models. In particular, the linearized Schwarzschild metric is obtained as a solution. It arises from a quasi-static would-be massive graviton mode, and slowly decreases during the cosmic expansion. Along with the propagating graviton modes, this strongly suggests that 3+1 dimensional (quantum) gravity emerges from the IKKT matrix model on this background. For the dynamical scalar modes, non-linear effects must be taken into account. We argue that they lead to non-Ricci-flat metric perturbations with very long wavelengths, which would be perceived as dark matter from the GR point of view.Since the notorious problems in attempts to quantize gravity arise primarily from the Einstein-Hilbert action, is is very desirable to find another framework for gravity, which is more suitable for quantization. String theory provides such a framework, but the traditional approach using compactifications leads to a host of issues, notably lack of predictivity. This suggests to use matrix models as a starting point, and in particular the IKKT or IIB model [2], which was originally proposed as a constructive definition of string theory. Remarkably, numerical studies in this nonperturbative formulation provide evidence [3-5] that 3+1-dimensional configurations arise at the non-perturbative level, tentatively interpreted as expanding universe. However, this requires a new mechanism for gravity on 3+1-dimensional non-commutative backgrounds as in [1], which does not rely on compactification. The present paper provides further evidence and insights for this mechanism. The (linearized) Schwarzschild metric is clearly the benchmark for any viable theory of gravity. There has been considerable effort to find noncommutative analogs of the Schwarzschild metric from various approaches, leading to a number of proposals [6-9] and references therein, cf. also [10]; however, none is truly satisfactory. The proposals are typically obtained by some ad-hoc modification of the classical solution, without any intrinsic role of noncommutativity, which is put in by hand. In contrast, the quantum structure (or its semi-classical limit) plays a central role in the present framework. Our solution is a deformation of the noncommutative background which respects an exact SO(3) rotation symmetry, even though there are only finitely many d.o.f. per unit volume. The solution has a good asymptotics at large distances, allowing superpositions corresponding to arbitrary mass distributions. In fact we obtain generic quasi-static Ricci-flat linearized perturbations, which complement the Ricci-flat propagating gravitons found in [1]. This realization of the (linearized) Schwarzschild solution is remarkable and may seem surprising, because the action is of Yang-Mills type, and no Einstein-Hilbert-like action is required 1 . This means that the theory has a good chance to survive upon quantization, which is naturally defined via integration over the s...