シンプソンズルール計算機

シンプソンの公式は、二次多項式を用いて定積分の値を近似する数値積分法です。長方形や台形で曲線下の面積を見積もるのではなく、連続する3つの標本点の各組に放物線を当てはめ、各放物線の下の正確な面積を計算します。これらの放物線の面積の合計が近似値となります。

この方法は18世紀のイギリスの数学者トーマス・シンプソンにちなんで名付けられましたが、その根本的なアイデアはそれ以前に他の人々の研究にも現れていました。その強さは単純な観察から生まれます。放物線は直線よりも滑らかな曲線にずっとよく沿うことができるのです。積分される関数が滑らかな場合、放物線の断片は真の曲線にぴったりと密着し、分割数が控えめでも驚くほど正確な近似が得られます。

なぜ3点が重要なのかを理解するために、直線は2点で完全に決まるのに対し、放物線（二次多項式）には3点が必要であることを思い出してください。シンプソンの公式は、積分区間を偶数個の部分区間に分割し、それらを2つずつ組みにし、各組に1本の放物線を当てはめることでこれを利用します。各組は2つの部分区間と3つの点（左端点、中点、右端点）にまたがります。

古典的な形は、式中に係数1/3が現れるためシンプソンの1/3公式と呼ばれます。3つの部分区間の組に三次多項式を用いるシンプソンの3/8公式もありますが、1/3公式が標準的な主力であり、このページの計算機を含むほとんどの計算機に実装されているバージョンです。

シンプソンの公式はニュートン・コーツの公式群に属します。これは、等間隔の点で被積分関数を評価し、それらの値を固定の重みで組み合わせることで積分を近似するものです。シンプソンの公式を際立たせるのはその重みパターンです。1, 4, 2, 4, 2, ..., 4, 1。端点は重み1、奇数位置の内点は重み4、偶数位置の内点は重み2を受けます。この交互パターンがこの方法の特徴であり、精度の源泉です。

直線区間の長さは簡単です。端点間の距離を測ればよいのです。曲線は違います。曲がったりねじれたりするため、単一の直線測定ではその真の広がりを捉えられません。曲線の長さを求めるために、数学者は曲線を微小な断片に分け、各断片を直線とみなして近似し、結果を足し合わせます。

曲がりくねった小道を歩くことを想像してください。とても小さな歩幅で歩けば、一歩一歩はほぼ直線です。道全体としては曲がっていてもです。これらの小さな直線の歩幅をすべて足し合わせれば、歩いた総距離のよい推定値が得られます。歩幅をどんどん小さくしていけば、推定値は正確な長さに収束します。これが弧長の核心的なアイデアです。

