N = 19;

fun = @(x, n) x.^n .* exp(-x);

z = zeros(N, 1);

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

z