%
% Main-Programm zur Aufgabe 4, Serie 2
%
% Vorlesung "Einfuerung in die Numerik", Uni Basel, FS 09
%
% T. Mitkova, Februar 2009
%

clear all
clf

xs = 0 : 0.2 : 0.8;

for i = 1 : length(xs)-1
    x = xs(i) : 0.01 : xs(i+1);
    
    l = LinInterpol(xs(i), xs(i+1), x);
    q = QuadInterpol(xs(i), xs(i+1), x);
    
    subplot(1,2,1), hold on
    plot(x, f(x), 'k-', x, l, 'k-.'), axis([-0.03 0.83 -0.5 1.5]), title('Stueckweise lin. Interpol.')
    
	subplot(1,2,2), hold on
    plot(x, f(x), 'k-', x, q, 'k-.'), axis([-0.03 0.83 -0.5 1.5]), title('Stueckweise quad. Interpol.')
end