subroutine arrout(a, m, n)
      implicit none
      integer m, n
      double precision a(m,n)
      integer i, j
           
      do i=1,m
        print *, (a(i,j), j=1,n)
      end do
      
      end