Calculation of the eigenvalues of a symmetric tridiagonal matrix by the method of bisection

The procedure bisect is designed to replace the procedures tridibi 1 and 2 given in [5]. All three procedures are based essentially on the following theorem.