显示数组溢出，实在找不到问题了，请大家帮忙

kif117

我觉得问题也许是在subroutine getgama中。一般我的经验是按步骤测试，首先去掉sub，运行主程序，然后一步一步加入call。加入getgamma后，print*, gama 没有正确显示数值(出现NaN)。

纯属个人意见，仅供参考。