C.2.3 プログラミングの仕方

初期値 $x_0$ が与えられたとき、漸化式 (C.3) によって、 数列 $\{x_j\}_{j=1,\cdots,N}$ を計算するプログラムはどう作ったらよいで しょうか？ここでは二つの素朴なやり方を紹介しましょう。

配列を使う方法

数列を配列で表現するのは、Fortran では自然な発想です。例えば

        integer MAXN
        parameter (MAXN = 1000)
        real x(0:MAXN)

のように配列 ``x'' を用意しておいて

        x(0) = x0
        t = a
        do j=0,N-1
          x(j+1) = x(j) + h * f(t,x(j))
          t = t + h
        end do

とするわけです。

配列を使わないですませる方法

漸化式 (C.3) を解くために、配列は絶対必要というわけで はありません^C.7。例えば、変数 ``x'' に各段階の $x_j$ の値を収めておくと して

        x = x0
        t = a
        do j = 0,N-1
          x = x + h * f(t,x)
          t = t + h
        end do

のようなプログラムで計算が出来ます。