It still takes a little ingenuity to compute thousands of digits, like in Crandall’s code, tweaked by me, found at
Try to beat these MRB constant records!
Timing[etaMM[m_, pr_] :=
Module[{a, d, s, k, b, c},
a[j_] := N[(-PolyLog[1, -j]/(j + 1))^m, pr];
n = Floor[1.32 pr];
d = Cos[n ArcCos[3]];
{b, c, s} = {-1, -d, 0};
Do[c = b - c;
s = s + c a[k];
b = N[(k + n) (k - n) b/((k + 1) (k + 1/2)), pr], {k, 0, n - 1}];
Return[N[s/d, pr] (-1)^m]];
eta[s_] := (1 - 2^(1 - s)) Zeta[s];
eta1 = Limit[D[eta[s], s], s -> 1];
MRBtrue = mm;
prec = 1500;(*number of digits. 1500 will take around 200 seconds.*)
MRBtest =
eta1 - Sum[(-1)^m etaMM[m, prec]/Gamma[m + 1], {m, 2,
Floor[.45 prec]}]; Print[MRBtest]]