ex0714.f90

fortran 95书的源程序全集

program ex0714
implicit none
  integer I,J
  integer, parameter :: size = 5
  integer :: a(size,size)

  forall ( I=1:size, J=1:size, I>J )  a(I,J)=1 ! 上半部分
  forall ( I=1:size, J=1:size, I==J ) a(I,J)=2 ! 对角线部分
  forall ( I=1:size, J=1:size, I<J )  a(I,J)=3 ! 下半部分

  write(*,"(5(5I5,/))") a

  stop
end