visual fortran常用数值算法集_从零开始的科研之路：数值积分

感谢师兄提供的题图。

（想学习代码的，直接看代码即可，文字大部分都是故事回忆。）

（重申，对数值积分感兴趣的，建议去阅读相关文献或教材，本文并不重点涉及讨论这些算法。）

Z师兄是个动漫迷，老二次元，刚进实验室的新人一般都会被他拿这个当测量工具，测一测你是不是同好。我刚进去的时候，没看过几部动漫，看过的都是热血漫，也没聊到一块去。后来混熟了发现，这个人看的动漫都是卖肉的（大胸妹子），卖肉这个词还是他教我的。很难想象，这个人从小到大都是看国学的。

刚开始科研那会儿，要算点场论的东西，从很简单的平均场开始，需要用到一些积分。当时我们组做的事比较新的方向，所以物理理论方面一般以新为主，不追求精确和深度，平均场就够用了。

虽然我本科学过数值方法这门课，讲过数值积分的方法，但是我都忘得差不多，想起来只是记得矩形积分。基本思想就是，积分就是算曲线下的面积，将面积细分，每一块用矩形近似，减小分的区间（增加分的次数），就可以逼近真正的面积。

实现如下：我们以