function [an]=esround1(N)
an=zeros(N+1,1);
e=exp(1);
an(1)=(e-1)/e;
for n=0:N-1
    an(n+2)=1-(n+1)*an(n+1);
end
figure(1);
plot((0:N),an,'*');
xlabel('n'); ylabel('a_n')