Monte Carlo simulations are undertaken for hard ellipsoids for a number of values of its semi-axes in the range c / a ∈ { 0.1 , 10 } . The isotropic phase results are compared to the Vega isotropic phase equation of state for hard biaxial ellipsoids [C. Vega, Virial coefficients and equation of state of hard ellipsoids, Mol. Phys. 92 (1997) 651]. The Vega equation of state is found to perform exceedingly well for all of the models studied. The biaxial phase is seen to form only after the previous formation of a discotic phase.