#include <stdio.h>
#include <math.h>
#define EPSILON 1.0e-8
int main() {
double sum, term;
int p, i;
p = 1;
sum = 1.0;
for (i = 1;; i++) {
/* break-condition */
if ((term = 1.0 / (p = p * i)) < EPSILON)
break;
sum += term;
}
printf(" +0.1234567\n");
printf("result: %+.12f\n", sum);
printf("exp(1) = %+.12f\n", exp(1.0));
printf("diff: %+.12f\n", sum - exp(1.0));
return 0;
}
/* end */