NevilleThetaS[z,m]
gives the Neville theta function
.
NevilleThetaS
NevilleThetaS[z,m]
gives the Neville theta function
.
Details
- Mathematical function, suitable for both symbolic and numerical manipulation.
-
- NevilleThetaS[z,m] is a meromorphic function of
and has a complicated branch cut structure in the complex
plane. - For certain special arguments, NevilleThetaS automatically evaluates to exact values.
- NevilleThetaS can be evaluated to arbitrary numerical precision.
- NevilleThetaS automatically threads over lists.
Examples
open all close allBasic Examples (3)
Scope (27)
Numerical Evaluation (6)
NevilleThetaS[-5., 1]NevilleThetaS[Pi / 6., 0]N[NevilleThetaS[5 / 7, 2 / 3], 50]N[NevilleThetaS[2Pi, 7 / 11], 50]The precision of the output tracks the precision of the input:
NevilleThetaS[2.211111100011100001, 2 / 3]N[NevilleThetaS[1 + I, 5 - I]]Evaluate efficiently at high precision:
NevilleThetaS[2, 1 / 2`100]//TimingNevilleThetaS[2, 17`1000];//TimingCompute average-case statistical intervals using Around:
NevilleThetaS[Around[2, 0.01], 1 / 2]Compute the elementwise values of an array:
NevilleThetaS[π I / 3, {{0, 1}, {1, 0}}]Or compute the matrix NevilleThetaS function using MatrixFunction:
MatrixFunction[NevilleThetaS[π I / 3, #]&, {{0, 1}, {1, 0}}]//NSpecific Values (3)
Values at corners of the fundamental cell:
Table[NevilleThetaS[u1 EllipticK[m] + u2 I EllipticK[1 - m], m], {u1, {0, 1}}, {u2, {0, 1}}]NevilleThetaS for special values of elliptic parameter:
NevilleThetaS[x, 0]NevilleThetaS[x, 1]Find the first positive maximum of NevilleThetaS[x,1/2]:
xmax = x /. FindRoot[D[NevilleThetaS[x, 1 / 2 ], x] == 0, {x, 2}]Plot[NevilleThetaS[x, 1 / 2 ], {x, -9, 9}, Epilog -> Style[Point[{xmax, NevilleThetaS[xmax, 1 / 2 ]}], PointSize[Large], Red]]Visualization (3)
Plot the NevilleThetaS functions for various values of the parameter:
Plot[{NevilleThetaS[x, -1 / 2], NevilleThetaS[x, 1 / 2], NevilleThetaS[x, 1 / 3]}, {x, -5, 5}]Plot NevilleThetaS as a function of its parameter
:
Plot[{NevilleThetaS[1, m], NevilleThetaS[2, m], NevilleThetaS[3, m]}, {m, -10, 5}, PlotRange -> Automatic]ComplexContourPlot[Re[NevilleThetaS[z, 1 / 2]], {z, -2 - 2I, 2 + 2I}, IconizedObject[«PlotOptions»]]ComplexContourPlot[Im[NevilleThetaS[z, 1 / 2]], {z, -2 - 2I, 2 + 2I}, IconizedObject[«PlotOptions»]]Function Properties (11)
The real domain of NevilleThetaS:
FunctionDomain[NevilleThetaS[x, 0], x]The complex domain of NevilleThetaS:
FunctionDomain[NevilleThetaS[z, 0], z, Complexes]FunctionRange[NevilleThetaS[x, 0], x, y]//QuietFunctionRange[NevilleThetaS[x, 1], x, y]//QuietNevilleThetaS threads elementwise over lists:
NevilleThetaS[{1, 2, 3, 4}, 0.4]
is an analytic function of
for
:
FunctionAnalytic[NevilleThetaS[x, m], x]
is neither non-decreasing nor non-increasing:
FunctionMonotonicity[Subscript[ϑ, s](x❘(1/3)), x]FunctionInjective[Subscript[ϑ, s](x❘(1/3)), x]Plot[{NevilleThetaS[x, 1 / 3], .5}, {x, -7, 7}]FunctionSurjective[Subscript[ϑ, s](x❘(1/3)), x]Plot[{Subscript[ϑ, s](x❘(1/3)), 1.3}, {x, -10, 10}]
is neither non-negative nor non-positive for noninteger
:
Table[FunctionSign[Subscript[ϑ, s](x❘(1/n)), x], {n, 4}]
has no singularities or discontinuities for noninteger
:
Table[FunctionSingularities[Subscript[ϑ, d](x❘(1/n)), x], {n, 4}]Table[FunctionDiscontinuities[Subscript[ϑ, d](x❘(1/n)), x], {n, 4}]
is affine only for
and otherwise it is neither convex nor concave:
Table[FunctionConvexity[Subscript[ϑ, d](x❘1 / n), x], {n, 5}]TraditionalForm formatting:
NevilleThetaS[z, m]//TraditionalFormDifferentiation (2)
D[NevilleThetaS[u, m], u]D[NevilleThetaS[u, m], m]Table[D[NevilleThetaS[u, 1 / 3], {u, k}], {k, 1, 3}]//Simplify//TraditionalFormPlot the higher-order derivatives:
Plot[%, {u, -10, 10}, PlotLegends -> {"First Derivative", "Second Derivative", "Third Derivative"}]Series Expansions (2)
Find the Taylor expansion using Series:
Series[NevilleThetaS[x, m], {x, 1, 2}]Plots of the first three approximations around
:
terms = Normal@Table[Series[NevilleThetaS[x, 1 / 2], {x, 1, m}], {m, 1, 3, 1}]//N;
Plot[{NevilleThetaS[x, 1 / 2], terms}, {x, -10, 10}]The Taylor expansion for small elliptic parameter
:
Series[NevilleThetaS[x, m], {m, 0, 2}]Series[NevilleThetaS[x, m], {m, 1, 1}]Generalizations & Extensions (1)
NevilleThetaS can be applied to power series:
NevilleThetaS[Sin[z] + O[z]^4, m]Applications (7)
Plot over the arguments' plane:
Plot3D[NevilleThetaS[x, m], {x, -3 / 2Pi, 3 / 2Pi}, {m, -1, 1}]Conformal map from a unit triangle to the unit disk:
w[z_] := With[{ζ = (z Gamma[(1/3)]^3/π 2^1 / 3 3^1 / 4), m = Sin[(π/12)]^2}, (2^1 / 3 (1 + Sqrt[3]) (NevilleThetaC[ζ, m] - NevilleThetaN[ζ, m]) ((-2 + Sqrt[3]) NevilleThetaC[ζ, m] - NevilleThetaN[ζ, m])) / (NevilleThetaC[ζ, m]^2 - 2 NevilleThetaC[ζ, m] NevilleThetaN[ζ, m] + NevilleThetaN[ζ, m]^2 + 2 3^1 / 4 NevilleThetaD[ζ, m] NevilleThetaS[ζ, m])]Show points before and after the map:
trianglePoints = Flatten[Table[α Table[N[#1 + j / 16(#2 - #1)], {j, 0, 16}]&@@@({{2, -1 + I Sqrt[3]}, {-1 + I Sqrt[3], -1 - I Sqrt[3]}, {-1 - I Sqrt[3], 2}} / (2Sqrt[3])), {α, 1 / 8, 1, 1 / 8}]];Row[{Graphics[Point[{Re[#], Im[#]}]& /@ trianglePoints, Frame -> True, ImageSize -> Small],
Graphics[Point[{Re[#], Im[#]}]& /@ (w /@ trianglePoints), Frame -> True, ImageSize -> Small]}, Spacer[20]]Uniformization of a Fermat cubic
:
{a[u_], b[u_]} = With[{m = (2 - Sqrt[3]) / 4}, {-1 + (4 3^1 / 4 NevilleThetaD[u, m] NevilleThetaS[u, m]) / ((NevilleThetaC[u, m] - NevilleThetaN[u, m])^2 + 2 3^1 / 4 NevilleThetaD[u, m] NevilleThetaS[u, m]), -(2^1 / 3 (NevilleThetaC[u, m] - NevilleThetaN[u, m]) ((-1 + Sqrt[3]) NevilleThetaC[u, m] + (1 + Sqrt[3]) NevilleThetaN[u, m])) / ((NevilleThetaC[u, m] - NevilleThetaN[u, m])^2 + 2 3^1 / 4 NevilleThetaD[u, m] NevilleThetaS[u, m])}];ListPlot[Table[{a[u], b[u]}, {u, -Pi, Pi, 0.1}], AspectRatio -> Automatic, ImageSize -> Small]Verify that points on the curve satisfy
:
Cases[%, Point[l__] :> Map[Total, l ^ 3], Infinity, 1]Current flow in a rectangular conducting sheet with voltage applied at a pair of opposite corners:
flow[m_] = NevilleThetaD[x, m] NevilleThetaD[y, 1 - m] NevilleThetaS[x, m] NevilleThetaS[y, 1 - m] / (NevilleThetaC[x, m] NevilleThetaC[y, 1 - m] NevilleThetaN[x, m] NevilleThetaN[y, 1 - m]);With[{m = 0.2},
ContourPlot[flow[m], {x, 0, EllipticK[m]}, {y, 0, EllipticK[1 - m]}, ContourShading -> False, Contours -> 50]]Parametrize a lemniscate by arc length:
lemniscate[s_] := With[{t = Sqrt[2]s}, {NevilleThetaC[t, (1/2)] NevilleThetaD[t, (1/2)], (NevilleThetaC[t, (1/2)] NevilleThetaS[t, (1/2)]/Sqrt[2])} / NevilleThetaN[t, (1/2)]^2]Show the classical and arc length parametrizations:
r[θ_] := If[Cos[2θ] > 0, Sqrt[Cos[2θ]], 0];Row[{ParametricPlot[{r[θ]Cos[θ], r[θ]Sin[θ]}, {θ, -Pi, Pi}, ImageSize -> Small],
ParametricPlot[Evaluate[lemniscate[s]], {s, 0, (Gamma[(1/4)]^2/Sqrt[2 π])}, ImageSize -> Small]}, Spacer[20]]Complex parametrization of a sphere:
ellipticSphere[{α_, β_}, m_] :=
{(Sqrt[m] NevilleThetaS[α, m] NevilleThetaS[β, m]/NevilleThetaN[α, m] NevilleThetaN[β, m]), (I Sqrt[m] NevilleThetaC[α, m] NevilleThetaC[β, m]/Sqrt[1 - m] NevilleThetaN[α, m] NevilleThetaN[β, m]), (NevilleThetaD[α, m] NevilleThetaD[β, m]/Sqrt[1 - m] NevilleThetaN[α, m] NevilleThetaN[β, m])}The square of all points on the complex sphere is 1:
N[ellipticSphere[{1 - 2I, 3 + 2I}, 1 / 2], 20]%.%Conformal map from an ellipse to the unit disk:
w[z_, {a_, b_}] := With[{m = InverseEllipticNomeQ[(a - b) ^ 2 / (a + b) ^ 2]}, (Sqrt[m]NevilleThetaS[(2ArcSin[z / Sqrt[a ^ 2 - b ^ 2]]
EllipticK[m]) / Pi, m]) / NevilleThetaN[
(2ArcSin[z / Sqrt[a ^ 2 - b ^ 2]]EllipticK[m]) / Pi, m]]ParametricPlot[{Re[w[ρ (2 Cos[φ] + I Sin[φ]), {2, 1}]], Im[w[ρ (2 Cos[φ] + I Sin[φ]), {2, 1}]]}, {ρ, 0, 1}, {φ, 0, 2Pi}, Mesh -> 15]Properties & Relations (4)
Basic simplifications are automatically carried out:
NevilleThetaS[-z, m]NevilleThetaS[z + 2 EllipticK[m], m]All Neville theta functions are a multiple of shifted NevilleThetaS:
NevilleThetaS[z + EllipticK[m], m]NevilleThetaS[z + I EllipticK[1 - m], m]NevilleThetaS[z + EllipticK[m] + I EllipticK[1 - m], m]//SimplifyUse FullSimplify for expressions containing Neville theta functions:
NevilleThetaS[z, m] == (NevilleThetaD[z, m]/JacobiDS[z, m])//FullSimplifyNumerically find a root of a transcendental equation:
FindRoot[NevilleThetaS[z, (1/3)]^3 + NevilleThetaS[z, (1/2)] + z == 3, {z, 1}]Related Guides
Related Links
History
Introduced in 1996 (3.0)
Text
Wolfram Research (1996), NevilleThetaS, Wolfram Language function, https://reference.wolfram.com/language/ref/NevilleThetaS.html.
CMS
Wolfram Language. 1996. "NevilleThetaS." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/NevilleThetaS.html.
APA
Wolfram Language. (1996). NevilleThetaS. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/NevilleThetaS.html
BibTeX
@misc{reference.wolfram_2026_nevillethetas, author="Wolfram Research", title="{NevilleThetaS}", year="1996", howpublished="\url{https://reference.wolfram.com/language/ref/NevilleThetaS.html}", note=[Accessed: 13-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_nevillethetas, organization={Wolfram Research}, title={NevilleThetaS}, year={1996}, url={https://reference.wolfram.com/language/ref/NevilleThetaS.html}, note=[Accessed: 13-June-2026]}