新人求教，程序当中的内循环都只执行了一次，不懂，求解

Rong · 发表于 2019-10-20 09:25:58

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

01MODULE  my_Constants
02implicit none
03real(8),PARAMETER :: df=1         
04real(8),PARAMETER :: fmax=196078   
05real(8),PARAMETER :: f00=191693   
06real(8),PARAMETER :: f01=191694
07real(8),PARAMETER :: f02=191696
08end MODULE my_Constants
09use my_Constants
10implicit none
11integer  i,j,k,l,m,n,s,z,a,b,c,d
12real(8),dimension(4385)::fijk,f
13open(11,file='useful.dat')
14do i=1,4385
15  f(i)=0.0
16enddo
17    f(1)=f00
18    f(2)=f01
19    f(3)=f02
20    n=4
21    z=1
22    write(11,*)f(1),z
23    write(11,*)f(2),z
24    write(11,*)f(3),z
25  do 10 i=4,4385           
26       m=n   
27   f(m)=f00+i*df
28    do  j=1,m
29        do   k=1,m
30            do  l=1,m
31                fijk(l)=f(j)+f(k)-f(l)
32                do 50 s=1,m
33                if((fijk(l).ne.f(s)).and.(fijk(l).le.fmax).and.(fijk(l).ge.f00)) then
34                write(11,*)f(m),z
35                n=m+1
36                goto 10
37                endif
3850              continue
39               enddo
40            enddo
41        enddo
4210    continue
43    end