function psound(c, s)
% Play the sound of frequencies s(j)*c, each for 1 second
Fs = 8192;
x=linspace(0,1,Fs);
y=[];
for j=1:length(s)
    freq=s(j)*c;
    y=[y sin(2*freq*pi*x)];
end
sound(y, Fs);