e49100b21987a1111e2c4d4b5715f4d386dad6b8
[physik/computational_physics.git] / integral-1_2.c
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <math.h>
4
5 int main(int argc,char **argv) {
6  double p_N;
7  double p;
8  int start;
9  int i,j;
10
11  if(argc!=3) { 
12   printf("usage: %s <start> <startwert>\n",argv[0]);
13   return 1;
14  }
15
16  start=atoi(argv[1]);
17  p_N=atof(argv[2]);
18
19  printf("debug: N=%d startwert=%f\n",start,p_N);
20
21  for(i=start;i>=0;i--) {
22   p=(1/((double)i+1))*(M_E-p_N);
23   printf("p_%d = %f\n",i,p);
24   p_N=p;
25  }
26
27   return 1;
28 }