parameter (n=100)
      real*8 a(n,n),v(n,n),d(n),t
      integer np,nrot,i,j
      
      np=n

      do i=1, n
      do j=1, n
         a(i,j)=0.0d0
         if(j.eq.i-1) a(i,j)=1.d0
         if(j.eq.i+1) a(i,j)=1.d0
      enddo
      enddo

      call jacobi(a,n,np,d,v,nrot)

      do i=1,n
      do j=1,n-1
         if(d(i).gt.d(j)) then
            t=d(j)
            d(j)=d(i)
            d(i)=t
          do k=1, n 
            t=v(k,j)
            v(k,j)=v(k,i)
            v(k,i)=t
          enddo
         end if
      enddo
      enddo

      do i=1,n
         write(*,'(i4,4f12.8)') i, d(i),v(i,1),v(i,2),v(i,3) 
      enddo

      stop
      end