sum ln(k)/n*((1-FractionalPart[n/k])(1-k/n*FractionalPart[n/k])-1/2) from k=1 to n