統計数値表編集委員会: 簡約統計数値表. 日本規格協会,東京,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
}
}
}