この発言は
#00348　はんばあぐ      さんのRE:数値積分用Unit
に対するコメントです

この発言に対し以下のコメントが寄せられています

#00359　はんばあぐ      さん　RE^3:数値積分用Unit




　はんばあぐ さん、こんにちは。 Fermion です。o(^-^)o

　o(^o^o)(o^-^)o ワクワク！ 大変興味深いですね！！これから詳しく
拝見させていただきます。m(_"_)m ﾍﾟｺﾘ

　ところで、完全なお節介モードですみませんが、ちょっと気になりま
したので・・・。 ((((((^^;)

　台形で近似されている部分ですが、両端を除き、あるセルの下底部分
が次のセルの上底にあたりますので、それぞれ２回ずつ計算に絡んで来
ることになります。それゆえ結果的に、

         b
　Ｓ ＝ ∫f(x)dx
         a

　　　　　 f(a)    n-1　　　　f(b)  
　　 ≒ h{ ── ＋ Σf(xi) ＋ ── }
　　　　　　２     i=1　　　　 ２   

　ただし、 h ＝ ( b − a )／n  ,  xi ＝ a + ih  ( i = 1, 2, ... , n-1 )

　ということになりますので、

//=====================================================================
begin
  S := 0;
  xdelta := ( Fb - Fa ) / Fn;  //台形の高さ
  for i := 1 to Fn - 1 do
    begin
      xi := Fa + i * xdelta;
       S := S  + Ff(xi);
    end;
  S := S + ( Ff( Fa ) + Ff( Fb ) ) / 2.0;
  Result := xdelta * S;
end;
//=====================================================================

　のようにされた方がよろしいかも知れません。
　私 → (＠_☆ばき！＼(｀´;) 大きなお世話だ！

＃　大きなお世話レス、すみません。（外していましたらご免なさい。）

　m(_"_)m


                                     98/01/17(土) 00:30 Fermion [KHF03264]

Original document by Fermion         氏 ID:(KHF03264)

ここにあるドキュメントは NIFTY SERVEの Delphi Users' Forum の16番会議室「玉石混淆みんなで作るSample蔵」に投稿されたサンプルです。これらのサンプルはボーランド株式会社がサポートする公式のものではありません。また、必ずしも動作が検証されているものではありません。これらのサンプルを使用したことに起因するいかなる損害も投稿者、およびフォーラムスタッフはその責めを負いません。使用者のリスクの範疇でご使用下さい。