harmonic function

Construct harmonic functions for fitting harmonic trend model