У меня проблемы с реализацией Gammainc of Matlab в Java. < /p>
Я попытался сделать это частично. Я использовал эту функцию, которая использует аппроксимацию Lanczos, которую я получил в Интернете для решения γ (a) < /strong>: < /p>
private static double logGamma(double x) {
double tmp = (x - 0.5) * Math.log(x + 4.5) - (x + 4.5);
double ser = 1.0 + 76.18009173 / (x + 0) - 86.50532033 / (x + 1)
+ 24.01409822 / (x + 2) - 1.231739516 / (x + 3)
+ 0.00120858003 / (x + 4) - 0.00000536382 / (x + 5);
return tmp + Math.log(ser * Math.sqrt(2 * Math.PI));
}
private static double gamma(double x) { return Math.exp(logGamma(x)); }
< /code>
Затем я использовал правило Симпсона при решении интегральной части, а затем я объединил их, чтобы сделать Gammainc, но вывод, который я получаю, не является разумным. < /p>
Интегральная часть также может рассматриваться как более низкая функция гаммы.>
Подробнее здесь: https://stackoverflow.com/questions/870 ... te-in-java
