Consider the Fourier transform $F(\omega)$ of the function $f(t)$. The magnitude of $F(\omega)$ depends on $\omega$ and thus also depends on the scale of the $t$-axis. For example, when $f_1(t)$ is a box function which is 1 between $[0,1]$ s and $f_2(t)$ is a box function which is 1 between $[0,1]$ ns = $[0,10^{-9}]$ s, then the amplitude of $F_1(\omega)$ is $10^9$ times larger than the amplitude of $F_2(\omega)$. This is because the time scale is reduced by a factor $10^9$ in the second case.
However, no difference would be obtained when you do a numerical Fourier transform. If you have 100 points of the function $f_1(t)$ equally spaced between $[0,10]$ and also 100 points of the function $f_2(t)$ equally spaced between $[0,10\times10^{-9}]$, then numerically the amplitude of the Fourier transform of both functions is the same (which is not correct). Numerically, you only consider the Fourier transform of the "y-values" $f(t_i)$, so the information about the time scale is lost...
You can solve this problem by multiplying the result with $10^{-9}$ in the second case because the Fourier transform of the rectangular function is proportional with $1/\omega$. However, in general, this proportionality is not present... How to solve this problem for a general function $F(\omega)$? For example, when the Fourier transform is proportional to $e^\omega$ or proportional to $1/\omega^2$ or ...? Is there a numerical approach to take the scale of the time axis into account?