program call_matrix implicit none integer m ,i ,j parameter (m=1000) real*8 a(m,m)
call matrix( a, m ) do i = 1, m , 1 do j = 1, m , 1 write(*,*)"a(",i,",",j,")=",a(i,j) end do end do stop end program call_matrix
subroutine matrix( a, m ) implicit none integer m, i,j real*8 a(1:m,1:m), b(1:m,1:m) do i = 1, m , 1 do j = 1, m , 1 a(i,j) = dble( i+j ) b(i,j) = dble( i*j ) end do end do return end subroutine matrix
