syms t f;                              
syms a positive
s = exp(-a*t)*heaviside(t);
S = fourierf(s)                 % Find S(f)