# 
restart;
S:=n->Sum(1/k!,k=0..n);
S(2);
value(S(2));
value(S(4));
S10:=value(S(10));
exp(1.0);
evalf(S10);
epsilon:=0.00000001;
2/epsilon;
for n from 7 to 17
do
    print(n,n!,n^n*exp(1.0-n));
od:
S12:=value(S(12));
epsilon;
Digits:=50;
evalf(exp(1)-S12);
T:=(x,n)->Sum(x^k/k!,k=0..n);
A1:=value(T(5,12));
exp(5.0);
evalf(A1);

#