restart;
f:=n->(1+1/n)^n;
f(1);
f(2);
evalf(f(3));
for k from 1 to 10
do
    print("k=",k,"f(k)=",f(k),evalf(f(k)));
od;
g:=n->(1+1/n)^(n+1);
for k from 1 to 10
do
    print("k=",k,"g(k)=",g(k),evalf(g(k)));
od;