function fft1(n) n dt=0.01 fr=0.5 w=2*pi*fr; t=dt*(0:n-1); tl=dt*n df=1/tl f=df*(0:n-1); d=sin(w*t); afd=abs(fft(d)); subplot(2,1,1); plot(t,d); grid on subplot(2,1,2); plot(f,afd); axis([0 4.5]); grid on endfunction