import sounddevice as sd                # conda install sounddevice 
                                        # conda install PortAudio
import numpy as np                      # conda install numpy
import matplotlib.pyplot as plt         # conda install matplotlib

pi=np.pi
# (time, amplitude) array
fs=10000                                        #sample rate
t_max=1                                         # sound duration
time_array = np.arange(0, t_max - 1/fs, 1/fs)

freq=300
sound1=np.sin(2*pi*freq*time_array)                 # tonique
sound2=np.sin(2*pi*1.5*freq*time_array)             # quinte
sound=sound1+sound2
sd.play(sound, fs)
sd.wait()  # Wait until the sound finishes playing
plt.plot(time_array,sound)
plt.xlim(0, 0.01)
plt.show()