program ch0811
  implicit none
  integer, parameter :: nrow = 5
  integer, parameter :: ncol = 6
  real, dimension (1:nrow*ncol) :: results = &
    (/ 50, 47, 28, 89, 30, 46, &
       37, 67, 34, 65, 68, 98, &
       25, 45, 26, 48, 10, 36, &
       89, 56, 33,     45, 30, &
       65, 68, 78, 38, 76, 98, 65 /)
  real, dimension (1:nrow, 1:ncol) :: exam_results = 0.0
  real, dimension (1:nrow) :: people_average = 0.0
  real, dimension (1:ncol) :: subject_average = 0.0

  exam_results = reshape(results, (/nrow,ncol/), (/0.0,0.0/), (/2,1/) )
  exam_results(1:nrow, 3) = 2.5*exam_results(1:nrow, 3)
  subject_average = sum(exam_results, dim=1)
  people_average = sum(exam_results, dim=2)
  people_average = people_average/ncol
  subject_average = subject_average/nrow
  print *, ' People averages'
  print *, people_average
  print *, ' Subject averages'
  print *, subject_average
end program ch0811