/*
 * Calculate integral \int_0^1 x^n e^(-x) dx
 */
#include <stdio.h>
#include <math.h>

int main() {
    int n;
    double I;

    n = 0;
    I = 1. - 1./M_E;
    printf("%3i   % 24.14f\n", n, I);
    
    for(n = 1; n < 24; n++) {
        I = -1./M_E + n*I;
        printf("%3d   % 24.14f\n", n, I);
    }

    return(0);
}