Parallel DO - Fortran

      program mm

      parameter (NRA =1000)
      parameter (NCA =1000)
      parameter (NCB =1000)

      integer 	cols,avecol,i,j,k
      real*8	a(NRA,NCA), b(NCA,NCB), c(NRA,NCB)


C     Initialize A and B 
      do 30 i=1, NRA
      do 30 j=1, NCA
        a(i,j) = 1.
 30   continue

      do 40 i=1, NCA
      do 40 j=1, NCB
	b(i,j) = 1.
 40   continue

C     Do matrix multiply

      do 50 k=1, NCB 
      do 50 i=1, NRA
      c(i,k) = 0.0
      do 50 j=1, NCA
           c(i,k) = c(i,k) + a(i,j) * b(j,k)
  50  continue

C     Print results 
       print*, 'Done '	
       write(*,*)'c(1,1)=',c(1,1)
       write(*,*)'c(NRA,NCB)=',c(NRA,NCB)

      end