Advertisement
MatsGranvik

von Mangoldt function Fourier transform Dirichlet series

Jun 19th, 2013
620
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.32 KB | None | 0 0
  1.  
  2. scale = 10000;
  3. Print["Counting to 60"]
  4. Monitor[g1 =
  5. ListLinePlot[
  6. Table[Re[
  7. Zeta[1/2 + I*t]*
  8. Total[Table[
  9. 1/N[Log[scale]]*1/n*
  10. Total[1/Divisors[n]^(1/2 + I*t - 1)*
  11. MoebiusMu[Divisors[n]]], {n, 1, scale}]]], {t, 0, 60,
  12. N[1/6]}], DataRange -> {0, 60}, PlotRange -> {-0.02, 1.6}],
  13. Floor[t]]
  14.  
  15. Clear[f]
  16. f = ConstantArray[0, scale];
  17. f[[1]] = N@HarmonicNumber[scale];
  18. Monitor[Do[
  19. f[[i]] = N@MangoldtLambda[i] + f[[i - 1]], {i, 2, scale}], i]
  20. xres = .002;
  21. x = Exp[Range[0, Log[scale], xres]];
  22. tmax = 60;
  23. tres = N[1/6];
  24. Monitor[errList =
  25. Table[t*(x^(-1/2 + I t).(f[[Floor[x]]] - x)), {t,
  26. Range[0, 60, tres]}];, t]
  27. g2 = ListLinePlot[Im[errList]/Length[x], DataRange -> {0, 60},
  28. PlotRange -> {-.02, 1.6}, PlotStyle -> Black]
  29.  
  30. scale = 600;
  31. Print["Counting to 60"]
  32. Monitor[g3 =
  33. ListLinePlot[
  34. Accumulate[+1.3 -
  35. Table[Re[
  36. Zeta[1/2 + I*t]*
  37. Total[Table[
  38. 1/2/Pi*1/n*
  39. Total[1/Divisors[n]^(1/2 + I*t - 1)*
  40. MoebiusMu[Divisors[n]]], {n, 1, scale}]]], {t, 0, 60,
  41. N[1/6]}]], DataRange -> {0, 60}, PlotRange -> {0, 80}],
  42. Floor[t]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement