这种算法使用c语言，简单易懂，高效求积算法龙贝格，具体思路请参考程序

已知积分立，所以我们可以通过对该被积函数的数值积分来计算的近似值。 （1）分别利用三点、五点和七点Gauss-Legendre和Gauss-Chebyshev求积公式来计算的近似值（首先求出Gauss点）； （2）将积分区间分为四等分，分别用复合三点、五点和七点Gauss-Legendre和Gauss-Chebyshev求积公式来计算的近似值（每个子区间分别利用Gauss-Legendre和Gauss-Chebyshev求积公式）； （3）比较及分析计算结果。

1. 用不同数值方法计算积分 (1) 取不同的步长h. 分别用复合梯形及复合辛普森求积计算积分, 给出误差中关于h的函数, 并与积分精确值比较两个公式的精度, 是否存在一个最小的h, 使得精度不能再被改善? (2) 用龙贝格求积计算完成问题(1). (3) 用自适应辛普森积分, 使其精度达到10−4. 附录1 复合梯形求积MATLAB程序 附录2复合辛普森求积MATLAB程序 附录3龙贝格求积MATLAB程序 附录4 自适应辛普森求积MATLAB程序

二维偏微分方程微分求积法matlab代码，仅供参考！！！