program arraymultiply
      real A(5,5), B(5,5), C(5,5)
      integer imax,jmax,kmax
      data ((A(i,j),j=1,2),i=1,3)/1,2,4,6,1,0/
      data ((B(i,j),j=1,3),i=1,2)/2,1,4,3,0,7/
      imax=3
      jmax=3
      kmax=2
      do i=1,imax
        do j=1,jmax
          c(i,j)=0.0
          do k=1,kmax
            c(i,j)=c(i,j)+a(i,k)*b(k,j)
          enddo
        enddo
      enddo
!
!   The next lines print the result using implied do loop
!
      print 10,((c(i,j),j=1,jmax),i=1,imax) 
10    format(' C(i,j) is',/,(3(F4.1,1x)))
      end