function s=son_pur(f,v,d,sr)
	n = 1:(d*sr);
	s = v*sin(2.*pi.*f*n/sr);
endfunction

sr=8000;
f=440;
d=1;
v=0.9;

s=son_pur(f,v,d,sr);
[y,c]=stft(s,128);
c
s2=synthesis(y,c);
sound(s2);