Python profiling

profiling.py

import pstats, cProfile

def recip_square(i):
    return 1./i**2

def approx_pi(n=10000000):
    val = 0.
    for k in range(1,n+1):
        val += recip_square(k)
    return (6 * val)**.5
 
if __name__ == '__main__':
  cProfile.runctx("approx_pi()", globals(), locals(), "Profile.prof")

  s = pstats.Stats("Profile.prof")
  s.strip_dirs().sort_stats("time").print_stats()