RootFinding
Analytic
zeros of an analytic function in one variable
AnalyticZerosFound
Calling Sequence
Parameters
Description
Examples
Analytic(f, z, a+c*I..b+d*I, ...)
Analytic(f, z, re=a..b, im=c..d, ...)
AnalyticZerosFound()
f
-
an analytic function of z (or an equation defining such)
z
(optional) an unknown
a, b, c, d
real constants
For an analytic function f of z, Analytic attempts to find all complex zeros of f⁡z within the rectangular region a<=Re⁡z<=b, c<=Im⁡z<=d in the complex plane. It may also find some zeros outside but close to the boundary of the region.
Setting of infolevel[RootFinding:-Analytic] to be between 1 and 7 will result in detailed information concerning the solving process being displayed.
The AnalyticZerosFound() calling sequence returns a sequence of the zeros which have been located. These may be accessed after Analytic returns, or if its computation is interrupted.
Options
The remaining arguments are interpreted as options.
digits=n
Set Digits := n; for the purpose of the calculation. Thus the accuracy of the zeros computed is usually less than n digits since only n digits are used for the calculation. The default is the value of Digits already set by the user; the minimum is 5. If hardware floats are in use, this default may be increased to a maximum of trunc(evalhf(Digits)).
iterations=n
The number of iterations of Newton's method to be applied for each starting point. The default is 50.
continue
Continue with an interrupted calculation, remembering already found zeros. This option may also be used to change the boundaries of the region without losing already found zeros. This option may only be used with the same function and variable as the original call.
plot
Return a plot of the zeros instead of the zeros. Note, the zeros are still available to the user through the command RootFinding:-AnalyticZerosFound(), and are available until the routine RootFinding:-Analytic is once again used.
modulo=a+I⁢b and modulo_s=a+I⁢b
When plotting, the zeros will be reduced modulo a in the real direction and modulo b in the imaginary direction. With the option 'modulo' they will be reduced to the region 0<=Re⁡z<=a and 0<=Im⁡z<=b; With 'modulo_s' they will be reduced to the region −a2<=Re⁡z<=a2, −b2<=Im⁡z<=b2.
extra
For plotting purposes, remember the rectangular subregions in which the zeros were isolated, as well as the traces of the Newton iterations. These will appear when the zeros are plotted.
viewall
For plotting purposes, view the entire specified region. Do not restrict the viewing region to the area where zeros are found.
with⁡RootFinding
Analytic,AnalyticZerosFound,BivariatePolynomial,EnclosingBox,EvaluateAtRoot,HasRealRoots,Homotopy,Isolate,NextZero,Parametric,RefineRoot,WitnessPoints
f≔tan⁡sin⁡x−1
Analytic⁡f,x,−I2..1+I
0.903339110766515
g≔23⁢x5+105⁢x4−10⁢x2+17⁢x
Analytic⁡g,x,re=−5..1,im=−1..1
0.+0.⁢I,−0.637181318531050,0.304066454284907−0.404061905751759⁢I,0.304066454284907+0.404061905751759⁢I,−4.53616898134312
h≔10−ln⁡v+v2−112−ln⁡3+32−112
h≔10−ln⁡v+v2−1+ln⁡3+8
Analytic⁡h,v,re=2000..100000,im=−100000..100000,digits=32
64189.825354267506015885528786917
Analytic⁡sin⁡x,x,re=−10..10,im=−10..10,digits=10
0.+0.⁢I,−6.283185305,−3.141592654,−9.424777960,6.283185305,3.141592654,9.424777960,0.+0.⁢I,0.+0.⁢I
f0≔sin⁡x−2x
Analytic⁡f0,x,re=−10..10,im=−10..10,digits=10
−9.426231485,−3.247111097,−6.270228880,0.6634376855−1.171931328⁢I,4.980776967−4.145343143⁢I,9.224685705−7.087212710⁢I,0.6634376855+1.171931328⁢I,4.980776967+4.145343143⁢I,9.224685705+7.087212710⁢I
numlist≔23,1+I5,−12−I2,7+I:
f1≔mul⁡z−z0,z0=numlist
f1≔z−23⁢z−15−I5⁢z+12+I2⁢z−7−I
Analytic⁡f1,z,re=−2..10,im=−2..2
0.666666666666665,0.200000000000000+0.200000000000000⁢I,7.00000000000000+I,−0.50000000000000−0.50000000000000⁢I
The following examples may be executed in Maple.
f2≔add⁡z0z,z0=numlist
RootFinding:-Analytic⁡f2,z,re=−10..10,im=−10..10,plot
RootFinding:-AnalyticZerosFound⁡
f3≔add⁡1nz,n=1..10
RootFinding:-Analytic⁡f3,z,re=−5..2,im=−1..50,plot,extra
f4≔−100⁢I⁢x2+78⁢x4+39⁢I⁢x4−94⁢x3−17⁢x4+19⁢x−73⁢I⁢x−25−4⁢I−59⁢x5
zeros≔RootFinding:-Analytic⁡f4,x,re=−10..10,im=−10..10
plotscomplexplot⁡zeros,style=point,axes=boxed
See Also
fsolve
Download Help Document