g2 = np.pi**4/15.0 + 6*sc.zeta(4, c)