program test_det use fvn_linear use Kind_Definition implicit none real(kind=dp_kind),dimension(3,3) :: a real(kind=dp_kind) :: deta integer :: status,i call init_random_seed() call random_number(a) a=a*100 deta=fvn_det(3,a,status) do i=1,3 write (*,'(3(f15.5))') a(i,:) end do write (*,*) write (*,*) "Det = ",deta,"Status = ",status end program