We optimize the first and second intrinsic hyperpolarizabilities for a 1D piecewise linear potential dressed with Dirac delta functions for N non-interacting electrons. The optimized values fall rapidly for N > 1, but approach constant values of βint = 0.40, γ + int = 0.16 and γ − int = −0.061 above N 8. These apparent bounds are achieved with only 2 parameters with more general potentials achieving no better value. In contrast to previous studies, analysis of the hessian matrices of βint and γint taken with respect to these parameters shows that the eigenvectors are well aligned with the basis vectors of the parameter space, indicating that the parametrization was well-chosen. The physical significance of the important parameters is also discussed.