因为,integer类型默认为四字节整数,最大值为2^31-1=2147483647。而13的阶乘为6227020800,已经超出integer上限。可以改用real*8双精度类型,这样n=33就不会溢出了。双精度类型的上限是10的308次方,最大可以计算到n=170,即170的阶乘。
http://www.fcode.cn/algorithm-50-1.html我也不太理解,我是问的大神,论坛那有专门的讲解,你看看吧