Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- (*start*)
- Clear[nn, s, c, z, x];
- nn = 16;
- s = 10000;
- c = 10000;
- z = 10000;
- x = 10000;
- (*Table[Limit[Zeta[ss] \
- Total[1/Divisors[n]^(ss-z-x)*MoebiusMu[Divisors[n]]]/n^c,ss->s],{n,1,\
- nn}];*)
- A = Table[
- Table[If[Mod[n, k] == 0, MoebiusMu[k]*k^z/n^s, 0], {k, 1, nn}], {n,
- 1, nn}];
- B = Table[
- Table[If[Mod[k, n] == 0, n^x/k^c, 0], {k, 1, nn}], {n, 1, nn}];
- TableForm[T = Chop[N[A.B]]]
- "Double sum:"
- Sum[Sum[N[T, 20][[n, k]], {k, 1, nn}], {n, 1, nn}]
- "Generating function for the double sum:"
- N[Zeta[s]*Zeta[c]/Zeta[s + c - z - x], 20]
- (*end*)
Advertisement
Add Comment
Please, Sign In to add comment