#include <stdio.h>
#include <math.h>
#include <float.h>

//
//   estimate machine epsilon
//   For double precision, expect eps = 2^-53

int main()
{
  double eps, res, res2;
  int i, n;

  eps = 1.0;
  n=60;

  printf("Largest 64-bit floating point number %E \n",DBL_MAX);
  printf("Smallest 64-bit (normalized) floating point number %E \n",DBL_MIN);

  for(i=0; i<n; i++)
  {
    // remember to use lf ("long float") if your data are double
    res = (1.0 + eps) ;
    res2 = res -1.0;
    printf("%i \t%E \t%23.18lf \t%E \n", i,  eps, res, res2);
    eps=eps*0.5;
  }

  return 0;
}