function y = synth_complex2 % y=synth_complex1 - sum of sines (odd partials of 100 Hz) global sr; % sampling rate D = 0.1; % s - duration f0 = 100; % Hz - fundamental frequency p = [5, 7, 9, 11, 13]; % partial ranks x = abscissa(D); y = x .* 0; for i = p y = y + sin(x * i * f0 * 2 * pi); end