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.