% 定义函数句柄
f = @(x, n) x.^n .* exp(x);
n_values = 0:20; % 计算 I_0 到 I_20
I = zeros(size(n_values)); % 存储积分结果
for n = n_values
integrand = @(x) f(x, n);
I(n+1) = exp(-1) * integral(integrand, 0, 1);
end
% 打印结果
disp('计算结果:');
for n = n_values
fprintf('I_%d = %.8f\n', n, I(n+1));
end
% 绘制结果图形
figure;
plot(n_values, I, '-o');
title('I_n 计算结果');
xlabel('n');
ylabel('I_n');
grid on;
09-14
2933
08-09
5万+