Intermodulation and harmonics are important in frequency analysis of nonlinear systems. In neuron research, most investigations are taken in studying synchronization between the external stimuli and the output of neuron, but harmonics and intermodulation are often ignored. In this paper, harmonics and intermodulation of the subthreshold FitzHugh-Nagumo neuron are investigated and their magnitudes are used to predict frequency response of the neuron. Furthermore, through analyzing the magnitudes of harmonics, the intrinsic frequencies of the neuron could be identified. A harmonic of a wave is a component frequency of the signal that is an integer multiple of the fundamental frequency, while intermodulation is the result of two or more signals of different frequencies being mixed together, forming additional signals at frequencies that are not, in general, at harmonic frequencies of either. Both harmonics and intermodulation are caused by nonlinear behaviors of systems. Recently, they are well investigated in biological applications such as improving human's hearing and vision system. However, in neuron models, harmonics and intermodulation are rarely concerned. Investigators always focus on the external synchronization between the input and output. In traditional study of the external synchronization, signals at continuous frequencies are introduced to stimulate neuron to get corresponding output. By calculating the spectrum or Fourier coefficient of the output at each input frequency, the frequency response is obtained. So harmonics and intermodulation of the applied input signal have been ignored. In this paper, a sinusoid input is applied to FitzHughNagumo (FHN) neuron to make it remain in subthreshold oscillations and the corresponding harmonics and intermodulation are studied. They are incorporated to predict frequency response of the neuron. Also by identifying the maximum magnitude of harmonics, we can get the special input frequencies where the neuron could be more easily excited. Moreover, it demonstrates that the output response of a FHN neuron may not be maximum at the input frequency. Some output frequencies generated from harmonics or intermodulation may have significant effects on the output response which cannot be ignored in frequency analysis of the neuron.