くだすれFORTRAN（超初心者用）

くだすれFORTRAN（超初心者用） at TECH

374:デフォルトの名無しさん
06/06/02 19:58:25
微妙に写し間違えているｗ
DIMENSION NJUN(3)
...
S=( FX(A) + FX(B) ) / 2.0
DO 20 J = 1, N - 1
X = A + H * J
S = S + FX(X)
20 CONTINUE
S = S * H
WRITE(6, *) 'BUNKATSU = ', NJUN(I), 'MENSEKI =', S
10 CONTINUE
STOP
END
こんな感じになるのではないかな．
台形公式を>>364の様に素直に小さな台形の和を足して求めると，
関数の値を２回同じ値に対して求めなければならない．

そこを改良すると，ここでの例題のように，両端だけを半分にして，
残りを一回計算して足し上げてゆけば良いことになる．
絵を描いて考えればすぐ分かると思う．

NJUNは分割数だろう．何個の台形に分けるかを指定している．
数学的には分割が細かいほど正確な値になるが（積分は分割→０の極限で定義
されている）現実には数値誤差や計算時間の関係からほどほどの分割数がいい．
いくつがいいかは時と場合によるので一概には言えない．

次ページ

続きを表示

1を表示