Direct Monte Carlo simulations of the Enskog-Boltzmann equation for a spatially uniform system of smooth inelastic spheres are performed. In order to reach a steady state, the particles are assumed to be under the action of an external driving force which does work to compensate for the collisional loss of energy. Three different types of external driving are considered: (a) a stochastic force, (b) a deterministic force proportional to the particle velocity and (c) a deterministic force parallel to the particle velocity but constant in magnitude. The Enskog-Boltzmann equation in case (b) is fully equivalent to that of the homogeneous cooling state (where the thermal velocity monotonically decreases with time) when expressed in terms of the particle velocity relative to the thermal velocity. Comparison of the simulation results for the fourth cumulant and the high energy tail with theoretical predictions derived in cases (a) and (b) [T. P. C. van Noije and M. H. Ernst, Gran. Matt. 1, 57 (1998)] shows a good agreement. In contrast to these two cases, the deviation from the Maxwell-Boltzmann distribution is not well represented by Sonine polynomials in case (c), even for low dissipation. In addition, the high energy tail exhibits an underpopulation effect in this case.