Steven Foster Clark primefourierseries without zeta zeros

(* Mathematica start, with modification received from Steven Foster \
Clark on 22.10.2016 URL: https://www.primefourierseries.com/?p=1# \
comment-56 *)
f = 5;
nn = 30;
delta = 1/20;
MatrixForm[
  A = Chop[N[
     Total[Transpose[
       Table[Table[
         If[Mod[n, k] == 0,
          Table[Sum[Exp[-a*b/(n/k)*2*Pi*I], {b, 1, f n/k}], {a, 1, nn,
              delta}]*MoebiusMu[n/k], 0]
         , {k, 1, nn}], {n, 1, nn}]]]]]];
Show[ListLinePlot[Re[Total[A/Range[nn]]],
  PlotRange -> {-f*2, f*Log[nn]}],
 ListLinePlot[Table[f*Log[n], {n, 1, nn, delta}]]]
(* end *)

(*Slightly improved version 23.10.2017*)
(* Mathematica start, with modification received from Steven Foster \
Clark on 22.10.2016 URL:https://www.primefourierseries.com/?p=1# \
comment-56 *)
f = 2;
nn = 20;
delta = 1/20;
MatrixForm[
  A = Chop[N[
     Sum[Sum[If[Mod[n, k] == 0,
         Table[Sum[Exp[-a*b/(n/k)*2*Pi*I], {b, 1, f*n/k}], {a, 1, nn,
            delta}]*MoebiusMu[n/k], 0]/n, {n, 1, nn}], {k, 1, nn}]]]];
Show[ListLinePlot[Re[A], PlotRange -> {-f*2, f*Log[nn]}],
 ListLinePlot[Table[f*Log[n], {n, 1, nn, delta}]]]
(*end*)

(* Better Mathematica start,with modification received from Steven Foster \
Clark on 23.10.2016 URL:https://www.primefourierseries.com/?p=1# \
comment-56 *)
f = 5;
nn = 20;
delta = 1/20;
MatrixForm[
  A = Re[Chop[
     N[Table[Sum[
        Sum[If[Mod[n, k] == 0,
           Sum[Exp[-a*b/(n/k)*2*Pi*I], {b, 1, f*n/k}]*MoebiusMu[n/k],
           0]/n, {n, 1, nn*2}], {k, 1, nn*2}], {a, 1, nn, delta}]]]]];
Show[ListLinePlot[A, PlotRange -> {-f*2, f*Log[nn]}],
 ListLinePlot[Table[f*Log[n], {n, 1, nn, delta}]]]
(*end*)