The method of perturbation used for calculating the electron distribution function is discribed. The small signal gain and saturation gain of harmonic free electron laser operating in Compton region with planar undulator are obtained. Analyses and discussions on the model and results are given in this paper.