numerical fourier transform of an error function

I am trying to evaluate a fourier transform of an error
function f(x)=erf(x+16)-erf(x-16);
I used subs(maple('fourier',erf(x+16)-erf(x-16),x,w),w,
[0:0.001:1]');
All I need is get the value of w at [0:0.001:1] region.
Because the fourier transfrm of error function is not
explicit,it only gives me something like
0
1/2*exp(2/125*i)*fourier(erf(x),x,1/1000)-
1/2*exp(-2/125*i)*fourier(erf(x),x,1/1000)
1/2*exp(4/125*i)*fourier(erf(x),x,1/500)-
1/2*exp(-4/125*i)*fourier(erf(x),x,1/500)
1/2*exp(6/125*i)*fourier(erf(x),x,3/1000)-
1/2*exp(-6/125*i)*fourier(erf(x),x,3/1000)
1/2*exp(8/125*i)*fourier(erf(x),x,1/250)-
1/2*exp(-8/125*i)*fourier(erf(x),x,1/250)
1/2*exp(2/25*i)*fourier(erf(x),x,1/200)-
1/2*exp(-2/25*i)*fourier(erf(x),x,1/200)
1/2*exp(12/125*i)*fourier(erf(x),x,3/500)-
1/2*exp(-12/125*i)*fourier(erf(x),x,3/500)
1/2*exp(14/125*i)*fourier(erf(x),x,7/1000)-1/2*exp
(-14/125*i)*fourier(erf(x),x,7/1000)
............................................................
.......................................................
Could anyone tell me how to obtain the nemerical value of
this problem? Thank you all in advance.
.