added integral-2_2.c aufwaerts rekursion (0...N)
[physik/computational_physics.git] / integral-2_2.c
diff --git a/integral-2_2.c b/integral-2_2.c
new file mode 100644 (file)
index 0000000..2c8a317
--- /dev/null
@@ -0,0 +1,26 @@
+#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;
+}