c     examples/mult/mult_seq
c     multiplying two matrices, a and b
c
      parameter (IDIM=20)
      real*8 a(IDIM,IDIM),b(IDIM,IDIM),c(IDIM,IDIM)
      real*8 aa(IDIM),bb(IDIM),cc
      open(8,file='in.dat',status='unknown')
      read(8,*) n
      do i=1,n
         read(8,*) (a(i,j),j=1,n)
      enddo
      do i=1,n
         read(8,*) (b(i,j),j=1,n)
      enddo
      do i=1,n
         do j=1,n
            do k=1,n
               aa(k)=a(i,k)
               bb(k)=b(k,j)
            enddo
            call comput(n,aa,bb,cc)
            c(i,j)=cc
         enddo
      enddo
      do i=1,n
         write(*,*) (c(i,j),j=1,n)
      enddo
      stop
      end

      subroutine comput(n,aa,bb,cc)
      parameter (IDIM=20)
      real*8 aa(IDIM),bb(IDIM),cc
      cc=0.0
      do k=1,n
         cc=cc+aa(k)*bb(k)
      enddo
      return
      end