subroutine testmat(n,a)
c
      integer n, i, j
      double precision a(n,n)

      do 10 i = 1, n
        do 20 j = 1, n
          print *, 'a(', i, ',', j, ') =', a(i,j)
   20   continue
   10 continue
      return
      end