Generating function identities for ζ(2n+2), zeta(2n+3) via the WZ method