MatsGranvik

Mertens function

Apr 10th, 2014
385
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.57 KB | None | 0 0
  1. (*Mathematica Mats Granvik 7.4.2014*)
  2. Clear[A, t, n, k, a, nn];
  3. nn = 11;
  4. a = Table[StringJoin[{"x", ToString[n - 1]}], {n, 1, nn}]
  5. Length[a]
  6. t[n_, 1] = If[n <= 3, 1, a[[n]]];
  7. t[n_, k_] :=
  8. t[n, k] =
  9. If[n >= k,
  10. If[n <= 3, 1,
  11. If[And[k > 1],
  12. If[Or[k == 2, k == 3],
  13. t[n, k - 1] - Sum[t[n - i, k], {i, 1, k - 1}],
  14. If[k >= 4,
  15. Sum[t[n - i, k - 1], {i, 1, k - 2}] -
  16. Sum[t[n - i, k], {i, 1, k - 1}], 0], 0], 0, 0], 0], 0];
  17. A = Table[Table[t[n, k], {k, 1, nn}], {n, 1, nn}];
  18. MatrixForm[A]
  19. Inverse[A][[2 ;; nn, 1]]
Advertisement
Add Comment
Please, Sign In to add comment