Fortran 的選擇式結構

這是 Fortran 90 所提供的 選擇性的 程式結構，
很重要。

當 係數 a,b,c= 1, 3, 1 的時候，驗算結果會產生 誤差，
這就是 浮點數的運算，一定會 產生誤差。

如何 正確處理 誤差，是一個很重要的 學問。

评论

1 楼 sjgau 2007-01-31  

! main()    
implicit none    
real a,b,c, d, x1,y1, x2,y2    
   
! y= a*(x^2) + b*x + c    
! D= (b^2) - 4*a*c    
! -----------------------------------------------    
   
a= 1.0    
b= 3.0    
! b= 2.0    
! b= 1.0    
   
c= 1.0    
! 修改上面的 a,b,c 係數的值，重新 compile    
! -----------------------------------------------    
   
d= (b*b) - 4.0*a*c    
if (d > 0.0) then    
    print *, 'D > 0.0, D= ', d    
    pause    
   
    x1= (-b + sqrt(d))/(2.0*a)    
    x2= (-b - sqrt(d))/(2.0*a)    
   
    ! 驗算    
    y1= a*(x1*x1) + b*(x1) + c    
    y2= a*(x2*x2) + b*(x2) + c    
    print *, 'x1= ', x1, ', y1= ', y1    
    print *, 'x2= ', x2, ', y2= ', y2    
else if (d < 0.0) then    
    print *, 'D < 0.0, D= ', d    
    pause    
else   
    print *, 'D = 0.0, D= ', d    
    pause    
   
    x1= (-b + sqrt(d))/(2.0*a)    
   
    ! 驗算    
    y1= a*(x1*x1) + b*(x1) + c    
    print *, 'x1= ', x1, ', y1= ', y1    
end if   
   
end !of main()    
! -----------------------------------------------