/** Machine epsilon eps.c **/

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

int main ( )
{
	double ex, g=1, eps;

	do {
		g = g/2;
		ex = g*0.98 + 1;
		ex = ex - 1;
		printf("g=%e  ex=%e\n", g, ex );
		if ( ex>0 ) eps = ex;
	} while ( ex>0 );
	printf("\nMachine epsilon = %e %e\n", eps, sqrt(eps) );
}