Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

NewtonZero never comes down

  • To: mathgroup at smc.vnet.net
  • Subject: [mg51416] NewtonZero never comes down
  • From: Roger Bagula <tftn at earthlink.net>
  • Date: Sat, 16 Oct 2004 04:21:07 -0400 (EDT)
  • Reply-to: tftn at earthlink.net
  • Sender: owner-wri-mathgroup at wolfram.com

I've tried several ways to find better values for
zeta zeros fast and easy...
FindRoot[] just doesn't give much accuracy.
NewtonZero[] tends to just spin it's wheels on my machine....
Here's what I came up with that gets about 20 places
and pretty fast relatively.
It's much better than FindRoot[] anyway!
The zeta zeros all seem to be transcendental irrational numbers
and are harder to calculate than most numbers in modern number theory.

(* finding Zeta Zero roots*)
(* using approximation from analyitic continuation on 0<=z<=1 strip 
where zeta dosesn't converge*)
Clear[RSzeta,t,g3,f,g,a]
RSzeta[t_]=RiemannSiegelZ[t]/Exp[I*RiemannSiegelTheta[t]]
Digits=10000
t=(2 Sqrt[2])*25*Pi^2*Sqrt[n]/(6*Sqrt[Digits])-16
(* gives the zeta zeros given in Jahnke-Emde Tables of Functions, Dover*)
g3=Delete[Union[Table[If[N[Abs[RSzeta[t]]]<0.01,N[t,20],0],{n,1,Digits}]],1]
Dimensions[g3]
ListPlot[g3,PlotJoined->True]
(* definition of Newton function*)
f[x_]:=N[x-Zeta[x]/Zeta'[x],20]
(* ten iterations toward the zero*)
g[x_]:=NestList[f,1/2+I*x,10][[10]]
(* values starting at the continuation values approximated*)
a=Table[g[g3[[n]]],{n,1,Dimensions[g3][[1]]}]
ListPlot[Im[a],PlotJoined->True]

Respectfully, Roger L. Bagula

tftn at earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
alternative email: rlbtftn at netscape.net
URL :  http://home.earthlink.net/~tftn


  • Prev by Date: Re: Mathematica Symbolic Toolbox for MATLAB--Version 2.0
  • Next by Date: Re: Re: RandomReplacement
  • Previous by thread: Re: LegendreP (Symbolic) is different in Mathematica5 than previous versions (M4, M3 ..)
  • Next by thread: Need code to calculate the Lower Envelope for a set of (non collinear) points.