clear

format compact

N = 19;
fun = @(x, i) x.^i.*exp(-x);

z = zeros(N,1);

tic();
for n = 1:N
    ff = @(x) fun(x, n);
    z(n) = integral(ff, 0, 1);
end
toc()