Describes an efficient algorithm for the calculation of specified eigenvalues of Ax = λBx with band symmetric A and B, the latter being positive definite. Each eigenvalue is isolated using the Sturm sequence property of leading principal minors of A − λB and is then computed accurately using a modified version of successive linear interpolation.