% % % % % % % % % % % % % % % % % % % % %
%                                       %
%  Programmieraufgabe A 22, SS2007      %
%  Höhere Numerische Mathematik         %
%                                       %
% % % % % % % % % % % % % % % % % % % % % 


% benoetigt folgende Funktionen: 
% romberg.m
% f1.m
% f2.m
% f3.m
% f4.m

clear all
format long

relerr = 10^-5;



a = 0; 
b = 1; 
[T,T1,fehler] = romberg(@f1,a,b,relerr);
disp('Loesung a:')
T1 = T1
[i,k]=size(T);
disp(['Es wurden i = ',num2str(i),' Zeilen und k = ',num2str(k),' Spalten '])
disp('des Romberg Schemas berechnet.')
disp(['Dabei wurde eine Genauigkeit von ',num2str(fehler),' erreicht.'])
disp(' ')

a = 0; 
b = 2/pi;
[T,T1,fehler] = romberg(@f2,a,b,relerr);
disp('Loesung b:')
T1 = T1
[i,k]=size(T);
disp(['Es wurden i = ',num2str(i),' Zeilen und k = ',num2str(k),' Spalten '])
disp('des Romberg Schemas berechnet.')
disp(['Dabei wurde eine Genauigkeit von ',num2str(fehler),' erreicht.'])
disp(' ')


a = 1; 
b = 2; 
[T,T1,fehler] = romberg(@f3,a,b,relerr);
disp('Loesung c:')
T1 = T1
[i,k]=size(T);
disp(['Es wurden i = ',num2str(i),' Zeilen und k = ',num2str(k),' Spalten '])
disp('des Romberg Schemas berechnet.')
disp(['Dabei wurde eine Genauigkeit von ',num2str(fehler),' erreicht.'])
disp(' ')


a = 0; 
b = 1; 
[T,T1,fehler] = romberg(@f4,a,b,relerr);
disp('Loesung d:')
T1 = T1
[i,k]=size(T);
disp(['Es wurden i = ',num2str(i),' Zeilen und k = ',num2str(k),' Spalten '])
disp('des Romberg Schemas berechnet.')
disp(['Dabei wurde eine Genauigkeit von ',num2str(fehler),' erreicht.'])
disp(' ')

