5.苦しみどころ 変数の個数が違う h(n)の計算には{u(n)-u(n-1)}/dxが必要なのに u(n)の計算には{h(n+1)-h(n)}/dxが必要!! どうやら同じ添字を持った物理量が同じ座標にあるわけではないことが問題 配列添字がついていない、もともとの差分式をそのままプログラムとして書くことができればよろしいのではなかろうか? #そもそもfortranはFORmula TRANslaterなのでありました。