Using four-dimensional symplectic maps as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to the single resonances. We measure an exponential dependence of the size of the lobes of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by
the Nekhoroshev theorem. The variation of the size of the lobes turns out to be correlated to the diffusion coefficient