print('freq = %f and thd = %f' % (i, ThdValue))
#Thd(samples, sampleRate)
#ThdN(samples, sampleRate)
plt.plot(np.abs(f))
#plt.legend(['test1', 'test2'], loc='best')
plt.grid()
plt.show()
elif c == 't':
#fs, dataI, dataQ = IOs.readWaveFile('D:/Documents/Samples/piano/UprightPianoKW-small-SFZ-20190703/samples/A5vH.wav')
#fs, dataI, dataQ = IOs.readWaveFile('D:/Documents/Samples/piano/Piano.A1.wav')
#print('sample rate : ', fs)
fs = 44100.0
freq = sd.ChromaticTone(48)
data = sd.KarplusStrong(int(2 * fs), fs, freq)
#data = np.sin(2*np.pi*freq/fs*np.arange(2*fs))
#print(f'freq of single tune = {freq}')
#plt.plot(data)
#plt.grid()
#plt.show()
#sdev.play(data, fs)
#status = sdev.wait()
#b = fd.RemezFilter(301, 250, 400, fs)
#data_flt = np.convolve(b, dataI, 'same')
#data_shift = data*np.sin(np.arange(data.size)*2*np.pi*420.0/fs)
