N=4 a=-1im*2*pi/N*(0:N-1) b=0:N-1 A=exp.(a*b') Ainv=conj(A)/N x=[1,2,3,4] println("N=",N) @time xhat=A*x println("x=",x) println("xhat=",xhat)