Mathematica Stack Exchange is a question and answer site for users of Mathematica. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

An earlier post was migrated to Physics.SE and we determined that if this converged, we would have a constant close to the mass of the proton in kilograms. This constant is constructed by summing the reciprocals of the products of all the numbers between consecutive squares.

Per this:

SumConvergence[10^-26/Pochhammer[m^2 + 1, 2 m], m]  

(* True *)

we have convergence and thus a constant. However, the following:

Sum[10^-26/Pochhammer[m^2 + 1, 2 m], {m, \[Infinity]}]  

outruns the precision.

My question: Is there a way to use this sum to create the constant to a specified accuracy?

share|improve this question
To be sure, the common opinion on of the meaning of this sum is that it is a sheer coincidence that it is close to the proton mass in kg. – Sjoerd C. de Vries Nov 13 '12 at 11:17
@SjoerdC.deVries, I just got the Popular Question badge --- 1000 views and only 2 votes. Yes, they don't like it. – Fred Kline Nov 13 '12 at 17:32
up vote 2 down vote accepted

Clearly, you can factor out the tiny portion of your constant, so:

NSum[1/Pochhammer[m^2 + 1, 2 m], {m, 1, ∞}, WorkingPrecision -> 50] 1*^-26


...products of all the numbers between consecutive squares.

Product[j, {j, (k - 1)^2, k^2}]
   Pochhammer[(-1 + k)^2, 2 k]

This doesn't resemble what you have.

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.