fsampling = 44100;        // sampling frequency
t = 0 : 1/fsampling : 4;  // 4 seconds at this frequency
N=size(t,'*');            //number of samples

f1=440;
sinwav1 = sin (2 * %pi * f1 * t);    				// a wave at 440 Hz
f2=460;
sinwav2 = sin (2 * %pi * f2 * t);    				// a wave at 460 Hz

waves = [sinwav1; sinwav2]; // Construct a 2*N matrix

savewave ('beep_L2.wav',waves,fsampling);

// Plot the power spectrum
s=waves(1,:)+waves(2,:); // We want a 1D signal (simpler for FFT)
y=abs(fft(s));
f=fsampling*(0:N/2)/N;   // Associated frequency vector
n=size(f,'*');
plot2d(f(1:n),y(1:n));