Fortran epsilon test

gistfile1.txt

program test
implicit none
integer, parameter :: dp = 8
real(dp) :: a, b
a = 1.0_dp
b = a + epsilon(1.0_dp)
print *, "First we show, that we have two different 'a' and 'b':"
print *, "a == b:", a == b, "b-a:", b-a
print *, "using (es22.15)"
print "(es22.15)", a
print "(es22.15)", b
print *, "using (es23.16)"
print "(es23.16)", a
print "(es23.16)", b
end program