program subtract
  implicit none
  real :: a = 1.0002
  real :: b = 1.0001
  real :: c

  c = a - b
  print *, a
  print *, b
  print *, c
end program subtract