§3.5 Summation Formulas: The Real Case

An exponential summation formula (S. 77).

ϕ[
          u_] := (E^((2/(((1 - u))^2)) - (1/(2 (u^2)))));
n = 60; 
    k = N[Range[n - 1, 1, -1]];
c = 
        Append[k + ϕ[(k/n)] + 1, 1];
w=((Append[1+(n^(-1))ϕ'[(k/n)],1])^(1/2));
A = 
        Table[((w[[j]] w[[k]])/((c[[j]] + c[[k]] - 1) (c[[j]] + c[[k]])/2 - (c[[k]] - 1))), {j, n}, {k, n}];
Norm[A]