% geometrie du modèle
Coord (1,:) = [0  0  0];
n_fix = 1;

% element flex
Coord (2,:) = [ Lf 0 0 ];
Elem(1,:) = [ 1 2 ];
eFlex = 1;

% elements poutre
n = size(Coord,1);  ne = size(Elem,1);
for i = 1:np1
    Coord (n+i,:) = [ Lf+Lp1*i/np1     0    0 ];
    Elem(ne+i,:) = [ n+i-1  n+i ];
end
n_act = n+np1;

n = size(Coord,1);  ne = size(Elem,1);
for i = 1:np2
    Coord (n+i,:) = [ Lf+Lp1+Lp2*i/np2     0    0 ];
    Elem(ne+i,:) = [ n+i-1  n+i ];
end
n_mass = n+np2;

n = size(Coord,1);  ne = size(Elem,1);
ePoutre = [ 2:ne ];