  We address a detailed non-perturbative numerical study of the scalar theory
on the fuzzy sphere. We use a novel algorithm which strongly reduces the
correlation problems in the matrix update process, and allows the investigation
of different regimes of the model in a precise and reliable way. We study the
modes associated to different momenta and the role they play in the ``striped
phase'', pointing out a consistent interpretation which is corroborated by our
data, and which sheds further light on the results obtained in some previous
works. Next, we test a quantitative, non-trivial theoretical prediction for
this model, which has been formulated in the literature: The existence of an
eigenvalue sector characterised by a precise probability density, and the
emergence of the phase transition associated with the opening of a gap around
the origin in the eigenvalue distribution. The theoretical predictions are
confirmed by our numerical results. Finally, we propose a possible method to
detect numerically the non-commutative anomaly predicted in a one-loop
perturbative analysis of the model, which is expected to induce a distortion of
the dispersion relation on the fuzzy sphere.
