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