#include void printReal(REAL x) { if (sizeof(REAL) == sizeof(double)) printf("%22.16e",x); else printf("%15.9e",x); } int main() { REAL m=1.0, s=0.0, mf=1.0; int n=1; while (s+mf != s) { s += mf; m *= (REAL) n; mf = 1.0/m; printf("%d ",n); printReal(s); printf("\n"); n++; } }