for n=2:12 N=2^n a=-1im*2*pi/N*(0:N-1) b=0:N-1 A=exp.(a*b') Ainv=conj(A)/N x=rand(N) println("N=",N) @time xhat=A*x # println("x=",x) # println("xhat=",xhat) end