+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+
+int main(int argc,char **argv) {
+ double p_N;
+ double p;
+ int steps;
+ int i,j;
+
+ if(argc!=2) {
+ printf("usage: %s <steps>\n",argv[0]);
+ return 1;
+ }
+
+ steps=atoi(argv[1]);
+ p_N=M_E-1;
+
+ for(i=0;i<=steps;i++) {
+ printf("p_%d = %f\n",i,p_N);
+ p=M_E-i*p_N;
+ p_N=p;
+ }
+
+ return 1;
+}