% file rekursion.m
% Berechne Integrale der Form p_n = \int_0^1 dx x^n \exp(x)
% mit Hilfe der Rekursionsformel p_{n+1} = e - (n+1)*p_n
e=exp(1);
p(1)=1; % Startwert
for n=1:19
  p(n+1)=e-(n+1)*p(n);
end
[[1:20]' p']