統計数値表編集委員会: 簡約統計数値表. 日本規格協会,東京,1984. にも,n = 50 までの値は出ていますが,代表的な有意水準に対応する統計量という形です。
以下のスクリプトでは,あっという間に結果が出るでしょう。
別の所に示した u_table.awk は,この延長線上にあった。いわば,mp_table.awk が1次元とすると,u_table.awk は2次元的な広がりを持つということ。
# mp_table.awk BEGIN { LIMIT = 50 denom = table[0] = 1 sum = 0 for (n = 1; n <= LIMIT; sum += n++) { for (i = sum; i >= 0; i--) { table[i+n] += table[i] } denom /= 2 lim = (n+1)*n/2+1 cum2 = 1 cum1 = 0 printf "\n\nn = %i\n\n", n printf "%4s:%15s %13s%13s%13s\n", "S", "度数", "確率", "累積1", "累積2" for (i = 0; i < lim; i++) { x = table[i] y = x*denom cum1 += y printf "%4i:%15.15g %13.10f%13.10f%13.10f\n", i, x, y, cum1, cum2 cum2 -= y } } }