fs = 8192;  % Sampling rate (Hz)
timescale = 2*pi /fs; % Radians per sample
duration=2.0; %in seconds
t = timescale * (0:(duration*fs)); % scaled time for sinusoids

freqs=240*[1 2 3 4 5]; %overtone series, in Hz
initamp=[1 0 0 0 0];
initamp=initamp/sum(abs(initamp)); % limit amplitude to 1
Q=1000; % Quality of oscillators
damp=freqs/Q; % oscillators all with same Q
waves=exp(-damp'*t).*sin(freqs'*t); %each Fourier term, with its damping
stringwave=initamp*waves; %add them up with appropriate amplitudes
sound(stringwave,fs); % play