Differentiation in N Dimensions with the Fast Fourier Transform (FFT)

This is a quick follow-up to my recent article on taking derivatives using the FFT Fourier transform. As promised, I give the generalization from 1 to 𝑁 dimensions. The approach is quite the same, just the numpy functions are different and we have to use meshed grids.