function [x_bar, y_bar, M] = matrice_inertie(E_x, E_y, G_norme)

    PI = sum(G_norme);
    
    x_bar = (G_norme.'*E_x)/PI;
    y_bar = (G_norme.'*E_y)/PI;
    
    M = [ (G_norme'*(E_x-x_bar).^2)/PI, (G_norme'*((E_x-x_bar).*(E_y-y_bar)))/PI
          (G_norme'*((E_x-x_bar).*(E_y-y_bar)))/PI, (G_norme'*(E_y-y_bar).^2)/PI ];

end