ビュフォンの針でとおもったんだけど、なぜかarcsinのテーラー展開が目に留まったので。

def factorial(i):
   return 1 if i ==0 else i * factorial(i-1)

def calc_pi(i):
   pi = 0
   for n in range(i+1):
       pi += 3.0*factorial(2*n)/(2**(4*n)*((factorial(n))**2)*(2*n+1))
       print pi
   return pi


if __name__ == "__main__":
   print calc_pi(16)

「無限と連続」の数学―微分積分学の基礎理論案内
瀬山 士郎
東京図書 / ￥ 2,625 ()
在庫あり。