We have mapped out the complete phase diagram of hard spherocylinders as a function of the shape anisotropy L/D. Special computational techniques were required to locate phase transitions in the limit L/D→ϱ and in the close-packing limit for L/D→0. The phase boundaries of five different phases were established: the isotropic fluid, the liquid crystalline smectic A and nematic phases, the orientationally ordered solids-in AAA and ABC stacking-and the plastic or rotator solid. The rotator phase is unstable for L/Dу0.35 and the AAA crystal becomes unstable for lengths smaller than L/DϷ7. The triple points isotropic-smectic-A-solid and isotropic-nematic-smectic-A are estimated to occur at L/D ϭ 3.1 and L/D ϭ 3.7, respectively. For the low L/D region, a modified version of the Gibbs-Duhem integration method was used to calculate the isotropic-solid coexistence curves. This method was also applied to the I-N transition for L/DϾ10. For large L/D the simulation results approach the predictions of the Onsager theory. In the limit L/D→ϱ simulations were performed by application of a scaling technique. The nematic-smectic-A transition for L/D→ϱ appears to be continuous. As the nematic-smectic-A transition is certainly of first order nature for L/Dр5, the tri-critical point is presumably located between L/Dϭ5 and L/Dϭϱ. In the small L/D region, the plastic solid to aligned solid transition is first order. Using a mapping of the dense spherocylinder system on a lattice model, the initial slope of the coexistence curve could even be computed in the close-packing limit.