情報処理演習 第2回 連立 1 次方程式 「Fortran での 1 次元配列表示」 (講義ノート)

計算機が苦手な細かい数字を取り扱うことと人間の目から見て分かりやすい形との違いを見るために、 整数型と倍精度型で、いくつか配列の定義の仕方を見てみる。 gfortran は、Fortran 90/95 に対応しており、 新しい記法が可能となっている。

	
      program vector
      implicit none
      integer :: i
      real(8) x(2)
      !integer x(2)
      x(0) = 1.0d0
      x(1) = 2.0d0
      !x(0) = 1
      !x(1) = 2
      write(*,*) (x(i), i = 0, 1)
      do i = 0, 1
         x(i) = x(i) * 2
      enddo 
      write(*,*) (x(i), i = 0, 1)
      do i = 0, 1
         x(i) = x(i) / 4
      enddo 
      write(*,*) (x(i), i = 0, 1)
      end program vector

また、下のように :: を導入することにより、初期値を代入することもできる。

	
      program vector
      implicit none
      real(8) :: x(0:1) = (/ 1.0d0, 2.0d0 /)
      !integer :: x(0:1) = (/ 1, 2 /)
      integer i
      write(*,*) (x(i), i = 0, 1)
         : 
         : 
      end program vector

コメントを変更し、倍精度の実数と整数型の違いを調べておこう。