Greg Fee, Catalan's Constant [Ramanujan's Formula…
“ The implementation of the algorithm is: inline void const_catalan (bigfloat & y) { bigfloat p; bigfloat t; int i = 1, j = 3; // j = 2*i+1 // y = t = p = 1/2 divide (y, 1, 2) ; t.assign (y) ; p.assign (y) ; // while t is greater than the desired accuracy while (!t.is_approx_zero () ) { // do // p = p * (i/j) ; // t = (t * i + p) / j; // y = y + t; // i++; j+=2; multiply (p, p, i) ; divide (p, p, j) ; multiply (t, t, i) ; add (t, t, p) ; divide (t, t, j) ; add (y, y, t) ”
