function y = synth_complex3(m) % m: mistuning in Hz % y=synth_complex1 - sum of mistuned harmonics of 100 Hz global sr; % sampling rate D = 0.1; % s - duration f0 = 100; % Hz - fundamental frequency p = [1, 2, 3, 4, 5, 6]; % partial ranks x = abscissa(D); y = x .* 0; for i = p y = y + sin(x * (i * f0 + m) * 2 * pi); end