/*
 * Calculate I_n = \int_0^1 e^{x-1} x^n dx
 *
 * integrating by parts, I_n = 1 - n*I_{n-1}
 */

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

#define J 30

int main(void) {
    double i;
    int j;

    i = 1 - 1/exp(1.);  /* I_0 */

    for (j = 1; j <= J; j++) {
        printf("%3d    %25.6f\n", j-1, i);
        i = 1. - j*i;   /* I_n = 1 - n*I_{n-1} */
    }
    
    return(0);
}