function y = clavecin(f,T,Fe,inharm)
%        y = clavecin(f,T,Fe,inharm)
% imitation clavecin avec inharmonicite
%--------------------------------------

h = 1/Fe;
th = 0:h:T;

n = 42; N = 1:n;
% poids des partiels
am = N.^(-.8);

% calcul inharmonicite
cent = 2^(1/1200);
harm = N.*cent.^(inharm*N.^2);
f = harm*f;

% somme des partiels
y = zeros(1,length(th));
deb = .01;
for k = 1:n
    z = sin(2*pi*f(k)*th);
    % avec enveloppe variable par harmonique
    tr = [0 deb deb*[1.001 1.002 1.003 1.004]+116*abs(sin(pi*k/10))^(1/3)/(k+8)^1.2*[.3 .6 .8 1] T];
    yr = [0 1 .4 .1 .02 .01 0];
    env = envelop(tr,yr,Fe);
    z = z.*env;
    y = y + am(k)*z;
end

y = .98*y/(max(abs(y)));
