function M = MassMatrix2D(p,t,h);

  n = size(p,2);    % Antal noder. (= antal kolonner i p)
  ntri = size(t,2); % Antal trianglar. (= antal kolonner i t)
  M = zeros(n,n);   % Initiera massmatris.

  for el = 1:ntri

    nodes = t(1:3,el);
    coords = p(:,nodes); 
    Me = ElementmMassMatrix2D(h);
    M(nodes,nodes) = M(nodes,nodes) + Me;

  end
  
% subroutines ----------------------------------------------------------------

function Me = ElementmMassMatrix2D(h)
%
% Komplettera med de rätta värdena på elementmassmatrisen, Me.
%
Me = [ 0.0 , 0.0 , 0.0;
       0.0 , 0.0 , 0.0;
       0.0 , 0.0 , 0.0 ];