過去にここのブログにて、C++言語を利用してGauss-Legendre quadrature の積分点を25桁まで求めたことを示すURLを紹介したが、改めて紹介する。上図は64個の積分点のデータを計算した結果です。下記URLに記載があります。作者は日本（埼玉県）にて暮らしていたことのある方です。
当方ではExcl_VBAを利用して15桁までの範囲であれば200個程度の積分点を求めるプログラムを公開済みです。必要な方は下記URLより入手されてください。
下図は当方作成の画面、関心のある方は上図の25桁と出力値を比較されてみてください。上図で紹介したURLのそれは積分点数が２^n(n=1,2,3....)に限られますが、当方のソフト(15桁出力）は積分点数は1024点以上の自然数でも正確な数値を出力します。（1024点数まで精度確認済み）
先日海外から C++Gauss-Legendre quadrature の積分点を求めるプログラムを作成中だが当方のVBAコードにとても興味があるとのメールを未知の方から頂いた。返信に上記URL（http://www.holoborodko.com/pavel/numerical-methods/numerical-integration/）
を紹介してあげた。
0 件のコメント:
コメントを投稿