fork download
  1. PROGRAM est_e
  2.  
  3. IMPLICIT NONE
  4.  
  5. INTEGER n, lcv
  6. DOUBLE PRECISION este, reldif, actdif, factor, actval
  7.  
  8.  
  9.  
  10. actval = 2.718281828D0
  11. n=1000
  12. factor = 1
  13. este = 1
  14. DO lcv = 1, n, 1
  15. factor = factor * lcv
  16. este = este + 1/factor
  17. END DO
  18.  
  19. actdif = actval - este
  20. reldif = actdif/actval
  21.  
  22.  
  23. print *, 'The estimate of e is: ', este
  24. print *, 'Actual difference is:', actdif
  25. print *, 'Relative difference is:', reldif
  26.  
  27.  
  28. END PROGRAM est_e
Success #stdin #stdout 0s 3984KB
stdin
Standard input is empty
stdout
 The estimate of e is:    2.7182818284590455     
 Actual difference is:  -4.5904569034860288E-010
 Relative difference is:  -1.6887347206612122E-010