文摘
We present a novel theoretical approach to the analysis of adaptive quadratures and adaptive Simpson quadratures in particular which leads to the construction of a new algorithm for automatic integration. For a given function \(f\in C^4\) with \(f^{(4)}\ge 0\) and possible endpoint singularities the algorithm produces an approximation to \(\int _a^bf(x)\,{\mathrm d}x\) within a given \(\varepsilon \) asymptotically as \(\varepsilon \rightarrow 0\). Moreover, it is optimal among all adaptive Simpson quadratures, i.e., needs the minimal number \(n(f,\varepsilon )\) of function evaluations to obtain an \(\varepsilon \)-approximation and runs in time proportional to \(n(f,\varepsilon )\).