next up previous contents
Next: C.2.3.0.1 重箱の隅をつつく注意: Up: C.2 数値解法(1) Previous: C.2.2 Euler(オイラー)法の紹介

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
のようなプログラムで計算が出来ます。




next up previous contents
Next: C.2.3.0.1 重箱の隅をつつく注意: Up: C.2 数値解法(1) Previous: C.2.2 Euler(オイラー)法の紹介
Masashi Katsurada
平成18年4月28日