  The asymptotic expansion method is extended by using currently available
accurate values for the first ten virial coefficients for hard sphere fluids.
It is then used to yield an equation of state for hard sphere fluids, which
accurately represents the currently accepted values for the first sixteen
virial coefficients and compressibility factor data in both the stable and the
metastable regions of the phase diagram.
