program mxm
	parameter (n=512)
	double precision a(n,n), b(n,n), c(n,n)
C$unroll 4
	do i=1,n
	do j=1,n
	do k=1,n
	c(i,j) = c(i,j) + a(i,k)*b(k,j)
	enddo
	enddo
	enddo
	print *, c(2,3)
	end