%
% Berechnung des Integrals einer Funktion f von a bis b naeherungsweise
% mit der summierten Gauss-Quadratur vierter Ordnung
%
% Vorlesung "Einfuerung in die Numerik", Uni Basel, FS 11
%

function Int = Quad_SumGauss4(f, a, b, h)

% Quadraturpunkte
xi_1 = ... ;
xi_2 = ... ;

% Gewichte
omega_1 = ... ;
omega_2 = ... ;

N = ... ;

x1 = ... ; % Argumente: a + h*xi_1 + j*h, j = 0, ..., N-1
x2 = ... ; % Argumente: a + h*xi_2 + j*h, j = 0, ..., N-1

f1 = f(x1);
f2 = f(x2);

% Summ. Gauss-Quadratur vierter Ordnung
Int = h*(omega_1* ... + omega_2* ...);

end
