The cubic response equations governing nonequilibrium and equilibrium solvation in the continuum approach are derived and implemented. We consider solvent effects on the frequency dependent second hyperpolarizability. We illustrate the effects of optical and static dielectric constants through sample calculations on H2O and H2CO.