(*Mathematica start*)nn = 200;
f[t_] = D[RiemannSiegelTheta[t], t];
Plot[Sign[RiemannSiegelZ[t]]*
Abs[Zeta[1/2 + I*t]*
Total[Table[
Total[MoebiusMu[Divisors[n]]/Divisors[n]^(1/2 + I*t - 1)]/
n, {n, 1, nn}]]]/(f[t] + HarmonicNumber[nn]), {t, 0, 60},
PlotStyle -> Thickness[0.004], ImageSize -> Large,
PlotRange -> {-4, 4}]
(*end*)