numeric_methods_1_2

% Вычисление самого интеграла
for n = 1:30
    syms x;
    f = x^n * exp(x - 1);
    I = int(f, 'x', 0, 1);
    fprintf('Шаг %d. Значение n: %d. Значение интеграла: %f.\n', n, n, I)
end

% Вычисление по частям
I = exp(-1);
   fprintf('Шаг 1. Значение n: 1. Значение интеграла: %f.\n', I)

for n = 2:30
   I_prev = I;
   I = 1 - n * I_prev;
   fprintf('Шаг %d. Значение n: %d. Значение интеграла: %f.\n', n, n, I)
end