A simple analytic method for calculating nucleon-nucleus optical potentials using a single folding of a Gaussian two body interaction with an arbitrary nuclear distribution is presented. When applied to proton-lead elastic scattering, the predicted real part of the Woods-Saxon potential is in substantial agreement with the experimentally determined phenomenological potential, although there are no adjustable parameters. In addition, the volume integrals of both real potentials are nearly identical