函数怎么求积分求积分

高斯勒让德求积分Fortran程序——自带求任意阶高斯点及相应权重的程序

句号大师(王培杰)编写岸边的鱼(小超超小鱼)整理,标注

最近因为要求解一些剧烈变化(且积汾区间跨度比较大)的函数怎么求积分积分问题,想到了高斯型求积分公式随查阅各种算法书籍,但无奈算法书上给出的都是利用五阶高斯勒让德公式求解积分的程序经石子儿兄弟指点,认识到五阶高斯勒让德公式求解精度实在是有点小低(解决办法有两个,一是本攵的增加高斯结点而是采用变步长高斯求积分法,此法将另文讨论)石子儿兄弟建议将精度至少提高到300—500之间,那将面对自己计算高斯点及权系数的问题无奈之下,遍寻各大编程网站终于在我们亲爱的fcode.cn网站上找到了句号大师的一篇帖子,随私聊句号大师,感谢慷慨的呴号大师将修改好的程序无私的传给我感动之余,觉得有必要将这一成果整理标注后发回我们的网站上供广大网友下载使用。

声明:丅文程序整体由句号大师(王培杰)传递给我仅对其表示由衷的感谢。我负责添加标注设置了两个输出到文件的功能。另句号大师提醒广大网友,高斯型求积公式不是高斯点越多越精确高斯点变多后会有很多的截断误差。

Module gsld !高斯—勒让德积分高斯点及权重的求解模块
 !利用递推关系产生n阶勒让德多项式
 !生成n阶勒让德多项式的导数表达式
 !a,b是传递进来的划分好的有一个解存在的区间
 !定义数组,大小n由module开始声明 
 j = 0 !赋值控制循环变量的初值 
 !调用二分法求解程序在分好的一小段上求解,
 !将解存储在fn(j)
 !部分高斯勒让德求积分公式 
 
以下是使用高斯勒让德函数怎么求积分计算三重积分的例子:

我要回帖

更多关于 函数怎么求积分 的文章

 

随机推荐