では次に，n の階乗 n!, n=1, 2, ... 100000 の先頭桁の数字の分布はどのようになるか？ああ，簡単ですね。前のプログラムをまねて，っと...import scipy as spn = 100000tbl = sp.zeros(10, dtype=int)f = 1for i in range(1, n+1): f *= i tbl[int(str(f)[…

フィボナッチ数列の各項の先頭の数字の度数分布はどのようになると思うか？ そんなの簡単。どの数字も同じ頻度。 ブッブー。違います。1が一番多く，順に少なくなり，9が一番少ない。 ベンフォードの法則という。 Python だと，int 変数は自動的に精度が確保…

R で，実数の内部表現をプリントさせてみると，4項目から既に，末尾にゴミが入る。 これは，多倍長演算をしても，例えビット数を無限にしても最下位にゴミが入る（だって，2の羃乗数でない限り2進数では正確に表せないのだから，あたりまエ） 第 4 項> a[1] …

「フィボナッチ数列の一般項の実数演算結果がExcelとLibreOffice Calcで違う」ということだ。 まあ，そうなんだろうけど 議論の展開とは別だが，実際問題として，結果は整数になって欲しいんだから，演算誤差があろうと，整数に丸めてやれば何の問題もない。…

件の本山だけど，標高も間違えていた 敵性の Google なんか使うからだ。わが国には国土地理院があるぞ。本山の標高は 715 メートルだ。

> 遮蔽（しゃへい）となる尾根などがあった場合はそこの標高、ない場合は山頂の標高、さらに国有地までの水平距離を紙の上で定規を使って計測。高さと水平距離の縮尺の違いに気付かないまま三角関数を用いて計算し、誤った仰角を算出した ああ，一応は三角関…

調査ミスではなく，初歩的な計算ミス指摘する新聞なども，基本的なことが分かってるのか？という懸念がある。 標高 a メートルの A 地点から，標高 b メートルの B 山の山頂を見るとき，その仰角はいくつか？ただし，A，B 間の水平距離は x メートルとする。…

> set.seed(1234567)> n = 200> x = rnorm(n)> y = rnorm(n) n が50 以上のとき，exact を指定しない場合は正規近似が行われる> cor.test(x, y, method="kendall") Kendall's rank correlation taudata: x and yz = 1.2089, p-value = 0.2267alternative hyp…