% Skript zur Musterlösung
%
% Numerische Mathematik, SS 2013
%
% rekursive Berechnung von I_n = \int\limits_0^1 x^ne^x\,dx
%
% Ingolf Busch, Juni 2013

[ind,int] = my_integral(19);

fprintf(1,'Vorwaertsrekursion bis 19\n\n');
fprintf(1,'  i       Int     \n')
fprintf(1,'------------------\n')
for i=1:length(ind)
    fprintf(1,'  %2d     %10.7f\n',ind(i),int(i))
end
fprintf(1,'\n')

[ind,int] = my_integral(19,40);

fprintf(1,'Rueckwaertsrekursion von 40bis 19 \n')
fprintf(1,'  i       Int     \n')
fprintf(1,'------------------\n')
for i=1:length(ind)
    fprintf(1,' %2d     %10.7f\n',ind(i),int(i))
end