function hsound(c, s)
% Play the sound of frequencies s(1)*c, s(2)*c, ..., all together for 1 second
Fs = 8192;
x=linspace(0,1,Fs);
y=sin(2*s(1)*c*pi*x);
for j=2:length(s)
    freq=s(j)*c;
    y=y+s(j)*sin(2*freq*pi*x);
end
y=y/max(abs(y));
sound(y, Fs);