program wmatrix
      implicit none
      integer*4 LDA,LWORK,N,i,j,INFO
      parameter (LDA=20,LWORK=LDA)
      COMPLEX*16 a(LDA,LDA),WORK(LDA*LDA)
      REAL*8 RWORK(3*LDA),W(LDA)


      n=5
      a(1,1)=(1.0d0,0.0d0)
      a(2,2)=(2.0d0,0.0d0)
      a(3,3)=(3.0d0,0.0d0)
      a(4,4)=(4.0d0,0.0d0)
      a(5,5)=(5.0d0,0.0d0)

      a(1,2)=(0.0d0, 1.0d0)
      a(1,3)=(0.0d0, 2.0d0)
      a(1,4)=(0.0d0, 5.0d0)
      a(1,5)=(0.0d0,-1.0d0)
      a(2,3)=(0.0d0, 3.0d0)
      a(2,4)=(0.0d0, 6.0d0)
      a(2,5)=(0.0d0, 0.0d0)
      a(3,4)=(0.0d0, 7.0d0)
      a(3,5)=(0.0d0, 3.0d0)
      a(4,5)=(0.0d0,37.0d0)

      do 2 i=1,n
      do 2 j=1,i-1
2     a(i,j)=-a(j,i)
      
      call ZHEEV('V', 'U', N, A, LDA, W, WORK, LWORK, RWORK,
     $                  INFO )

      if(INFO.EQ.0)then
       write(6,*)'Eigenvalues:'
       do 1 j=1,n
1      write(6,*)j,w(j)
      endif
      if(INFO.lt.0)then
       write(6,*)INFO,' invalid argument number',abs(INFO)
       stop
      endif
      if(INFO.gt.0)then
       write(6,*)INFO,' convergence problem'
       stop
      endif
      end

      include 'blas.f'
      include 'lapackdiag.f'