function v = LoadVector(p)

%------------------------------------------------------------------------------
%                                                                            
%                                                                    
% Syfte:    Assemblerar högerledet, v, till ekvationen Mc = v,
%           givet en partition av [a,b] i N stycken delintervall, 
%           och en funktionsfil 'f.m' som definierar den funktion, 
%           f(x), som vi vill beräkna L2-projektionen av.         
%           
% Indata:   p:  Nodkoordinater.       (Vektor med dimension 1 x N+1)   
%                                                                      
%              ex.: p = 0:0.2:1 ger en likformig partition av [0,1]           
%              
% Anropade funktioner: f.m (du måste skapa f.m själv)                  
%                                                                      
%            ex. på hur f.m kan se ut:  function y = f(x)              
%                                       y = sin(x);                    
%                                                                      
%------------------------------------------------------------------------------

N = length(p) - 1;      % Antal delintervall.
v = zeros(N+1, 1);      % Initiera högerledsvektorn.

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

  % Trapetsmetoden:

  v(i)   = v(i)   + f(p(i))*h/2; 
  v(i+1) = v(i+1) + f(p(i+1))*h/2; 
  
end