X-Git-Url: https://www.hackdaworld.org/gitweb/?p=physik%2Fcomputational_physics.git;a=blobdiff_plain;f=integral-1_2.c;fp=integral-1_2.c;h=e49100b21987a1111e2c4d4b5715f4d386dad6b8;hp=0000000000000000000000000000000000000000;hb=e0f00c13eac7c9f9ae735bb7928057fc1445b056;hpb=48ea798c247dfff5f6bcb8c6b51b115e829571f8

diff --git a/integral-1_2.c b/integral-1_2.c
new file mode 100644
index 0000000..e49100b
--- /dev/null
+++ b/integral-1_2.c
@@ -0,0 +1,28 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+
+int main(int argc,char **argv) {
+ double p_N;
+ double p;
+ int start;
+ int i,j;
+
+ if(argc!=3) { 
+  printf("usage: %s <start> <startwert>\n",argv[0]);
+  return 1;
+ }
+
+ start=atoi(argv[1]);
+ p_N=atof(argv[2]);
+
+ printf("debug: N=%d startwert=%f\n",start,p_N);
+
+ for(i=start;i>=0;i--) {
+  p=(1/((double)i+1))*(M_E-p_N);
+  printf("p_%d = %f\n",i,p);
+  p_N=p;
+ }
+
+  return 1;
+}