QuickFFTSurface(m, anObject) displays the modulus of the Fourier transform of a 2-d array of real numbers as the z coordinate of a surface z(x,y). The phase of the Fourier transform is mapped into the color palette of the surface. The 2-d array m has to be a matrix.

set f to "sin(64*x*y)+sin(x*y*x*y)"
set {x, y, z} to QuickFormulaMatrix({0, pi / 2}, {0, pi / 2}, f)
QuickFFTSurface(z, 0)

Import script