Please write a Python program that reads in a wave file as a numpy array, Fourier transforms it (with the correct frequency axis omega), applies a quadratic phase term exp(i A omega^2) in  the spectral domain (A can be chosen by the user), Fouier transforms back and writes the real part of the resulting numpy array to an output wave file.