単位行列

RESHAPEで作ろうとして失敗していたが・・
それにつけても、宣言時に初期値を与えると自動的に SAVE 属性になるとは知らなかった。(ソースは"Fortran95/2003 Explained" chapter 6。)

PROGRAM test
IMPLICIT NONE
INTEGER :: i, n(10, 10) = 0
FORALL (i = 1:10) n(i, i) = 1
PRINT '(10i4)', n
STOP
END PROGRAM test