ベキ分布乱数: プログラム: 篠田孝祐ベキ分布乱数という単語が存在するかはわからないが，任意の確率密度関数に従った乱数をほしいと思うことはある．今回は，ベキ分布に従った乱数を生成する関数がないかと探したが見当たらなかったので，なんとなく作成した．それっぽいデータは生成できるが，正確とは限らないのであしからず． ガウス分布の用に計算で出力できたら良いのだが，アルゴリズムは見つからない． 指数ベキ分布の密度関数はめんどくさそうなので，Zipf-Mandelbrot lawに従う密度分布を以下の関数で作成． public static double powerlaw(double x, int N, double q, double s) { double h = 0; for(int i = 1; i <= N; i++) { h += 1/Math.pow(i + q, s); } return (1/Math.pow(x +
