The non-linear Boltzmann equation has been solved for shock waves in a gas of elastic spheres. The solutions were made possible by the use of Nordsieck's Monte Carlo method of evaluation of the collision integral in the equation. Accurate solutions were obtained by the same method for the whole range of upstream Mach numbers M^ from 1.1 to 10 even though the corresponding degree of departure from equilibrium varies by a factor greater than 1000. Many characteristics of the internal structure of the shock waves have been calcu lated from the solutions and compared with Navier-Stokes, Mott-Smith and Krook descriptions which, except for low Mach numbers, are not based upon the Boltz mann equation itself. Among our conclusions are the following: 1. The reciprocal shock thickness is in agreement with that of the 2 Mott-Smith shock (u-moment) from M^ of 2.5 to 8. The density profile is asymmetric with an upstream relaxation rate (measured as density change per mean free path) approximately twice as large as the downstream value for weak shocks and equal to the downstream value for strong shocks. 2. The temperature density relation is in agreement with that of the Navier-Stokes shocks for the lower Mach numbers in the range of 1"1 to 1.56. The Boltzmann reciprocal shock thickness is smaller than the Navier-Stokes value at this range of Mach number because the viscosity-temperature relation computed is not constant as predicted by the linearized theory. 2 3. The velocity moments of the distribution function are, like the Mott-Smith shock, approximately linear with respect to the number density; however, the deviations from linearity are statistically significant. The four functionals of the distribution function discussed show maxima within the shock. 4. The entropy is a good approximation to the Boltzmann function for all M^. The solutions obtained satisfy the Boltzmann theorem for all Mach numbers. The increase in total temperature within the shock is small, but the increase is significantly different from zero. 5. The ratio of total heat flux q to (associated with the longi tudinal degree of freedom) correlates well with local Mach number for all M^ in accord with a relation derived by Baganoff and Nathenson. The Chapman-Enskog linearized theory predicts that the ratio is constant. The (effective) transport coefficients are larger than the Chapman-Enskog equivalents by as much as a factor of three at the mid-shock position. 6. At M 1=4, and for 40% of velocity bins, the distribution function is different from the corresponding Mott-Smith value by more than three times the 90% confidence limit. The rms value of the percent difference, in distribu tion function is 15% for this Mach number. The halfwidth and several other characteristics of the function Jfdw^dvz differ from that of the Chapman-Enskog first iterate, and many of the deviations are in agreement with an experiment by Muntz and Harnett. 7. The ratio of the collision integral (found from our solution of the Boltzmann equation) to that calculated from Mott-Smit...