Advertisement
MatsGranvik

"Dirichlet series" for log(log(n))

Oct 22nd, 2022 (edited)
86
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.80 KB | None | 0 0
  1. (*Matheamtica start*)
  2. LL = 19;
  3. Monitor[TableForm[numerator = Table[Clear[nn, loglog, s];
  4. $MaxPrecision = 100;
  5. s = s;
  6. nn = LL;
  7. loglog = L;
  8. TableForm[
  9. Z = Table[
  10. Table[If[Mod[n, k] == 0, 1/(n/k)^s, 0], {k, 1, nn}], {n, 1,
  11. nn}]];
  12. TableForm[
  13. A = Table[
  14. Table[If[n/k == loglog, -loglog/(n/k)^s, 0], {k, 1, nn}], {n,
  15. 1, nn}] + IdentityMatrix[nn]];
  16. (*TableForm[Z.A];*)
  17. generatingFunction = -(Zeta[
  18. s] ((IdentityMatrix[nn] - MatrixPower[(Z.A), s - 1])[[All,
  19. 1]]));
  20. limit = Limit[generatingFunction, s -> 1];
  21. accumulate = Accumulate[limit];
  22. cesaro = N[Total[accumulate]/Length[accumulate], 30];
  23. Numerator[limit], {L, 1, LL}]], L]
  24.  
  25. LL = 19;
  26. Monitor[TableForm[denominator = Table[Clear[nn, loglog, s];
  27. $MaxPrecision = 100;
  28. s = s;
  29. nn = LL;
  30. loglog = L;
  31. TableForm[
  32. Z = Table[
  33. Table[If[Mod[n, k] == 0, 1/(n/k)^s, 0], {k, 1, nn}], {n, 1,
  34. nn}]];
  35. TableForm[
  36. A = Table[
  37. Table[If[n/k == loglog, -loglog/(n/k)^s, 0], {k, 1, nn}], {n,
  38. 1, nn}] + IdentityMatrix[nn]];
  39. (*TableForm[Z.A];*)
  40. generatingFunction = -(Zeta[
  41. s] ((IdentityMatrix[nn] - MatrixPower[(Z.A), s - 1])[[All,
  42. 1]]));
  43. limit = Limit[generatingFunction, s -> 1];
  44. accumulate = Accumulate[limit];
  45. cesaro = N[Total[accumulate]/Length[accumulate], 30];
  46. Denominator[limit], {L, 1, LL}]], L]
  47. TableForm[ratio = numerator/denominator]
  48. TableForm[
  49. logTable =
  50. Table[Table[If[Mod[k, n] == 0, 1 - n, 1]/k, {k, 1, nn}], {n, 1, nn}]]
  51. TableForm[ratio = Range[Length[denominator]]*numerator/denominator]
  52. TableForm[ratio + logTable]
  53. TableForm[Floor[ratio + logTable]]
  54. (*end*)
  55.  
  56. (*start*)
  57. LL = 19;
  58. Monitor[TableForm[numerator = Table[Clear[nn, loglog, s];
  59. $MaxPrecision = 100;
  60. s = s;
  61. nn = LL;
  62. loglog = L;
  63. TableForm[
  64. Z = Table[
  65. Table[If[Mod[n, k] == 0, 1/(n/k)^s, 0], {k, 1, nn}], {n, 1,
  66. nn}]];
  67. TableForm[
  68. A = Table[
  69. Table[If[n/k == loglog, -loglog/(n/k)^s, 0], {k, 1, nn}], {n,
  70. 1, nn}] + IdentityMatrix[nn]];
  71. (*TableForm[Z.A];*)
  72. generatingFunction = -(Zeta[
  73. s] ((IdentityMatrix[nn] - MatrixPower[(Z.A), s - 1])[[All,
  74. 1]]));
  75. limit = Limit[generatingFunction, s -> 1];
  76. accumulate = Accumulate[limit];
  77. cesaro = N[Total[accumulate]/Length[accumulate], 30];
  78. Numerator[limit], {L, 1, LL}]], L]
  79.  
  80. LL = 19;
  81. Monitor[TableForm[denominator = Table[Clear[nn, loglog, s];
  82. $MaxPrecision = 100;
  83. s = s;
  84. nn = LL;
  85. loglog = L;
  86. TableForm[
  87. Z = Table[
  88. Table[If[Mod[n, k] == 0, 1/(n/k)^s, 0], {k, 1, nn}], {n, 1,
  89. nn}]];
  90. TableForm[
  91. A = Table[
  92. Table[If[n/k == loglog, -loglog/(n/k)^s, 0], {k, 1, nn}], {n,
  93. 1, nn}] + IdentityMatrix[nn]];
  94. (*TableForm[Z.A];*)
  95. generatingFunction = -(Zeta[
  96. s] ((IdentityMatrix[nn] - MatrixPower[(Z.A), s - 1])[[All,
  97. 1]]));
  98. limit = Limit[generatingFunction, s -> 1];
  99. accumulate = Accumulate[limit];
  100. cesaro = N[Total[accumulate]/Length[accumulate], 30];
  101. Denominator[limit], {L, 1, LL}]], L]
  102. TableForm[ratio = numerator/denominator]
  103. TableForm[
  104. logTable =
  105. Table[Table[If[Mod[k, n] == 0, 1 - n, 1]/k, {k, 1, nn}], {n, 1, nn}]]
  106. TableForm[ratio = Range[Length[denominator]]*numerator/denominator]
  107. TableForm[ratio + logTable]
  108. TableForm[ttable = (ratio + logTable)*3]
  109. TableForm[
  110. divisorTable =
  111. Table[Table[If[Mod[k, n] == 0, n, 0], {k, 1, nn}], {n, 1, nn}]]
  112. TableForm[ttable - divisorTable]
  113. TableForm[endTable = (ttable - Transpose[divisorTable])]
  114. TableForm[Floor[1/2 + ttable - Transpose[divisorTable]]]
  115. endTable[[1, 1]] = 0;
  116. Total[endTable]
  117. Total[Transpose[endTable]]
  118. Total[Total[endTable]]
  119. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement