How many significant decimal digits for real single precision?

@certik suggested print *, x, y for checking the values of x and y. Different compilers choose different numbers of digits to print (the Fortran standard allows that), and so I recommend using an explicit format such as '(F0.8)' or even '(F0.9)' etc. when one is investigating the last digit.