A new method for the extraction of coefficients in an nth-order polynomial function is presented. The method may be used as an efficient tool to calculate any kth-order harmonic distortion in semiconductor devices for any input amplitude, without having to perform Fourier or AC analysis. To illustrate its utility, the calculation procedure has been applied to determine the harmonic distortion produced by the I-V output characteristics of an Independently Driven Double-Gate MOSFET. The accuracy of the results is analyzed as compared to the Weighted Least Squares algorithm. The calculated distortion harmonics are compared to those obtained by the use of traditional Fourier analysis.