function q=my_quad2Dgen( f, p, t)

nnodes = size(p, 2);  % Antal noder. (= antal kolonner i p)
ntri = size(t, 2);    % Antal trianglar. (= antal kolonner i t)

q = 0;                % Initiera q = integralens värde


for n = 1:ntri  % Loopa över alla trianglar i nätet
  
  % Ta reda på hörnens nodnummer:
  
  i =               % Nodnummer för hörn nummer 1 i triangel nummer n.
  j =               % Nodnummer för hörn nummer 2 i triangel nummer n.
  k =               % Nodnummer för hörn nummer 3 i triangel nummer n.
                
  % Ta reda på hörnens koordinater:
  
  xnodi =           % x-koordinat för nod nummer i.
  xnodj =           % x-koordinat för nod nummer j.
  xnodk =           % x-koordinat för nod nummer k.
  
  ynodi =           % y-koordinat för nod nummer i.
  ynodj =           % y-koordinat för nod nummer j.
  ynodk =           % y-koordinat för nod nummer k.
  
  % Räkna ut triangelns mittpunkt:  
  
  xmitt = 
  ymitt = 
  
  % Räkna ut triangelns area = 
  % halva arean av den parallellogram
  % som spänns av vektorerna mellan
  % nod j och nod i samt mellan nod k
  % och nod i
  
  triangelarea = ((xnodj - xnodi)*(ynodk - ynodi) - ...
                  (xnodk - xnodi)*(ynodj - ynodi))/2;
  
  % (... gör att man kan fortsätta skriva på nästa rad)
  
  % Räkna ut integralens bidrag från triangel
  % n och lägg detta till q:

  q = q + 
 
end
