non_parallel_pi.py の履歴(No.1)

  1
  2
  3
  4
  5
  6
  7
  8
  9

n=10000
d=1.0/n
d2=d/2.0
s=0.0
for i in range(0,n):
    x=i*d+d2
    s=s+d/(x*x+1.0)
s=s*4.0
print(s)