gam2k = np.where(abs(c) < 1e-7, np.pi**2.0/6.0,