harmonic analysis funtion