f[ϵ_]:=Integrate[((E^(-x))/(1+ϵ x)),{x,0,∞}]

N[f[1/20],40]

Table[N[Sum[(((-1))^n)n!(((1/40))^n),{n,0,d}],40],{d,0,120}]//Column