Advertisement
MatsGranvik

ZetaZeros

Jul 1st, 2012 (edited)
212
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.77 KB | None | 0 0
  1. In[543]:= Clear[n, s, a];
  2. s = 1/2;
  3. NSum[1/(n + 0)^s - 1/(n + 1)^s - 2/(n + 2)^s - 1/(n + 3)^s +
  4. 1/(n + 4)^s + 2/(n + 5)^s, {n, 1, Infinity}, WorkingPrecision -> 400]
  5. N[2.549127729379167407581967029267929878/(Im[ZetaZero[18]]/2), 20]^-1
  6. N[Im[ZetaZero[1]], 20]
  7. N[2.549127729379167407581967029267929878/(Im[ZetaZero[33]]/2), 20]^-1
  8. N[Im[ZetaZero[2]], 20]
  9. N[2.549127729379167407581967029267929878/(Im[ZetaZero[42]]/2), 20]^-1
  10. N[Im[ZetaZero[3]], 20]
  11.  
  12. During evaluation of In[543]:= NIntegrate::ncvb: NIntegrate failed to converge to prescribed accuracy after 9 recursive bisections in n near {n} = {<<406>>}. NIntegrate obtained -<<406>> and <<411>> for the integral and error estimates. >>
  13.  
  14. Out[545]= \
  15. -2.5491277293791674075819670292679298788388868361501383611646759945550\
  16. 2114281280979141892279032135882066112601786029735792734346777007416385\
  17. 1031714272022490694260100601769646158455686154700856477098022672879033\
  18. 4388858289379494842246720038280526184788343913638074859695887808089231\
  19. 7899988144451825952399138069988923352228536471316749487823829937101794\
  20. 3934028300079177278303471836107425673221464057551975
  21.  
  22. Out[546]= 14.135650568603255663
  23. Out[547]= 14.134725141734693790
  24. Out[548]= 21.020643640006420723
  25. Out[549]= 21.022039638771554993
  26. Out[550]= 25.011827067342131577
  27. Out[551]= 25.010857580145688763
  28.  
  29.  
  30. "start"
  31. Clear[s];
  32. s /. FindRoot[-(Zeta[s]*Zeta[s]/(2*Zeta[2*s - 1]) + (1 - 1/2^(s - 1))*
  33. Zeta[s] + Zeta[s]/(2*Zeta[s - 1]) + 1/Zeta[s - 1]) ==
  34. Im[ZetaZero[1]], {s, -1/3}, WorkingPrecision -> 50]
  35.  
  36. -0.14134053532941803058116067801572629
  37. "end"
  38.  
  39. "start"
  40. Clear[s];
  41. -Limit[Zeta[s]*Zeta[s]/(2*Zeta[2*s - 1]) + (1 - 1/2^(s - 1))*Zeta[s] +
  42. Zeta[s]/(2*Zeta[s - 1]) + 1/Zeta[s - 1], s -> 1]
  43. N[%, 30]
  44. "end"
  45.  
  46. "Output:"
  47. 2 + EulerGamma - Log[4 \[Pi]]
  48.  
  49.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement