Fourier Transform using code

Question

Can someone help me? I need to obtain the Fourier transform of the following function:

$$ g_s(t) = \sum_{n=-\infty}^{+\infty} \text{rect}(50(t - 20 \pi n)) $$

It could be solve by code, like in Mathematica of WolframAlpha or Scilab.