> restart;
> with(inttrans):
> assume(lambda>0,a>0,k>0):
> z:=exp(-a*x^2-I*b*x);
>
> F1:=factor((fourier(z,x,k)));