Fortran 95二重合計が正しく出力されない

Question

Fortran 95で二重合計を書こうとしていますが、正しい出力を得られないようです。常にFortranの非常に小さな実数を出力するようです。私のコードは以下の通りです：

program summation 
    implicit none 
    integer::i,j,m,n 
    real::s,s1 
    s1=0 
     do i=1,100 
     do j=1,100 
      s1=s1+(1/((i**2)+(j**2))) 
     end do 
     s=s1+s 
     end do 
    print*, s 
    end program summation 

sの最終値は約6.4でなければなりません。

Answer 1

最初にsを設定したことはありません。だからあなたは完全なゴミを得る。

s = s + s1を作成した後、すべての繰り返しでs1を0に再設定する必要があります。

また、1を大きな整数で割ると、結果は常に整数0になります。あなたはtye用語の1つを作成する必要があります。real。私の例では、1.を使って1を実際にします。

ので

s=0 

do i=1,100 
s1=0 
do j=1,100 
    s1=s1+(1./((i**2)+(j**2))) 
end do 
s=s1+s 
end do 

print*, s 

可能ですが、なぜちょうど1つのアキュムレータにすべてを合計していませんか？

s=0 

do i=1,100 
do j=1,100 
    s=s+(1./((i**2)+(j**2))) 
end do 
end do 

print*, s