fortran 代码中的循环不遵守循环规则

Question

首先，我是 fortran 的新手。我试图在 0.0000005 的容差范围内找到 1/e 的值。我使用了 1/e 的求和表示，即 (epsilon) n 从零到无穷大 ((-1)^n)/n!。我从n=2开始，当值为1/n！小于我的容差，程序将停止并打印计算出的总值。但是我的程序只进入 n=3 并且只打印 1/3 的值！这是 1.666666。 !程序被编辑。编辑后的表格是计算我想要的。在外部while的条件是（num3

   program Ecalculator
   implicit none

   integer :: mult,num1,n,num4,num5
   real :: summ,num3,fact

   mult=1
   n=2.0
   fact=1.0
   summ=0.0

      DO WHiLE(n<13)
      fact=n
      num1=n-1

             DO WHiLE(num1>0)
             fact=fact*num1;
             num1=num1-1;

             END DO

       fact=fact*mult;
       num3=1.0/fact;
       mult=mult*(-1);
       summ=summ+num3;
       n=n+1;

        END DO

    print *, summ
    read  *, num5

    end program Ecalculator

Answer 1

看起来阶乘计算不正确：fact=num4*num1 行可能没有执行您想要的操作，因为 num4 只是分配为 n 并且在整个内部循环中永远不会更改。我认为您不需要num4 和fact 这两个变量；您可以将它们组合成一个变量。