%
% Implementation des Fehlers in der L2-Norm der numerischen Loesung 
% des Problems -u'' + u = 1 in (0, 1) mittels stueckweise
% QUADRATISCHEN Finiten-Elementen
%
% (zur Serie 9 / Aufgabe 3)
%
% T. Mitkova, November 2009
%

function [L2_err] = L2_Error_Quad_FE_1D(x, u, uh)

% Formfunktionen
N1 = inline('...');
N2 = inline('...');
N3 = inline('...');

N = ... ; % Anzahl der Gitterpunkte
N_El = ... ; % Anzahl der Elemente

[quad_p, quad_w] = Quad_GaussLeg(); % Quadratur

L2_err = 0;

for el = ...
    
    h = ... ; % Elementlaenge
                         
    for i = 1 : size(quad_p, 2) % Schleife ueber den Quadraturpunkten
        
        % Auswertung der numerischen Loesung im Quadraturpunkt
        Uh = ... ;
         
         % Auswertung der exakten Loesung
         U = ... ;
        
        L2_err = L2_err + ... ;
    end
end;

L2_err = ... ;