#include <stdio.h>
#include <math.h>
int main(int argc, char **argv)
{
printf("expm1(-1) = %.6f\n", expm1(-1));
printf("expm1(1) = %.6f\n", expm1(1));
printf("expm1(0) = %.6f\n", expm1(0));
printf("expm1(-0) = %.6f\n", expm1(-0));
printf("expm1(5) = %.6f\n", expm1(5));
printf("expm1(0.5) = %.6f\n", expm1(0.5));
printf("expm1(-INFINITY) = %f\n", expm1(-INFINITY));
printf("expm1(100) = %f\n", expm1(100));
printf("expm1(710) = %f\n", expm1(710));
return 0;
}
// 1.#INF00 - inf - Numerical result out of range
/*
run:
expm1(-1) = -0.632121
expm1(1) = 1.718282
expm1(0) = 0.000000
expm1(-0) = 0.000000
expm1(5) = 147.413159
expm1(0.5) = 0.648721
expm1(-INFINITY) = -1.000000
expm1(100) = 26881171418161356000000000000000000000000000.000000
expm1(710) = 1.#INF00
*/