return np.pi**2 / 2 + (np.log(2) + np.euler_gamma)**2