anchor.diff:  program EMCg ! version 1.0
   implicit none
   real(kind=8) :: eigenval, get_next_eigeval, get_2nd_deriv, get_mixed_deriv1, ft(3,3), ev(3,3)
  -real(kind=8) :: kp(3), dk, E(-2:2,-2:2,-2:2), A(4), m(3,3), b(3), bi(3), WORK(12), DUMMY(1,1)
  +real(kind=8) :: kp(3), dk, E(-2:2,-2:2,-2:2), A(4), m(3,3), b(3), bi(3), WORK(100), DUMMY(1,1)
   real(kind=8) :: trash1(3), trash2(3), tmp(3)
   
   real(kind=8), parameter :: a2b = 1.0D0/0.52917721092D0
  
positive.code: m(2,3) = m(3,2)
  
  write(*,*) "Original matrix:"
  do i=1,3
      write(*,fmt="(A1,F10.5,A2,F10.5,A2,F10.5,A2)") "{", m(1,i), ", ", m(2,i), ", ", m(3,i), "},"
  end do
  write(*,*)
  
  
  call inverse(m, 3)
  call DGEEV('N', 'V', 3, m, 3, b, bi, DUMMY, 1, ev, 3, WORK, 12, ok)
  
  write(*,*) "Eigensystem:"
  if (ok .eq. 0) then
      do i=1,3
          write(*,fmt="(F15.10,A2,F15.10,F15.10,F15.10)") b(i), ": ", ev(1,i), ev(2,i), ev(3,i)
      enddo
  else
      write (*,*) "An error occured while diagonalizing matrix in DGEEV()"
  endif
  
negative.code: character*20 wc1
  character(len=80) :: cha
  character(len=1) :: prg
  external DSYEV
  
  version_number='1.0'
  nkpoints = 61
  