cosx = (exp(ix) + exp(-ix))/2
sinx = (exp(ix) - exp(-ix))/2i