In this article, we consider the optical second harmonic nonlinearity (corresponding to x(2)) associated with intersubband transitions in a doped parabolic multiquantum well structure, due to the anharmonic correction. The calculation indicates that the x(2)(2ω) is larger than that of the bulk GaAs by 1 to 2 orders of magnitude, and the resonant transition frequencies are in the middle infrared region