I am trying to use the result that $p(omega)^n$ = Gamma(n, beta) where $p(omega)$ is the fourier transform of an $Exp(beta)$ distribution. My code does well for small n but as n increases the gamma starts to move across the domain faster. Here I have just set beta=1 for simplicity. I can not find the error in my code.