We analyze the algorithm used to generate the value for the machine
epsilon in the Eispack suite of routines and show that it can fail on a
binary floating-point system. The comments in the code describing the
conditions under which this method will work are not restrictive enough
and we provide a replacement set of assumptions.  We conclude by
suggesting how the algorithm may be modified to overcome most of the
shortcomings.