基于微积分理论近似计算定积分∫₋₁² e^(-x²)dx的方法探讨

嘿，这个积分挺有意思的——毕竟e^(-x²)没有初等原函数，所以只能靠微积分里的纯理论方法来近似。我来逐个拆解你提到的思路，再补充些实用的细节：

• 黎曼和（Riemann sums）：
  你说它看起来粗糙，但其实选对采样点的话精度能提上来。比如把[-1,2]分成n等份，每个小区间长度Δx=3/n，取每个区间的中点作为采样点，计算$\sum e^{-x_i2} \cdot \Delta x$。n取大一点（比如n=20），就能得到还不错的近似值。虽然它比不上后面的方法，但绝对是微积分里最基础、最贴合“纯理论”要求的近似思路。

• 利用全实数域积分结果√π：
  这个思路特别巧妙！我们知道$\int_{-\infty}^\infty e^{-x2}dx = \sqrt{\pi} \approx 1.77245$，那可以把目标积分拆成：
  $$\int_{-1}^2 e^{-x2}dx = \int_{-\infty}^\infty e^{-x2}dx - \int_{-\infty}^{-1} e^{-x2}dx - \int_{2}^\infty e^{-x2}dx$$
  剩下的两个无穷积分可以用分部积分做渐近展开来近似：
  比如对$\int_{a}^\infty e^{-x2}dx$，分部积分一次得：
  $$\int_{a}^\infty e^{-x2}dx = \frac{e^{-a2}}{2a} - \int_{a}^\infty \frac{e^{-x2}}{2x^2}dx$$
  再分部一次能得到更精确的展开式：
  $$\int_{a}^\infty e^{-x2}dx \approx \frac{e^{-a2}}{2a} \left(1 - \frac{1}{2a^2} + \frac{3}{4a^4} - \dots\right)$$
  代入a=1，得$\int_{-\infty}^{-1}e{-x^2}dx = \int_{1}^\infty e^{-x2}dx \approx \frac{e^{-1}}{2}(1-\frac{1}{2}) \approx 0.09197$；代入a=2，得$\int_{2}^\infty e^{-x2}dx \approx \frac{e^{-4}}{4}(1-\frac{1}{8}) \approx 0.00377$。这样目标积分≈1.77245 - 0.09197 - 0.00377≈1.6767，这个精度已经相当高了！

• 积分中值定理：
  积分中值定理告诉我们存在$c \in [-1,2]$，使得积分$=3 \cdot e^{-c2}$，但直接用的话范围太宽（e^{(-x²)在区间上的最大值是1，最小值是$e}{-4}$，所以积分在0.055到3之间）。不过可以拆分区间：把[-1,2]分成[-1,0]和[0,2]，分别用中值定理，得到的范围会窄一些，但还是不够精确。如果结合函数的单调性，用上下和来估计，能进一步缩小范围，但整体来说，这个方法的精度不如前面两个，更适合做粗略估计。

• 泰勒级数展开：
  这是最常用的纯微积分近似方法！我们知道$e^{t$的泰勒级数是$\sum_{n=0}}\infty \frac{t^{n}{n!}$，收敛于全体实数，代入$t=-x}2$得：
  $$e^{-x2} = \sum_{n=0}^\infty \frac{(-1)^n x^{2n}}{n!}$$
  因为这个级数在任何有限区间上一致收敛，所以可以逐项积分：
  $$\int_{-1}^2 e^{-x2}dx = \sum_{n=0}^\infty \frac{(-1)^n}{n!} \int_{-1}^2 x^{2n}dx$$
  计算积分$\int_{-1}^2 x^{2n}dx = \frac{2^{2n+1} + 1}{2n+1}$，所以最终积分的级数表达式是：
  $$\sum_{n=0}^\infty \frac{(-1)^n (2^{2n+1} + 1)}{n!(2n+1)}$$
  这是个交错级数，通项绝对值递减，所以前n项和的误差不超过下一项的绝对值。比如取前10项左右，就能得到非常接近真实值的结果（真实值约为1.6289）。如果觉得在[1,2]上收敛慢，还可以把区间拆成[-1,1]和[1,2]，在[1,2]上把$e^{-x2}$写成$e^{-1} \cdot e^{-(x2-1)}$再展开，收敛会更快。

备注：内容来源于stack exchange，提问作者ted