main.f90 209 B

123456789101112131415
  1. program main
  2. use m1, only : pi
  3. use m4, only : halfpi
  4. implicit none
  5. real :: rpi, rhalfpi
  6. rpi = pi() / 2
  7. rhalfpi = halfpi()
  8. print '(a,ES15.8)', 'floating point precision loss: ', rpi - rhalfpi
  9. end program