x=2;
a=1;
max=3;
F=0;
R=0;
f=0;

for n=0:1:max
    if n == 0
        F=exp(a); 
        R=(exp(2)-F)/exp(2)
    else if n==1
            F=F+(((x-a).^1)/factorial(1))*(exp(a));
            R=(exp(2)-F)/exp(2)
        else if n==2
                F=F+(((x-a).^2)/factorial(2))*(exp(a));
                R=(exp(2)-F)/exp(2)
            else if n==3
                    F=F+(((x-a).^3)/factorial(3))*(exp(a));
                    R=(exp(2)-F)/exp(2)
                end
            end
        end
    end
    plot(F,n,'o')
        hold on
        plot(R,n,'x')
        hold on
end