1 subroutine quindet2(c,b,d,n,mu,no)
3 real (kind=8) ,dimension(n) :: c,b,d
4 real (kind=8) ,dimension(3,5) :: r
5 integer ::i,j,k,l,ll,rr,w,n,no
6 ! the arrays c,b i d contain the diagonal subdiagonal and subsubdiagonal elements of a simmetric quindiagonalmatrix. the value of the output parameter no is the number of eigenvalues greater than mu;
22 ! k count the major stages
28 !interchange row w and row 3 if necessary
29 if (dabs(r(3,1)).gt. dabs(r(w,1))) then
35 if (r(3,1).ge.0.eqv.r(w,1).ge.0) rr=-rr
37 !elimination of subdiagonal elements
44 r(3,j-1)=r(3,j)-x*r(w,j)
50 if (r(3,1).lt.0) rr=-rr
52 ! update elements of array r