#include <iostream>
#include <boost/math/distributions.hpp>
using namespace std;
using namespace boost::math;
int main()
{
// 二項分布の例 (40% の確率で表になるコインを 5 回投げた場合の確率分布)
binomial b( 5, 0.4 );
// 確率密度関数 (ちょうど x 回表が出る確率)
for(int x=0; x<=5; ++x)
cout << x << ": " << pdf(b, x) << endl;
cout << "----" << endl;
// 累積分布関数 (x 回以下表が出る確率)
for(int x=0; x<=5; ++x)
cout << x << ": " << cdf(b, x) << endl;
}
0: 0.07776 1: 0.2592 2: 0.3456 3: 0.2304 4: 0.0768 5: 0.01024 ---- 0: 0.07776 1: 0.33696 2: 0.68256 3: 0.91296 4: 0.98976 5: 1
確率分布・統計分布に関する計算あれこれ。