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

double          d1mach(long i)
{
    switch (i) {
    case 1:
        return DBL_MIN;
    case 2:
        return DBL_MAX;
    case 3:
        return DBL_EPSILON / FLT_RADIX;
    case 4:
        return DBL_EPSILON;
    case 5:
        return log10((double) FLT_RADIX);
    default:
        break;
    }
    fprintf(stderr, "invalid argument: d1mach(%ld)\n", i);
    exit(EXIT_FAILURE);
}