function y = synthe(f0,t,Fe,f1,beta)
%        y = synthe(f0,t,Fe,f1,beta)
% synthetiseur FM de Chowning (profondeur = 1)
%---------------------------------------------

% quelques valeurs de f1 interesantes :
% f1 = [1, 1.5, 2, 2.5...]*f
% le cyber boomerang : f0 = 330; f1 = f0/60; beta = 220;
% le fender : f0 = 220; f1 = 1.01*f0; beta = 244;

h = 1/Fe;
th = 0:h:t;
tr = [0 .015 .3 t*[.3 .95 1]];
yr = [0  1  .5 .3 .1  0];
env = envelop(tr,yr,Fe);

yg = sin(2*pi*f0*th+beta*sin(2*pi*f1*th)/f1);
yg = yg.*env;
yg = yg/(max(abs(yg)));

f0  = f0+.8;
yd = sin(2*pi*f0*th+beta*sin(2*pi*f1*th)/f1);
yd = yd.*env;
yd = yd/(max(abs(yd)));

y = [yg;yd]';

