After a brief presentation of random matrices as a random surface QFT approach to 2D quantum gravity, we focus on two crucial mathematical physics results: the implementation of the large N limit (N being here the size of the matrix) and of the double-scaling mechanism for matrix models. It is worth emphasizing that, in the large N limit, it is the planar surfaces which dominate. In the third section of the chapter we introduce tensor models, seen as a natural generalization, in dimension higher than two, of matrix models. The last section of the chapter presents a potential generalisation of the Bollobás–Riordan polynomial for tensor graphs (which are the Feynman graphs of the perturbative expansion of QFT tensor models).