blob: 7a8680046f8123319a47776202255bf4e21caa3c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
/*
* l.cpp
*
* (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
*/
#include <cstdlib>
#include <cmath>
#include <cstdio>
int main(int argc, char *argv[]) {
double x = 0.5;
double g = tgamma(x);
printf("limit: %20.16f\n", g);
double p = 1;
long long N = 100000000000;
long long n = 10;
for (long long k = 1; k <= N; k++) {
p = p * k / (x + k - 1);
if (0 == k % n) {
double gval = p * pow(k, x-1);
printf("%12ld %20.16f %20.16f\n", k, gval, gval - g);
n = n * 10;
}
}
return EXIT_SUCCESS;
}
|
