program test_matev
use fvn_linear
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