function myEch = ech(N)  

sexe = binornd(1,0.3,1,N);

## on tire la taille des hommes
m0= 1.75
s0= 0.12;
tailleH = normrnd(m0, s0,1, N);

m1=1.68;
s1=0.08;
tailleF = normrnd(m1, s1,1, N);

%on constitue le vecteur de taille
taille = tailleH;
indexFemmes = find(sexe);

taille(indexFemmes) = tailleF(indexFemmes)

myEch = [sexe; taille];

return 