2008-03-21から1日間の記事一覧
ちょっと計算したらすごく簡単な式になった。 これを基に再実装。 open Num let pi length = let min = Int 1 // Int 10 **/ Int length in let zero = Int 0 and one = Int 1 and two = Int 2 and three = Int 3 and four = Int 4 in let rec sum n an x = …
なんとなくやりたくなったからやってみた。 意外と素直にあっさり書けた。やっぱり有理数型があるとこういうの楽だなぁ。 使ったアルゴリズムは、円周率 - Wikipediaで見つけたこの式。 arcsin(x)のx=1/2におけるテイラー展開で、収束が速いらしい。へぇー。…