I also used EMS,but I used the one from wiki,and it leads you to the answer in just 1 step. You have to calculate an integral and that's it.

]]>First we attempt to approximate the sum as accurately as we can:

with confidence in all the digits. Next a PSLQ was done to try to identify the number in terms of simple constants. This turned out to be fruitless. So the Euler Mclaurin formula was used next.

Basically the EMS is a formula that relates sums, integrals and derivatives together. It has many forms but the one we are interested in looks like this.

Often this form is of use for tough sums because it is often easier to integrate and differentiate. Plugging the above sum into that equation and asking mathematica to evaluate it produces a big mess. With some work you can get this out of it.

We can take the limit as n approaches infinity term by term and we are left with

so

and we are done.

]]>