w=50;          % częstotliwość próbkowania [Hz]
T=1;            % rozważany okres [s]
tt=1/w;         % odległość pomiędzy dwoma punktami [s]
t=0:tt:T-tt;    % punkty na osi OX [s]
n=length(t);    % liczba rozważanych punktów

x=sin(2*pi*t);  % rozważana funkcja
subplot(2,1,1);
plot(t,x,'+');

x1=fft(x);      
x1=abs(x1);     % moduł 

subplot(2,1,2);
stem(x1);

waitforbuttonpress;