区間[a, b]上の関数y = f(x)について、微小な水平ステップdxを考えます。そのステップの間に、関数は約dy = f'(x) dxだけ上昇します。ここでf'(x)は導関数（傾き）です。曲線の微小区間の長さはピタゴラスの定理から、ds = √((dx)² + (dy)²) = √(1 + (f'(x))²) dxとなります。これらの無限小の断片をすべて積分で足し合わせると、弧長の公式が得られます。

L = ∫_a^b √(1 + (f'(x))²) dx

この公式の各部分には明確な意味があります。記号∫_a^bは、x = aからx = bまでの寄与を累積するよう指示します。式f'(x)は導関数で、各点における曲線の上昇の急さを測ります。これを二乗し、1を加えて平方根を取ることで、傾きは局所的な伸び率に変換されます。つまり、その位置で曲線が平らな水平線よりどれだけ長いかです。曲線が平らなところではf'(x) = 0で係数は1です。曲線が急なところでは係数は大きくなります。

この公式では関数が十分滑らかである必要があります。fは区間上で連続な導関数を持ち、被積分関数√(1 + (f'(x))²)は積分可能でなければなりません。微積分の授業や工学の仕事で出会うほとんどの関数では、これらの条件は満たされます。鋭い角、垂直な接線、不連続点では問題が生じます。そこでは導関数の振る舞いがおかしくなるからです。

難しいのは、この積分がしばしば初等的な原始関数を持たないことです。y = x²のような単純な関数でも、原始関数には逆双曲線関数が含まれます。多くの現実的な曲線には、閉形式の原始関数はまったく存在しません。まさにここでシンプソンの公式のような数値解法が不可欠になります。原始関数を必要とせず、積分を直接見積もるからです。

シンプソンの公式を弧長に適用するには、明確な手順の順序があります。重要な洞察は、シンプソンの公式が元の関数f(x)に適用されるのではないことです。弧長被積分関数g(x) = √(1 + (f'(x))²)に適用されるのです。この区別は重要であり、よくある誤りの原因でもあります。

完全な手順は次のとおりです。

1. 関数から始める。 y = f(x)を特定し、区間上で微分可能であることを確認します。
2. 導関数を求める。 f'(x)を解析的に求めるか、利用可能であることを確認します。
3. 被積分関数を構築する。 g(x) = √(1 + (f'(x))²)を作ります。これがシンプソンの公式で積分される関数です。
4. 区間を特定する。 測定したい曲線区間の下限aと上限bを確認します。
5. 分割数を選ぶ。 偶数n個の部分区間を選びます。シンプソンの公式はペアで動作するため、数は偶数でなければなりません。
6. ステップ幅を計算する。 h = (b - a) / nを計算します。これは各部分区間の幅です。
7. 標本点を求める。 点はx₀ = a, x₁ = a + h, x₂ = a + 2h、以下同様にx_n = bまでとなります。
8. 被積分関数を評価する。 各標本点でg(x₀), g(x₁), ..., g(x_n)を計算します。
9. 重みを適用する。 シンプソンのパターンで掛け合わせます。端点は1、奇数インデックスの内点は4、偶数インデックスの内点は2です。
10. 組み合わせてスケーリングする。 重み付きの値を合計し、h/3を掛けて弧長の推定値を得ます。

弧長の複合シンプソン1/3公式は次のとおりです。

L ≈ (h/3)[g(x₀) + 4g(x₁) + 2g(x₂) + 4g(x₃) + ... + 4g(x_n-1) + g(x_n)]

ここでh = (b - a)/n、g(x) = √(1 + (f'(x))²)です。

nが偶数でなければならないという要件は、ペア構造に由来します。各放物線は2つの部分区間（3点）にまたがります。nが奇数だと、最後の部分区間にパートナーがなくなり、方法は破綻します。何らかの理由で奇数個のパネルが必要な場合は、代わりに台形公式や混合アプローチが使われます。

重みパターンを注意深く見てください。最初と最後の値（g(x₀)とg(x_n)）はそれぞれ重み1で一度だけ現れます。奇数インデックスの内点（x₁, x₃, x₅, ...）は重み4を受けます。各点が放物線パネルの中点として機能し、放物線の影響が最も強いからです。偶数インデックスの内点（x₂, x₄, ...）は重み2を受けます。各点が隣接する2つのパネル間の共有端点として機能するからです。

複合シンプソンの公式にはいくつかの要素があり、それぞれが特定の役割を果たします。それらを個別に理解すれば、この方法ははるかに身近なものになります。

ステップ幅（h）。 各部分区間の幅であり、h = (b - a)/nとして計算されます。hが小さいほどパネルが多くなり、関数評価の回数が増える代わりに一般に精度が向上します。弧長の場合、hは曲線の標本化の細かさを制御します。

積分範囲（aとb）。 測定される曲線区間を定義します。区間の選択は結果に劇的に影響します。曲線に沿ってbをさらに延ばせば長さが加わります。常にa < bであり、両方が関数の定義域内にあることを再確認してください。

標本点。 x₀からx_nまでのn+1点は[a, b]上で等間隔です。シンプソンの公式はこれらすべての点で被積分関数を評価します。点は一つも飛ばされません。

端点の寄与。 g(x₀)とg(x_n)の値はそれぞれ重み1で一度だけ合計に入ります。境界で近似を固定します。

奇数インデックスの寄与。 x₁, x₃, x₅などの点は重み4を受けます。これらは各放物線パネルの中点です。重み4は、中点が放物線下の面積に最も大きな影響を与えるという事実を反映しています。

偶数インデックスの内部寄与。 x₂, x₄, x₆など（端点を除く）の点は重み2を受けます。これらの各点は隣接する2つのパネル間の端点として共有されるため、その総寄与は2倍になります。

重み付き合計。 各関数値にその重みを掛け、すべてを合計し、h/3を掛けます。係数h/3は重み付き合計を正しい単位と大きさにスケーリングします。

次の表はn = 6の場合の重みパターンを示しています。

このページの計算機は、いくつかの入力モードでシンプソンの公式を使って弧長を見積もります。各モードは曲線の異なる記述方法を扱います。どのモードが自分の問題に合うかを理解することが第一歩です。

関数モード。 x^2、sin(x)、exp(x)のような式と、区間の境界aとbを入力します。計算機は数値的に微分するか、提供された構造を使用し、弧長被積分関数を構築してシンプソンの公式を適用します。微積分の問題で最も一般的なモードです。

半径と角度モード。 円弧の場合、半径と中心角を入力します。弧長は幾何学から直接得られ、円の場合自体にはシンプソンの公式は不要ですが、比較にはこのモードが便利です。

パラメトリックモード。 パラメータtの関数としてx(t)とy(t)、およびパラメータ区間を入力します。計算機はパラメトリック弧長の公式を使用し、√((dx/dt)² + (dy/dt)²)にシンプソンの公式を適用します。

極座標モード。 r(θ)と角度区間を入力します。シンプソンの公式は角度範囲にわたって√(r² + (dr/dθ)²)を積分します。

点モード。 座標ペアのリストを入力します。計算機は連続する点間の直線距離を合計します。これはシンプソンの公式ではなく多角形近似で、式ではなく測定データがある場合に便利です。

3Dパラメトリックモード。 空間曲線のx(t)、y(t)、z(t)を入力します。シンプソンの公式は√((dx/dt)² + (dy/dt)² + (dz/dt)²)を積分します。

関数モードでシンプソンの公式による弧長計算を実行するには、次の手順に従います。

1. 関数を入力する。 標準的な記法で式を入力します。例：x^2またはsin(x)。
2. 下限を設定する。 曲線区間の開始x値aを入力します。
3. 上限を設定する。 終了x値bを入力します。a < bであることを確認してください。
4. 分割数を選ぶ。 偶数nを選びます。20や50などの中程度の値から始めましょう。
5. 方法を確認する。 このページでは数値解法がシンプソンの公式に固定されているため、選択は不要です。
6. 計算を実行する。 ツールはn+1点で被積分関数を評価し、推定値を返します。
7. 収束を確認する。 より大きなn（例えば2倍）で再実行し、比較します。結果がほとんど変わらなければ信頼できます。

現実的な例として、[0, 1]上でf(x) = x²をn = 20、次にn = 40で試してみてください。2つの結果は小数点以下数桁まで一致するはずで、推定値が安定していることが確認できます。

このセクションでは、各中間ステップを示しながら完全な例を解いていきます。関数は区間[0, 1]上のy = x²で、n = 4部分区間を使用します。すべての値は独立に計算されており、どの計算機でも検証できます。

ステップ1：関数。 y = x²という単純な放物線から始めます。

ステップ2：導関数。 微分するとy' = 2xとなります。

ステップ3：弧長被積分関数。 式に代入すると、g(x) = √(1 + (2x)²) = √(1 + 4x²)となります。

ステップ4：積分。 厳密な弧長はL = ∫₀¹ √(1 + 4x²) dxです。この積分は逆双曲線関数を含む閉形式を持ち、その値は約1.4789428575です。シンプソンの公式がどれだけ近づくかを見ていきましょう。

ステップ5：ステップ幅。 [0, 1]上でn = 4の場合、h = (1 - 0)/4 = 0.25です。

ステップ6：標本点。 等間隔の5点は次のとおりです。

• x₀ = 0.00
• x₁ = 0.25
• x₂ = 0.50
• x₃ = 0.75
• x₄ = 1.00

ステップ7：関数値。 各点でg(x) = √(1 + 4x²)を評価します。

• g(0.00) = √(1 + 0) = 1.0000000000
• g(0.25) = √(1 + 0.25) = √1.25 ≈ 1.1180339887
• g(0.50) = √(1 + 1) = √2 ≈ 1.4142135624
• g(0.75) = √(1 + 2.25) = √3.25 ≈ 1.8027756377
• g(1.00) = √(1 + 4) = √5 ≈ 2.2360679775

ステップ8：重みを適用する。 n = 4のシンプソンパターンは1, 4, 2, 4, 1です。

• 1 × g(x₀) = 1 × 1.0000000000 = 1.0000000000
• 4 × g(x₁) = 4 × 1.1180339887 = 4.4721359550
• 2 × g(x₂) = 2 × 1.4142135624 = 2.8284271247
• 4 × g(x₃) = 4 × 1.8027756377 = 7.2111025508
• 1 × g(x₄) = 1 × 2.2360679775 = 2.2360679775

ステップ9：重み付き合計。 これらを足すと1.0000000000 + 4.4721359550 + 2.8284271247 + 7.2111025508 + 2.2360679775 = 17.7477336082となります。

ステップ10：最終結果。 h/3 = 0.25/3 ≈ 0.0833333333を掛けます。L ≈ 0.0833333333 × 17.7477336082 ≈ 1.4789778007。

厳密値1.4789428575と比較すると、誤差は約0.0000349、つまりおよそ3.5 × 10^-5です。わずか4部分区間（5回の関数評価）で、シンプソンの公式はすでに小数点以下4桁の精度を達成しています。これがこの方法の強みです。滑らかな被積分関数に対する急速な収束です。

比較のために、n = 6に倍増すると1.4789405084（誤差 ≈ 2.4 × 10^-6）となり、粗いn = 2でも1.4821537045（誤差 ≈ 0.0032）となります。nを倍増するたびに誤差は劇的に減少します。これがシンプソンの4次収束の特徴です。

どの数値解法も厳密な答えではなく近似を生み出します。近似と真の値との差を誤差と呼びます。この誤差がどのように振る舞うかを理解することで、適切な設定を選び、結果を信頼できるようになります。

シンプソンの公式は4次法です。つまり、ステップ幅を小さくするにつれて誤差はh⁴（または等価に1/n⁴）に比例して減少します。実用的には、部分区間の数を2倍にすると誤差は通常約1/16になります。これが、上記の計算例がn = 2からn = 4、n = 6へと劇的に改善した理由です。

複合シンプソンの公式の古典的な誤差上界は次のとおりです。

|E_s| ≤ M(b - a)⁵ / (180n⁴)

ここで、Mは区間[a, b]上での被積分関数g(x)の4階導関数の絶対値の上界です。すなわち、[a, b]内のすべてのxに対してM ≥ |g⁽⁴⁾(x)|です。4階導関数は、被積分関数自体の曲率がどれほど急速に変化しているかを測ります。gが非常に滑らかな場合（4階導関数が小さい場合）、この上界は小さな誤差を保証します。gが急速に振動する場合、Mは大きく、上界はあまり安心できるものではありません。

いくつか重要な注意点があります。第一に、この上界ではgが区間上で連続な4階導関数を持つ必要があります。被積分関数に高階導関数の折れ点、カスプ、不連続点がある場合、この上界は適用されません。第二に、実際にMを計算することはしばしば困難であるため、この上界は日常の誤差推定よりも理論的理解に役立ちます。第三に、この上界は非常に悲観的になる可能性があります。実際の誤差は、式の予測よりもはるかに小さいことが多いのです。

実際には、最も信頼できる精度チェックは収束テストです。n分割で弧長を計算し、次に2n、4nで計算します。必要な小数点以下の桁数まで結果が一致すれば、推定値は信頼できます。結果が大きく変わり続ける場合は、nをさらに増やすか、被積分関数に問題のある振る舞い（急な曲がり、準特異点、急速な振動）がないか調べてください。

元の関数fの滑らかさは、gを通じて間接的に重要です。fに鋭い角があると、f'が不連続に跳び、g = √(1 + (f')²)は折れ点を継承します。シンプソンの公式は依然として機能しますが、収束は遅くなります。そのような場合は、問題のある点で区間を分割し、各滑らかな断片にシンプソンの公式を個別に適用することで、急速な収束が回復します。

台形公式はシンプソンの公式の最も近い親戚であり、最も一般的な代替手段です。どちらも等間隔の点を用いるニュートン・コーツの公式ですが、標本点の間で被積分関数をどのように近似するかが根本的に異なります。

台形公式は連続する標本点を直線で結び、台形を作ります。シンプソンの公式は3点の組に放物線を当てはめます。放物線は被積分関数に沿って曲がることができますが直線はできないため、シンプソンの公式は台形公式が見逃す曲率を捉えます。この幾何学的な違いは直接精度に反映されます。シンプソンの公式は4次（誤差は1/n⁴に比例）であるのに対し、台形公式は2次（誤差は1/n²に比例）です。

実用的には、これはシンプソンの公式がn = 20で達成できることを、台形公式が同等になるにはn = 200以上必要とすることが多いことを意味します。少なくとも滑らかな被積分関数については、です。代償として、シンプソンの公式は偶数のnと交互重みの少し慎重な帳簿管理を必要とします。台形公式は任意のnで動作し、より単純な重みパターン（1, 2, 2, ..., 2, 1）を持ちます。

どちらの方法も普遍的に優れているわけではありません。台形公式には、完全な周期にわたって積分される周期関数に対する驚くべき利点があり、指数関数的に速く収束できます。シンプソンの公式は、4次理論が適用されるほど滑らかでない関数では苦労することがあります。粗いデータやノイズの多いデータには、より単純な台形公式の方が頑健な場合もあります。

健全なワークフローは、シンプソンの公式を主要な方法として使用し、被積分関数の振る舞いが不確かな場合に高解像度で台形公式と照合することです。2つの方法が一致すれば、信頼性は高いです。大きく食い違えば、被積分関数にはより多くの分割数または区間の分割が必要である可能性が高いです。

曲線は複数の座標系で記述でき、それぞれに独自の弧長公式があります。シンプソンの公式はすべてに適用できます。各公式が滑らかな被積分関数の定積分だからです。

直交座標の弧長

[a, b]上のy = f(x)に対して、公式はL = ∫_a^b √(1 + (dy/dx)²) dxです。これは微積分の授業で教えられる標準形です。シンプソンの公式は、このガイド全体で説明したようにg(x) = √(1 + (f'(x))²)を積分します。fが連続な導関数を持って微分可能である限り、うまく機能します。

パラメトリック弧長

tがaからbまでのx = x(t)、y = y(t)で記述される曲線に対して、公式はL = ∫_a^b √((dx/dt)² + (dy/dt)²) dtです。パラメトリック形式は、y = f(x)として書けない円、楕円、サイクロイドなど、垂直線テストに失敗する曲線を扱えます。シンプソンの公式はパラメトリック速度関数√((x'(t))² + (y'(t))²)を積分します。これはtが増加するにつれて点が曲線に沿ってどれほど速く移動するかを表します。

極座標の弧長

θがαからβまでのr = r(θ)に対して、公式はL = ∫_α^β √(r² + (dr/dθ)²) dθです。極座標は渦巻き、カージオイド、バラ曲線を自然に記述します。被積分関数は動径距離rと角変化率dr/dθを組み合わせます。シンプソンの公式は直接適用でき、θについて積分します。

三次元の弧長

空間曲線x = x(t)、y = y(t)、z = z(t)に対して、公式は自然に拡張されます。L = ∫_a^b √((dx/dt)² + (dy/dt)² + (dz/dt)²) dtです。各座標はその二乗速度成分を寄与します。これは螺旋経路、3D軌道、コンピュータグラフィックスの曲線に不可欠です。シンプソンの公式は3D被積分関数を2D版と同じように扱います。

いずれの場合も、戦略は同じです。被積分関数を特定し、区間上で滑らかであることを確認し、偶数のnを選び、重み付き合計を適用します。このページの計算機は関数、パラメトリック、極座標、3Dパラメトリックの各モードをサポートし、それぞれが適切にシンプソンの公式を経由します。

弧長に対するシンプソンの公式は、数学モデルや測定データから曲線距離を見積もらなければならない多くの分野に登場します。

機械工学。 カムの輪郭、歯車の歯の曲線、ガイドレールは滑らかな数学的曲線として設計されます。その長さを見積もることは、材料の所要量、接触距離、摩耗パターンの決定に役立ちます。シンプソンの公式は設計の反復サイクル中に迅速な数値チェックを提供します。

道路と交通の幾何学。 高速道路のカーブ、鉄道の緩和曲線、滑走路の配置は、慎重に設計された幾何学的プロファイル（クロソイド、円弧、放物線の縦断曲線）に従います。これらの線形の真の長さを計算することは、コスト見積、標識の配置、建設計画にとって重要です。

CADとコンピュータ支援設計。 設計ソフトウェアは曲線をスプラインやベジェ曲線の断片として表現します。数値積分は、製造のためにこれらの複合曲線の長さを見積もります。切断経路や材料の使用は正確な長さデータに依存するからです。

ロボティクスと動作計画。 ロボットアームや自動運転車は、空間内の計画された軌道に従います。全経路長は移動時間、エネルギー消費、ミッション計画に影響します。シンプソンの公式はパラメータ化された軌道に沿った速度プロファイルを積分します。

科学計算。 粒子の軌道から場の線まで、物理現象をモデル化する研究者は日常的に曲線の長さを必要とします。支配方程式から閉形式の弧長が得られない場合、数値積分が標準的なアプローチです。

構造工学。 アーチ橋、曲線梁、シェル構造物の長さや表面の測定値は、その幾何学的定義から導かれます。正確な弧長計算は荷重解析と材料の積算を支えます。

数学モデリングと研究。 純粋数学と応用数学において、弧長は微分幾何学、変分問題、最適化に登場します。数値推定は理論的結果に対する直観と検証を提供します。

これらすべての文脈で、シンプソンの公式は比較的少ない関数評価で高い精度を提供するため評価されており、対話型ツールと大規模計算の両方に効率的です。

学生や実務家を陥れる繰り返し起こる誤りがいくつかあります。これらの落とし穴を認識することで、ほとんどの悪い結果を防げます。

奇数個の分割を使用する。 シンプソン1/3公式は根本的に偶数のnを必要とします。奇数のnではペア構造が崩れます。計算前に必ずnが偶数であることを確認してください。問題が奇数のカウントを示唆する場合は、次の偶数に切り上げます。

間違った導関数。 被積分関数はf'(x)を使用します。f(x)ではありません。よくある間違いは、誤って微分すること（例えば、x²の導関数を2xではなくxと書くこと）や、合成関数の連鎖律を忘れることです。すべての導関数を再確認してください。

gではなくfにシンプソンを適用する。 上記に関連して、一部のユーザーは弧長被積分関数g(x_i) = √(1 + (f'(x_i))²)ではなく、元の関数値f(x_i)に直接重みを適用します。これはまったく異なる積分を計算してしまいます。

誤った区間境界。 a > bを入力したり、関数の定義域外の境界を入力したりすると、意味のない結果が生じます。定義域が制限された関数（√(x)にはx ≥ 0が必要など）では、区間が有効であることを確認してください。

分割数が少なすぎる。 シンプソンの公式は速く収束しますが、非常に曲がった被積分関数にはn = 2やn = 4では不十分な場合があります。常により大きなnで収束チェックを行ってください。

度とラジアンの混同。 fに三角関数が含まれる場合、計算機がラジアンと度のどちらを期待しているか確認してください。混ぜると導関数と結果が静かに壊れます。ほとんどの数学ツールはラジアンを使用します。

厳密さを仮定する。 シンプソンの公式が厳密なのは3次までの多項式に対してのみです。一般の被積分関数では、結果は近似です。収束の証拠なしにそれを厳密とみなすのは間違いです。

特異点を無視する。 [a, b]内でf'が発散する場合（例えば、垂直な接線の近く）、被積分関数gは大きくなり、シンプソンの公式は遅く収束したり信頼できなくなったりすることがあります。区間を分割するか、特殊な方法を使用してください。

シンプソンの公式は数値弧長計算の優れたデフォルト選択ですが、それが輝くときと他の方法を探すべきときを知ることで、より効果的な問題解決者になれます。

シンプソンの公式を使うべきとき：

• 被積分関数g(x) = √(1 + (f'(x))²)が区間上で滑らかで、折れ点や特異点がない。
• 初等的な原始関数が存在しないか、見つけることが数値推定に見合う以上の労力を要する。
• 中程度の関数評価回数で良好な精度が必要。
• 授業の課題、工学のチェック、設計の検証のために弧長を計算している。
• 曲線がパラメトリックまたは極座標で記述され、結果の被積分関数が滑らか。

代替手段を検討するか、特に注意すべきとき：

• 被積分関数に区間内の不連続点や鋭い角がある。まず問題のある点で区間を分割する。
• 関数が急速に振動する。非常に大きなnが必要になるか、振動積分用に設計された方法が必要かもしれない。
• 滑らかな式ではなく、ノイズの多い測定データしかない。台形公式やスプラインに基づくアプローチの方が適切かもしれない。
• 極めて高い精度が必要。ガウス求積や適応的シンプソン公式（必要に応じて分割を自動的に細分化する）の方がよいかもしれない。
• 区間が無限または半無限。まず積分を有限区間に変換するか、特殊な求積を使用する。

ほとんどの微積分の学生や現役のエンジニアにとって、収束チェック（n、2n、4nの比較）を伴うシンプソンの公式は、必要なすべての精度を提供します。実装が簡単で、ソフトウェアに広く利用可能で、理論的にもよく理解されています。