高斯勒让德求解复杂函数问题

pasuka

kerb 发表于 2016-8-17 21:07
http://bbs.fcode.cn/forum.php?mod=viewthread&tid=814&page=1&extra=#pid5135
建议参考这个帖子，
如果 ...

lz利用多项式求根方法，是否比较过60阶之后的结果？印象中十多年前有老外在mathworks网站上贴出对比，多项式求根方法阶次过60后，积分权重出现负值。一般而言，低于10阶，查手册，高的话构造伴随矩阵算特征值，类似的matlab代码很多，改回fortran不难。*************************************************************************************************************************
补充一点干货
ACM的CALGO下面的Fortran现成代码
http://people.sc.fsu.edu/~jburkardt/f_src/toms655/toms655.html
某老外几年前的代码且附带高斯积分点结果
http://www.holoborodko.com/pavel ... erical-integration/

pasuka

kerb 发表于 2016-8-23 22:05
我不明白你的意思

多项式直接二分法求根的办法，稳定性随着阶次增加会下降
构造伴随矩阵求解能够避免上述问题，这个是计算方法课程的基本内容。。。

pasuka

算法研究版面不讨论算法的稳定性、收敛性、鲁棒性等等有趣的问题，那就不好玩了。。。

pasuka

kerb 发表于 2016-8-24 12:41
说句刺激一点的话，叫做自以为是！

原来不是lz本尊啊~~~
灌水灌多了。。。。