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

//
//   estimate machine epsilon
//   For single precision, expect eps = 2^-24

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

  eps = 1.0;
  n=30;

  printf("Largest 32-bit floating point number %E \n",FLT_MAX);
  printf("Smallest 32-bit (normalized) floating point number %E \n",FLT_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%23.18lf \n", i,  eps, res, res2);
    eps=eps*0.5;
  }

  return 0;
}