A gas of long thin rods undergoes an order-disorder phase transition as a function of rod concentration. We have evaluated the critical concentration at which this first-order transition occurs using Onsager's hard-core interaction model. We obtain the nematic angular distribution function of rods in the ordered phase expanded in a series of Legendre polynomials.