Please write a short Python script that reads in the first 5 seconds of a wav file as a float numpy array, Fourier transforms the data, multiplies a spectral phase term exp(i A f^2) with f the frequency axis obained from the sample rate with the fftfreq function, inverse Fourier transforms and writes the real part of the resulting numpy array to an output wav file.