clc
clear all
close all

x = linspace(-3, 3, 1000);
n = input("Insert the degree: ");
figure
hold on
y = exp(x);
plot(x, y, 'm')
grid on
pause
y = 0;
for k = 0:n
    y = y + x.^k/factorial(k);
    plot(x, y, 'c')
    fprintf("k = %d\n", k)
    pause
end