[2.1]配列にデータを格納
1.D0
2.D0
3.D0
4.D0
5.D0
と縦に並んだデータをファイル（retsu.datとする）から読み, 配列a(5)に格納するプログラム例：
implicit none
integer i
real*8 a(5)
open(1,file='retsu.dat')
do i=1,5
read(1,*) a(i)
write(*,*) a(i)
enddo
end
1.D0 2.D0 3.D0 4.D0 5.D0
と横に並んだデータをファイル（gyou.datとする）から読み, 配列a(5)に格納するプログラム例：
implicit none
real*8 a(5)
integer i
open(1,FILE='gyou.dat')
read(1,*) (a(i),i=1,5)
write(*,*) (a(i),i=1,5)
end
[2.2]二次元配列
配列の要素を指定する番号は複数あってもよい. 二つになると二次元配列となり, 三つになると三次元配列となる. FORTRANでは7次元の配列まで使用可能である. 配列，変数の名前は，英字から始まる31文字以内の任意の英数字を用いてつける．英字については，大文字と小文字の区別は行わない．そのため，ABC(10)とabc(10)は同じ配列となる．
二次元配列の宣言:
宣言文 配列名（寸法の下限：寸法の上限, 寸法の下限：寸法の上限）
例1
real*8 a(10,10)
a(1,1)からa(10,10)まで要素を100個持つ倍精度実数型二次元配列を宣言
例2
integer i(2,-1:8)
i(1,-1)からi(2,8)まで要素を20個もつ実数型二次元配列を宣言
二次元配列を使ったプログラム例：
implicit none
integer a(5,5),b(5,5),c(5,5),n,i,j
n=0
do i=1,5
do j=1,5
n=n+1
a(i,j)=0
b(i,j)=n
enddo
enddo
do i=1,5
do j=1,5
c(i,j)=b(i,j)
enddo
enddo
do i=1,5
write(*,*) (a(i,j),j=1,5)
enddo
write(*,*)
do i=1,5
write(*,*) (b(i,j),j=1,5)
enddo
write(*,*)
do i=1,5
write(*,*) (c(i,j),j=1,5)
enddo
end
[2.3]二次元データの読み取り
以下の数値をファイル matrix.dat にコピーし, それぞれの数値を二次元配列a(2,3)に格納するプログラムを作成せよ. 確認のため画面に値を表示させること．
1.D0 2.D0 3.D0
4.D0 5.D0 6.D0
配列に初期値を代入しておくには，data文を用いるやりかたと，配列式を用いるやりかたがある．すべての要素に同じ値を代入するときは，配列式が便利である．部分的に異なる値を代入しておく場合は，data文を用いる．
[3.1]DATA文
配列，変数を宣言した後に記述する．構文は以下のとおり．
data 変数や配列などの並び / 定数の並び /
例：
implicit none
integer i,j
real*8 a(3,3),b
data a,b / 4*0.D0, 5*1.D0, 10.D0 /
do i=1,3
write(6,*) i,j,a(i,j)
enddo
write(6,*) b
end
4*0.D0と書くと，始めの４つの要素にゼロを代入する．引き続いて5*1.D0で残りの５つの要素に1.D0が代入される．上の例で，どの要素に0.D0が，どの要素に1.D0が格納されているか注意すること．また，変数や配列の要素の数と，定数の数が異なるとエラーになるので注意すること．
[3.2]配列式
配列を用いて式を記述することができる．要素の数が同じであれば，加減乗除が可能．しかし，各要素についての演算を行うだけなので，次節で取り扱う通常の意味での行列の積にはなっていないので注意．
例：
implicit none
integer i,j
real*8 a(3,3),b(3,3),c(3,3)
a=1.D0
b=2.D0
c=a+b
do i=1,3
write(*,*) (c(i,j),j=1,3)
enddo
write(*,*)
c=a*b
do i=1,3
write(*,*) (c(i,j),j=1,3)
enddo
write(*,*)
c=a/b
do i=1,3
write(*,*) (c(i,j),j=1,3)
enddo
end
上の例では，配列a,bの全ての要素に1.D0と2.D0がそれぞれ代入され，配列の要素ごとに足し算，掛け算，割り算が行われる．通常の意味での行列の積や逆行列との積とはまったく異なるので注意すること．
[4.1]行列を取り扱う際のポイント
- [超重要] 数式とみかけが同じになるようにプログラムを組むこと！
- 変数，配列の名前は式と同じ記号にする
- 添字についての和はループに対応する．
例：行列とベクトルの積
例えば，3×3の行列Aと，3次元ベクトルBの積を考える．この積はベクトルであり，それをCとする．数式で書くと以下になる．
 |
(1) |
式通りにプログラムを書くと次になる．
do j = 1,3
C(i) = A(i,j) * B(j)
end do
[4.2]行列の掛け算
先ほどの例と同様に，行列AとBとの積を考える．この結果は行列となり，それをCとすると次の式になる．この式を参考に行列の掛け算を行う．
 |
(2) |
hai1.dat
2.D0 3.D0 4.D0
5.D0 6.D0 7.D0
8.D0 9.D0 10.D0
hai2.dat
1.D0 2.D0 3.D0
4.D0 5.D0 6.D0
2.D0 4.D0 6.D0
上の二つの3行3列の行列をそれぞれ hai1.dat，hai2.dat とする．
練習問題
- ファイルから数値を読み込み，それぞれ配列a(3,3),b(3,3)に値を格納しなさい．
- それぞれの行列の対角和を求めなさい．
- この行列の掛け算を行い，結果を画面に出力するプログラムをつくりなさい.
- 配列aに格納された行列（hai1.dat)を10乗して結果を画面に出力するプログラムをつくりなさい.