Exact volterra-series computation of nonlinear polarization in optical media