PROGRAM matriisi
  IMPLICIT NONE
  INTEGER, PARAMETER :: n = 5
  CHARACTER(LEN=*), PARAMETER :: form = '(A,5F7.2)'
  REAL, DIMENSION(n,n) :: matrix
  INTEGER :: i, j

  matrix = RESHAPE( (/ ( (i-j, i = 1, n), j = 1, n) /), &
      (/ n, n /) )
  WRITE(*,form) 'Matriisin alkioiden summa on', &
      SUM(matrix)
  WRITE(*,form) 'Matriisin sarakesummat ovat', &
      SUM(matrix,1)
  WRITE(*,form) 'Matriisin rivisummat ovat', &
      SUM(matrix,2)
END PROGRAM matriisi