对同一个常数，函数前后输出值不同

不知道 · 发表于 2016-10-20 21:17:05

主程序中首尾两个is_prime(32)输出值不同，这是怎么回事？
把do循环删除后又变得正常了。

[Fortran] 纯文本查看 复制代码

01program test
02  implicit none
03  logical,external :: sum_is_5,is_prime
04  integer :: n
05  print *,is_prime(32)
06  do n=2,100
07    if(sum_is_5(n).and.is_prime(n)) print *,n
08  end do
09  print *,is_prime(32)
10end program test
11logical function sum_is_5(i)
12  implicit none
13  integer :: i,digit(3)=0
14  sum_is_5=.false.
15  digit(1)=mod(i,10)
16  digit(2)=mod(i/10,10)
17  digit(3)=i/100
18  if(sum(digit)==5) sum_is_5=.true.
19end function sum_is_5
20logical function is_prime(j)
21  implicit none
22  integer :: j,k=2
23  is_prime=.true.
24  do while(k<j)
25    if(mod(j,k)==0) then
26      is_prime=.false.
27      exit
28    end if
29    k=k+1
30  end do
31end function is_prime