its very good - explains finally why we need fftshift twice. However, it should be stressed that the method works only for data of even size. For odd size data, the amplitude is correct but spectral phase will be wrong.
Excellent description of the theory of the FFT.
But my code threw up errors...cannot define functions at this point?? Could not run it now. I will have to look into it a little.