c      program eigen valuesq
***************************************************************************
      parameter (ipn=_)
      parameter (lwork=ipn*ipn)

      double precision a (ipn,ipn)
      double precision eval(ipn)
      double precision work(lwork)

      do 20 i=1,ipn
         do 20 j=1,ipn
            a(i,j)=0
 20   continue

      a(1,1)=____D0
      a(1,2)=____D0
      a(1,3)=____D0
      a(2,1)=____D0
      a(2,2)=____D0
      a(2,3)=____D0
      a(3,1)=____D0
      a(3,2)=____D0
      a(3,3)=____D0
c------------------------------------eigen value decomposition
      n   =ipn
      lda1=ipn
      call dsyev(___,___,_,_,____,____,____,_____,____)

      write(6,*) info
c-------------------------------------------------------------
c------------------------------------explained variance of each mode
      do 10 j=1,ipn
         write(6,'(4f10.5)') eval(j), (a(i,j),i=1,ipn)
 10   continue

      stop
      end