def brouncker(n): S = 0 for k in range(n): S += 1/(2*k+1)/(2*k+2) return S