clear all
%exp(1)
tic
n=1000;
format long
ee=exp(1);

ee1=(1+1/n)^n;
[ee,ee1]
norm(ee-ee1)/norm(ee)*100


%1+1+1/2!+1/3!+...
fac=1;
S(1)=2;

Niter=100;
E(1)=(1+1)^1;

for ia=1:Niter
    fac=fac/(ia+1);
    S(ia+1)=S(ia)+fac;
    E(ia+1)=(1+1/(ia+1))^(ia+1);
   
end;

x=[1:1:Niter+1]
figure(10)
plot(x,S,'r',x,E,'b',x,ee*exp(0*x),'--r')



toc