1      subroutine prev(v,m,n,ndim)
2c
3c     ----- print out matrices
4c     vectors density etc.
5c
6      implicit double precision (a-h,o-z)
7      common /rdwr/ iun1,iun2,iun3,iun4,iun5
8      dimension v(ndim,*)
9
10      max  = 12
11      imax = 0
12
13100      imin = imax + 1
14         imax = imax + max
15         if (imax .gt. m) imax = m
16         write (iun3,9008)
17         write (iun3,8028) (i,i = imin,imax)
18         write (iun3,9008)
19         do j = 1,n
20            write (iun3,8048) j,(v(j,i),i = imin,imax)
21         end do
22         if (imax .lt. m) go to 100
23
24      return
25 9008 format(/)
26 8028 format(7x,12(3x,i3,3x))
27 8048 format(i5,2x,12f9.5)
28
29      end
30