Note: This article is available in Chinese only. 本文暂无英文版本。
View original
16沈阳的阴影还在orz，来学习一下辛普森积分。
参考资料：梯形多步法和辛普森积分
辛普森积分是一种数值积分方法（然后现在只记得教计算方法的是一个小姐姐，并不记得当时学了什么orz
大概就是用二次抛物线近似计算曲边梯形面积，辛普森积分公式如下：
$$ \int_a^b f(x) \, dx \approx \frac{b-a}{6} \left(f(a) + 4 f\left(\frac{a+b}{2}\right) + f(b)\right) $$下面放代码：
1double f(double x){return sin(x)*x;}//这是被积函数
2double simpson(double l,double r){return (r-l)*(f(l)+f(r)+4*f((l+r)/2))/6;}
3double di(double l,double r){//越二分以得到更精确的结果
4 double m=(l+r)/2;
5 double ans=simpson(l,r);
6 if(sgn(ans-simpson(l,m)-simpson(m,r))==0)return ans; //在误差之内
7 return di(l,m)+di(m,r); //不再误差之内就继续将区间缩小
8}切了个练手题，慢慢补充总结。
需要注意的是，simpson对精度要求比较高。。。eps开到1E-10才过。。
所以问题就在于写出积分公式（如果是多重积分要变成累次积分？orz)