Steady simple shear flows of smooth inelastic spheres are studied by means of a model kinetic equation and also of a direct Monte Carlo simulation method. Both approaches are based on the Enskog equation and provide for each other a test of consistency. The dependence of the granular temperature and of the shear and normal stresses on both the solid fraction and the coefficient of restitution is analysed. Quite a good agreement is found between theory and simulations in all cases. Also, simplified expressions based on the analytical solution of the model for small dissipation are shown to describe fairly well the simulation results even for not small inelasticity. A critical comparison with previous theories is carried out.