fortranの四則演算

Question

標準偏差を求めるプログラムを作りたいのですが、ルートの計算の出力が「NaN」とでます。
解決方法を教えてください。
以下がプログラムです
c234567
dimension a(100)
*
read(,) n
do i = 1, n
read(,) a(i)
end do
*
sum = 0.0
sum2 = 0.0
do i = 1, n
sum = sum + a(i)
sum2 = sum2 + a(i)a(i)
write(,) i, a(i), sum, sum2
end do
*
ave = sum/real(n)
sig = sqrt(1/nsum2-sumsum)
*
write(,) 'sum =' , sum
write(,) 'ave =', ave
write(,*) 'sig =', sig
*
stop
end

Answer 1

NaNになるのは、sqrt(1/n*sum2-sum*sum)の1/nが整数除算なのでゼロになり、sqrtの引き数が負になるからです。

そもそも、この部分は計算式がおかしいけど。

Comments

[kawasann]

回答ありがとうございます。nを実数にすることは分かったのですが、正しい計算式がわからないです。教えてもらえませんか？

[otn]

標準偏差： https://ja.wikipedia.org/wiki/%E6%A8%99%E6%BA%96%E5%81%8F%E5%B7%AE

[kawasann]

解決しました！