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におけるテイラー展開で、収束が速いらしい。へぇー。…