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

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

768:デフォルトの名無しさん
06/07/25 15:12:27
>>766
（２）　
解析的には8/3+2｜^2_0=4.6666・・・・・・・だが、１００分割で計算すると4.6668000
台形公式だしｈ＾２に誤差が比例しているとするとまぁこんなところか。

文関数を使っているがあまりお勧めできない。積分をサブルー朕に、被積分関数をｆｕｎｃｔｉｏｎにすると
もっと通っぽいだろう。あと、台形の面積を正直に求めているが、これは同じ値を２回計算しているので
ウンコである。両端を０．５倍、その他を１倍で足すほうがまとも。

初心者っぽさを出すために、マジ台形で計算してみた。

PROGRAM okki
F(x) = x**2 + 1.0
n = 100
s = 0.0
x0 = 0.0
x1 = 2.0
h = (x1 - x0) / n
DO 10 i = 1, n
x = x0 + h * REAL(i - 1)
s = s + h * ( F(x) + F(x + h) ) / 2.0
10 CONTINUE
WRITE(*, *) ' AREA =', s
STOP
END

次ページ

続きを表示

1を表示