如何翻译使用超几何的C ++程序 来自GNU Scientific Library
的分发函数http://www.gnu.org/software/gsl/manual/html_node/The-Hypergeometric-Distribution.html
进入一个C ++程序,而是使用Boost library中的类似函数?
#include <cstdlib>
#include <iostream>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_cdf.h>
int main(int argc,char *argv[]) {
unsigned int n1 = 256;
unsigned int n2 = 1583;
unsigned int t = 300;
unsigned int k = 40;
std::cout << gsl_ran_hypergeometric_pdf(k, n1, n2, t)
<< std::endl
<< gsl_cdf_hypergeometric_P(k, n1, n2, t)
<< std::endl;
return EXIT_SUCCESS;
}
答案 0 :(得分:1)
以下是Boost功能的翻译:
#include <cstdlib>
#include <iostream>
#include <boost/math/distributions/hypergeometric.hpp>
#include <boost/math/policies/policy.hpp>
int main(int argc, char *argv[]) {
unsigned int n1 = 256;
unsigned int n2 = 1583;
unsigned int t = 300;
unsigned int k = 40;
boost::math::hypergeometric_distribution<double> hg_dist(n1, t, n1 + n2);
std::cout << boost::math::pdf<double>(hg_dist, k)
<< std::endl
<< boost::math::cdf<double>(hg_dist, k)
<< std::endl;
return EXIT_SUCCESS;
}