function s = synthadl(a,f,p,T,Fe)
%        s = synthadl(a,f,p,T,Fe)
% synthese additive
% cette fonction cree un son de duree T, 
% compose des partiels f(n), d'amplitude a(n)
% et de phase a l'origine p(n).
% Fe est la frequence d'echantillonnage
% effet Leslie rajoute, son stereo
%---------------------------------------------

mu = .01;
eta = .3;
fL = 4;

% creation du vecteur temps discret
dt = 1/Fe;
t = 0:dt:T;
n = length(t);

% canal gauche
sg = zeros(1,n);
K = length(f);
for k = 1:K
    z = (1 - eta*sin(2*pi*fL*t)).*sin(f(k)*(2*pi*t - mu*sin(2*pi*fL*t)/fL));
    sg = sg+a(k)*z;
end
sg = .99*sg/max(abs(sg));

% canal droit
sd = zeros(1,n);
for k = 1:K
    z = (1 + eta*sin(2*pi*fL*t)).*sin(f(k)*(2*pi*t + mu*sin(2*pi*fL*t)/fL));
    sd = sd+a(k)*z;
end
sd = .99*sd/max(abs(sd));

% enveloppe
t = T*[0 .09 .98 1];
a = [0 1 .9 0];
env = envelop(t,a,Fe);

sg = sg.*env;
sg = sg/(max(abs(sg)));
sd = sd.*env;
sd = sd/(max(abs(sd)));

s = [sg;sd]';

