blob: dab8d238da60e7fa7c1600ae6e36a5180e046e0a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
|
# cython: profile=True
# calc_pi.pyx
def recip_square(int i):
return 1. / i ** 2
def approx_pi(int n=10000000):
cdef double val = 0.
cdef int k
for k in range(1, n + 1):
val += recip_square(k)
return (6 * val) ** .5
|