function s = flute(f1,T,Fe)
%        s = flute(f1,T,Fe)
% son flute
%----------------------------------------------------

h = 1/Fe; th = 0:h:T;
nt = length(th);

a = [1000 50 80 10 5 2 .1 1];
nh = length(a); N = 1:nh;
f = N*f1;

s = synthad(a,f,0*f,T,Fe);

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

% ajout de souffle
br = bruit(T,Fe);
fn = f1/Fe*2; wn = (fn*[.8 1]).^(2);
b = fir1(50,wn);
br = filter(b,1,br);
t  = T*[0 .05 .8 1];
a  =   [0  1 .5 0];
env = envelop(t,a,Fe);
br = br.*env;

s = s+br/7;
s = s/(max(abs(s)));
