Riemann zeta function on the critical line

scale = 1000000;
 xres = .00001;
 x = Exp[Range[0, Log[scale], xres]];
 RealPart = Log[x]*FourierDST[-(SawtoothWave[x] – 1)*x^(-1/2)];
 ImaginaryPart = Log[x]*FourierDCT[-(SawtoothWave[x] + 1)*x^(-1/2)];
 datapointsdisplayed = 300;
 ymin = -15;
 ymax = 15;
 g1 = ListLinePlot[
 Sqrt[scale]*{RealPart[[1 ;; datapointsdisplayed]],
 ImaginaryPart[[1 ;; datapointsdisplayed]]},
 PlotRange -> {ymin, ymax}, DataRange -> {0, 68.00226987379779},
 Filling -> Axis];
 Show[Flatten[{g1,
 Table[Graphics[{PointSize[0.013],
 Point[{N[Im[ZetaZero[n]]], 0}]}], {n, 1, 16}]}],
 ImageSize -> Large]