| 32 | double factorial(int n) { return std::tgamma(n + 1); } |
| 33 | |
| 34 | double Binomial(int n, int i) { |
| 35 | return factorial(n) / (factorial(i) * factorial(n - i)); |
| 36 | } |
| 37 | |
| 38 | double Bernstein(double x, int i, int n) { |
| 39 | return Binomial(n, i) * pow(x, i) * pow(1 - x, n - i); |