matthiasm@92: f = 50:800;
matthiasm@92: fs = 44100;
matthiasm@92: t = (0:(2*fs))/fs;
matthiasm@92: nFreq = length(f);
matthiasm@92: 
matthiasm@92: for iFreq = 1:nFreq
matthiasm@92:     f(iFreq)
matthiasm@92:     x = zeros(size(t));
matthiasm@92:     for iAdd = 1:100
matthiasm@92:         x = x + cos(2*pi*f(iFreq)*t*iAdd)*0.6^(iAdd-1);
matthiasm@92:         if f(iFreq)*iAdd > fs/2
matthiasm@92:             break
matthiasm@92:         end
matthiasm@92:     end
matthiasm@92:     x = x / max(abs(x));
matthiasm@92:     wavwrite(x, fs, sprintf('/Users/matthiasm/data/pyin/simplewavs/%iHz.wav', f(iFreq)));
matthiasm@92: end