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

y=addsynth(N,A,F0*(1:size(A,1))',phi);
