samer@34: function y=harmsynth(N,A,F0,phi)
samer@34: % harmsynth - Additive synthesis with harmonic frequencies
samer@34: %
samer@34: % harmsynth ::
samer@34: %    N:natural    ~'size of buffers to produce',
samer@34: %    seq([[K]])   ~'sequence of amplitudes for each component',
samer@34: %    seq(nonneg)  ~'sequence of fundamental frequencies',
samer@34: %    [[K]]        ~'initial phases'
samer@34: % -> seq([[1,N]]) ~'sum of components'.
samer@34: %
samer@34: % Note: the number of components must remain constant
samer@34: 
samer@34: y=addsynth(N,A,F0*(1:size(A,1))',phi);