Solving the dynamics of equilibrium dense liquids is a notoriously difficult problem [1]. In the low-density limit, or at short times, the solution is obtained via kinetic theory, while at any density, hydrodynamics provides a description of the dynamics at large length and time scales. However, at high densities or low temperatures, kinetic theory breaks down, and the hydrodynamic regime is pushed to scales that are much larger than any experimentally relevant scale. In this strongly interacting "supercooled liquid" regime, dynamics is slow, viscosity is high, and both correlation functions and transport coefficients display non-trivial behavior that is neither captured by kinetic theory nor by hydrodynamics. The only microscopic description of this regime is obtained by the Mode-Coupling Theory (MCT) [2,3], which is a set of closed equations derived from a series of poorly controlled approximations of the true dynamics. Despite its non-systematic nature [4][5][6][7], MCT accurately describes the initial slowing down of liquid dynamics upon supercooling, including the wavevector dependence of correlation functions [8].Interestingly, in the formal limit in which the spatial dimension d goes to infinity, liquid thermodynamics reduces to the calculation of the second virial coefficient [9][10][11][12]. Based on this observation, a first attempt to solve exactly liquid dynamics for d → ∞ was presented in [13] (see also [14]), but the full dynamical mean field theory (DMFT) that describes exactly the equilibrium dynamics in d → ∞ was only derived recently, via a second-order virial expansion on trajectories [15] or via a dynamic cavity method [16,17]. The DMFT provides a set of closed one-dimensional integro-differential equations, which exactly describe the many-body liquid dynamics in the thermodynamic limit, and are similar in structure to those obtained for quantum systems in the same limit [18]. We refer the reader to [17,[19][20][21] for a detailed review of the solution of liquid dynamics in infinite dimensions, including its extension to the out-of-equilibrium setting.Unfortunately, the analytical solution of the DMFT equations is out of reach. In this work, we present their numerical solution, obtained through an iterative method and a straightforward discretization of time. This strategy, however, only works for differentiable interaction potentials. We thus discuss how to derive the DMFT of hard spheres via a non-trivial limit of a soft sphere interaction. We present numerical results for soft and hard spheres, supported by analytical computations at low densities, in the short time limit, and at long times in the glass phase.The paper is organized as follows. In section II we review the basic equations of the DMFT of liquids. In section III, we discuss the Hard Sphere limit of DMFT. In section IV, we discuss the discretization and convergence algorithms used in this work. In section V, we present the numerical solution for Soft and Hard Spheres. Finally, in section VI we draw our conclusions and present s...