MatsGranvik

von Mangoldt Dirac comb mimic 0.500000000000000000000000000000 + 0.220635600152651593396456432118 I

Jun 28th, 2022 (edited)
96
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.99 KB | None | 0 0
  1. TableForm[Table[Clear[nn, t, n, k, M, x];
  2. nn = mm;
  3. t[n_, 1] = 1;
  4. t[1, k_] = 1;
  5. t[n_, k_] :=
  6. t[n, k] =
  7. If[n < k,
  8. If[And[n > 1, k > 1],
  9. I*Exp[-I*x*Pi] - Sum[t[k - i, n], {i, 1, n - 1}], 0],
  10. If[And[n > 1, k > 1],
  11. I*Exp[-I*x*Pi] - Sum[t[n - i, k], {i, 1, k - 1}], 0]];
  12. M = Table[Table[t[n, k], {k, 1, nn}], {n, 1, nn}];
  13. MatrixForm[M];
  14. Det[M]; x /. Solve[Det[M] == 0, x], {mm, 2, 12}]];
  15. MatrixForm[M];
  16. N[%%, 30]
  17.  
  18.  
  19.  
  20. (*Mathematica 8*)
  21. Clear[x, xx]
  22. scale = 200;
  23. xres = .001;
  24. x = Exp[Range[0, Log[scale], xres]];
  25. xx = Flatten[{0, Differences[Floor[Exp[Range[0, Log[scale], xres]]]]}];
  26. yy = Accumulate[xx]*xx;
  27. zz = Table[
  28. If[yy[[i]] == 0, 0, MangoldtLambda[yy[[i]]]], {i, 1, Length[yy]}];
  29. ListLinePlot[zz*x^(-1), PlotRange -> {-0.1, 0.8}, ImageSize -> Large]
  30. ListLinePlot[Accumulate[zz*x^(-1)], ImageSize -> Large]
  31. Graphics[Line[{{Log[2], Min[Accumulate[zz*x^(-1)]]}, {Log[scale],
  32. Max[Accumulate[zz*x^(-1)]]}}]]
  33. (*end*)
Add Comment
Please, Sign In to add comment