clearvars;

x = 1;
y = exp(x);

tol = 1e-10;
deg = 0;
[h, member] = tay(x,deg);

error = [];
while member > tol
    deg = deg + 1;
    [h, member] = tay(x,deg);
    plot([x], [y], 'bo', [x], [h], 'ro');
    error = [error; y-h];
end

format long;
error
format short;

function [y, member] = tay(x,n)
    y = 1;
    member = 1;
    for i = 1:n
        member = ((x.^i)./factorial(i));
        y = y + member;
    end
end