% Simpson, Aufgabe 15
%
%Thomas Dettbarn	5610641
%Marina Chanuchowa	5622011
%Bibiana Gluski		5600816
%
%% Eingabe:	a: linke Intervallgrenze
%%		b: rechte Intervallgrenze
%%		eps: Epsilon, die Genauigkeit
%% Ausgabe:	Das Simpsonintegral auf diesem Intervall mit hinreichender
%%		Genauigkeit
function Ads=Simpson(a,b,eps);
global n	%% verwende n als globale Variable
L=(b-a);	%% Berechne die Laenge des Ausgangsintervalles
delta=(15*eps)/L;%% Bereche delta
Sab=((b-a)/6)*(f(a)+4*f(a+((b-a)/2))+f(b));
		%% Berechne das Simpsonintervall auf dem Ausgangsintervall
n=0;		%% Starte mit n=0
erg=fopen('AdSimpson-Anwendung.txt','w+');
fprintf(erg,'Ergebnisse der adaptiven Simpsonregel\n');
fprintf(erg,'Rekursions                 Intervalle               Intervall-\nAufrufe Tiefe           a             b               Laenge\n');
fclose(erg);
		%% Schreibe die Ueberschriften in die Ausgabedatei
Ads=AdSimpson(1,a,b,Sab,delta); %% Starte das Adaptive Simpsonverfahren