%
% Test der Implementierung der Gauss-Quadratur der Ordnung 4
% (zur Serie 6, Aufgabe 3)
%
% Vorlesung "Einfuerung in die Numerik", Uni Basel, FS 09
%
% T. Mitkova, Maerz 2009
%

clear all

f_test1 = inline('exp(x)', 'x');
f_test2 = inline('1./sqrt(x)', 'x');

I_exact = ..... ; % exakter Wert des Integrals I1

for i = 1 : 10
    h(i)=2^(-i); % Schrittweite
    
    % Numerische Approximation von I1
    QuadGauss_Test1(i) = Quad_SumGauss4(..., ..., ..., ...);
        
    % Numerische Approximation von I2 
    QuadGauss_Test2(i) = Quad_SumGauss4(..., ..., ..., ...);
    
    if i > 1
        % Fehler der Berechnung von I2
        error_test2(i-1) = abs(....);
    end
end

figure
loglog(h, ... , 'ko-', h, h.^4)
legend('Gauss-Quadratur', 'Ordnung 4', 'Location','SouthEast')
title('Test 1: f(x) = exp(x)')

figure
loglog(h(2:10), ... , 'ko-', h, h.^(1/3))
legend('Gauss-Quadratur', 'Ordnung 1/3', 'Location','SouthEast')
title('Test 2: f(x) = 1/sqrt(x)')