/*
 * Estimate machine epsilon
 */

#include <stdio.h>

int main() {
    double eps = 1.;

    while ( eps + 1. > 1.) {
        eps /= 2.;
    } 

    eps *= 2;

    printf("machine eps: %g\n", eps);

    return(0);
}