复合辛普森求积公式、龙贝格方法、复合三点高斯-勒让德公式的实例应用编程

作者: Jonathan M. Borwein / Peter B. Borwein 出版社: Wiley-Interscience 副标题: A Study in Analytic Number Theory and Computational Complexity 出版年: 1998-06-29 页数: 432 定价: USD 137.95 装帧: Paperback ISBN: 9780471315155

函数公式说明： %%高斯型求积公式 %%Y是函数表达式，interval是求积区间，n是求积阶数 %%对于求一般形式的非反常积分，可用勒让德型， %%对于求形如f(x)/sqrt(1-x^2)的非反常积分，可用第一类切比雪夫型， %对于形如f(x)*sqrt(1-x^2)的非反常积分，可用第二类切比雪夫型，切比雪夫型积分应在[-1 1]上 %对于反常积分f(x)*exp(-x)或者f(x)*exp(-x^2)，且区间为[0,+inf]或[-inf,+inf]，可用拉盖尔或者埃尔米特型（注意Y是f(x)） %%正交多项式 %此函数包括勒让德正交多项式（定义区间[-1,1]），切比雪夫正交多项式（两类， %在这里，规定第一类切比雪夫多项式是以1/sqrt(1-x^2)作为权函数，第二类切比雪夫多项式以sqrt(1-x^2)作为权函数得到的）（定义区间[-1,1]） %拉盖尔正交多项式（定义区间[0 +inf]），埃尔米特正交多项式（定义区间[-inf +inf]），输入项数N应从1开始 %%n是多项式的项数，n>=0,type是类型，分为Legendre、Chebyshev1、等

高斯一维数值积分，可以任意均分积分区间，每个积分区间的高斯点为10.采用双精度浮点数。自己编的

我相信找到这里来的都已经知道啥玩意是高斯-勒让德积分了，这里就不哔哔了。 作为一个工科狗，离了matlab几乎不能生活，这可不行，要搞UI还是C++的好。 缺算法？所以我开始扒拉了。 从MATLAB quadgk抄过来的，说是抄，其实区别还是比较大的，核心计算基本上是重新排布了，毕竟人matlab矩阵计算多方便啊。烦就烦在以前找不到可靠的现成C++算法，从头写起我还得学一遍。。。 不多说，这玩意儿支持实数double的gk自适应积分，支持complex的gk自适应积分，算法和结果与matlab一模一样。有些删减和调整，你们自己看着办吧。 有需要进一步开发的话，网站联系我就好，我一周至少上线一次吧。

高斯-勒让德积分（Gauss-legendre积分）matlab程序 可直接在函数名前设置输出为 [x,w]= 其中x为积分点，w为权重

高斯—勒让德积分公式 (带有程序)，非常好的资源奥