program test_matev
  use fvn
  implicit none
  real(8),dimension(3,3) :: a
  complex(8),dimension(3) :: evala
  complex(8),dimension(3,3) :: eveca
  integer :: status,i,j
  call init_random_seed()
  call random_number(a)
  a=a*100
  call fvn_matev(3,a,evala,eveca,status)
  
  write(*,*) "The matrix :"
  write (*,'(3(e12.5))') a
  write (*,*)
  do i=1,3
  write(*,'("Eigenvalue ",I3," : (",e12.5,",",e12.5,")  ")') i,evala(i)
  write(*,*) "Associated Eigenvector :"
  do j=1,3
  write(*,'("(",e12.5,",",e12.5,")  ")') eveca(j,i)
  end do
  write(*,*)
  end do
  end program