ニュートン・コーツの公式(ニュートン・コーツのこうしき、)とは、等間隔の点における被積分関数の値に基づく数値積分法の総称である。名前はアイザック・ニュートンとロジャー・コーツに由来する。
ニュートン・コーツの公式は、等間隔の点での被積分関数の値が与えられた場合に有用である。もし他の点での値も求められるならば、ガウス求積やなどの他の方法の方が適している場合もある。
概要 ニュートン・コーツの公式は、端点を使う「閉じた」ものと、端点を使わない「開いた」ものの 2 種類に大別できる。
次の閉じたニュートン・コーツの公式は次のようになる。
ここで である。
は重みと呼ばれる。重みは以下のようにラグランジュ補間による補間多項式から導かれる。
また、以上の導出から重みは関数 によらず、 のみによって決まることがわかる。
次の開いたニュートン・コーツの公式の場合は、 とし、重みは閉じたものと同様である。
ニュートン・コーツの公式の一覧
+ 閉じたニュートン・コーツの公式 次数 名前 式 誤差項 1 台形公式 2 シンプソンの公式 3 シンプソンの3/8公式 4
+ 開いたニュートン・コーツの公式 次数 名前 式 誤差項 0 中点則 1 台形法 2 ミルンの公式 3
ここで、 は の略記である。
誤差項 は となる が存在することを意味する。 また、 の導関数の次数は、それ未満の次数の多項式が正確に積分できる(即ち、誤差が 0 になる)ことを示している。なお、 の次数と の導関数の階数は、1 つおきに 2 ずつ増加することに注意。
重みの計算 ニュートン・コーツの公式の重みは線形方程式系の解として求めることもできる。 これは補間多項式の一意性より が 次以下の多項式の場合 となることに基づく。係数行列はファンデルモンド行列である。
高次における不安定性
ニュートン・コーツの公式は、任意の次数で構築できる。しかし大きな次数 においてはルンゲ現象により誤差が の増加するにつれて指数関数的に大きくなる。そのため、通常は大きな次数ではガウス求積やなどの非等分点法の方が、安定してより正確な値を求められる。もしもそれらの方法を使えないならば、合成積分公式を使うことでルンゲ現象を避けることができる。高次の公式には積分の重みの中に負のものが含まれるなどの不自然さが伴う(Gaussの積分公式の重みは常に正である)。
合成積分公式 ニュートン・コーツの公式の精度を良くするには、ステップ長 は小さくする必要がある。つまり、積分区間 自体が小さくなければならない。このため、積分区間 を小さな部分区間に分割し、各部分区間ごとにニュートン・コーツの公式を使い、その結果を足し合わせるという方法が使われる。これは合成積分公式と呼ばれる。
関連項目 数値積分 内挿 台形公式 シンプソンの公式
参考文献
外部リンク
Newton–Cotes formulae on www.math-linux.com Newton–Cotes Formulae
Module for Newton–Cotes Integration, fullerton.edu Newton–Cotes Integration, numericalmathematics.com