|
MATHEMATICA
|
terms = 12;
E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}];
E4[x_] = 1 + 240*Sum[k^3*x^k/(1 - x^k), {k, 1, terms}];
E6[x_] = 1 - 504*Sum[k^5*x^k/(1 - x^k), {k, 1, terms}];
E8[x_] = 1 + 480*Sum[k^7*x^k/(1 - x^k), {k, 1, terms}];
a[n_] := (1/(24 n))*Sum[MoebiusMu[n/d]*SeriesCoefficient[E8[x]/E6[x] - E4[x]/E2[x], {x, 0, d}], {d, Divisors[n]}];
|