rmination criterion is reached.
    </p>

    <p>
      <b>Note</b>: the original reference lists the input variable <b>F</b>
      as an <i>integer</i> workspace array.  However, <b>F</b> is used in the
      CLUSTR routine exclusively as a <i>real</i> array.  Even in single
      precision, this causes the routine to compute incorrect results (try it,
      please!); in double precision it also causes memory overwrites.
      The code presented here has corrected this mistake.
    </p>

    <h3 align = "c