Sounds

// At first we create 0.5 seconds of sound parameters.
t=soundsec(0.5);
// Then we generate the sound.
s=sin(440*t)+sin(220*t)/2+sin(880*t)/2;
[nr,nc]=size(t);
s(nc/2:nc)=sin(330*t(nc/2:nc));
// We can easily make a Fourier analysis of it.
;analyze(s);

// Save the file in WAV format.
// we renormalize s in order to check that save+load is invariant
s=s-sum(s)/prod(size(s)); s=s/max(abs(s));
savewave("test.wav",s);
// Load it back.
s1=loadwave("test.wav");
if maxi(abs(s1-s)) < 1.e-4;end
// Now we can make a complete picture of the sound.
;mapsound(s);