在Fortran 95中将2行x 3列的矩阵数组写入输出文本文件

program write_matrix
  implicit none
  
  integer, parameter :: rows = 2
  integer, parameter :: cols = 3
  real :: matrix(rows, cols) = reshape([1.0, 2.0, 3.0, 4.0, 5.0, 6.0], [rows, cols])
  integer :: i, j
  character(len=20) :: filename = "output.txt"
  integer, parameter :: status = 0
  integer :: unit
  
  ! Open the file for writing
  open(newunit=unit, file=filename, status='replace', action='write', iostat=status)
  
  if (status /= 0) then
    write(*, *) "Error opening the file"
    stop
  end if
  
  ! Write the matrix to the file
  do i = 1, rows
    do j = 1, cols
      write(unit, '(f8.2)') matrix(i, j)
    end do
    write(unit, *)
  end do
  
  ! Close the file
  close(unit)
  
  write(*, *) "Matrix written to", trim(filename)
  
end program write_matrix