On the reduction of a symmetric matrix to tridiagonal form

A stable algorithm for reducing a symmetric, non-definite matrix of ordern to tridiagonal form, involving aboutn3/6 additions and multiplications is presented.