function M = MassMatrix(p)

%------------------------------------------------------------------------------
%                                                                      
% Syntax:   M = MinMassmatrisAssemblerare(p)                           
%                                                                      
% Syfte:    Assemblerar massmatrisen, M, givet en partition av [a,b]   
%           i N stycken delintervall.                                  
%                                                                      
% Indata:   p:  Nodkoordinater.       (Vektor med dimension 1 x N+1)   
%                                                                      
%              ex.: p = 0:0.2:1 ger en likformig partition av [0,1]    
%                               med fem delintervall                   
%                                                                      
% -----------------------------------------------------------------------------                                                                     
N = length(p) - 1;      % Antal delintervall.
M = zeros(N+1, N+1);    % Initiera massmatris.

for i = 1:N 		% Loopa över delintervallen.
 
  h = p(i+1) - p(i);         % Delintervallets längd.

  % Trapetsregeln:

  M(i, i)     = M(i, i)     + 1.0*h/2;
  M(i, i+1)   = M(i, i+1)   + 0.0*h/2;
  M(i+1, i)   = M(i+1, i)   + 0.0*h/2;
  M(i+1, i+1) = M(i+1, i+1) + 1.0*h/2;
  
end

