- 599 名前：デフォルトの名無しさん [2007/06/15(金) 14:20:39 ]
- 教えて下さい。
H(l)を以下のように計算したいのですが、走らせても答えが「NaN」と表示されます。
どこがおかしいか教えて下さい。お願いします。
program FILTER
c
implicit real*8(a-h,o-z)
dimension H(100000)
open(115,file='H.dat')
c
do 10 l=1,n/2+1
f=df*(l-1)
e=f/10.0
H(l)=((1.0/f)**0.5)*((1.0+(0.694*(e**2))+(0.241*(e**4))
&+(0.0557*(e**6))+(0.009664*(e**8))+(0.00134*(e**10))
&+(0.000155*(e**12)))**(-0.5))*((1.0-(dexp(-(f/0.5)**3)))**0.5)
write(115,*) H(l)
10 continue
c
pause
stop
end
（注）＆は数式が長くなったため、改行している事を意味しているようです。
