function s = chowning(f0,T,Fe)
%        s = chowning(f0,T,Fe)
% synthetiseur FM de Chowning, stereo

%------------------------------------------------------
%Preciser ici les valeurs de f1 et beta
% quelques valeurs de f1 interesantes :
% f1 = f0, 1.5*f0, 2*f0, 2.5*f0, ...
% le cyber boomerang : f0 = 330; f1 = f0/60; beta = 220;
% le fender (:~) : f0 = 220; f1 = 1.01*f0; beta = 244;

f1 = 1.5*f0; beta = f0;
%------------------------------------------------------

dt = 1/Fe;
t = 0:dt:T;
n = length(t);

te = T*[0 .015 .3 .95 1];
a = [0  1  .5 .1  0];
env = envelop(te,a,Fe);

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

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

s = [sg;sd]';

