作者sjgau (sjgau)
看板Fortran
標題Re: [問題] 關於誤差估計的問題
時間Thu May 20 11:25:58 2010
換成 real*8 以後，
(1.0D0 + 1.0D-8) 的結果，可以被計算出來，
有效位數，可以達到
2^52 --> 大概是 15位數
大家可以試試看
program VF1035
implicit none
real*8 x, y, z, ct1
x= 4.0D0*datan(1.0D0)
ct1= 0.0D0
print *, 'ct1, x= ', ct1, x
! pause
y= x + 1.0D0
! y > x
do while (y .GT. x)
x= x*(1.0D0 + 1.0D-8)
y= x + 1.0D0
ct1= ct1 + 1.0D0
end do
! y <= x
print *, 'x, y, ct1= ', x, y, ct1
end program VF1035
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.20.174.109