SpheroidalS1[n,m,γ,z]
gives the radial spheroidal function
of the first kind.
SpheroidalS1
SpheroidalS1[n,m,γ,z]
gives the radial spheroidal function
of the first kind.
Details
- Mathematical function, suitable for both symbolic and numerical manipulation.
- The radial spheroidal functions satisfy the differential equation
with the spheroidal eigenvalue
given by SpheroidalEigenvalue[n,m,γ]. - The
are normalized according to the Meixner–Schäfke scheme. - SpheroidalS1 can be evaluated to arbitrary numerical precision.
- SpheroidalS1 automatically threads over lists. »
Examples
open all close allBasic Examples (4)
SpheroidalS1[1, 0, 2.5, 1]Plot over a subset of the reals:
Plot[SpheroidalS1[1, 0, 2, x], {x, -6, 6}]Plot over a subset of the complexes:
ComplexPlot3D[SpheroidalS1[1, 0, 2, z], {z, -1 - I, 1 + I}, PlotLegends -> Automatic]Series expansion at a singular point:
Series[SpheroidalS1[1, 1, π / 2, x], {x, 1, 6}]//Normal//FullSimplifyScope (21)
Numerical Evaluation (5)
SpheroidalS1[5, 1, .1, .3]SpheroidalS1[1, 0, 1.5, 2]N[SpheroidalS1[7, 5, 1, 2], 50]//ChopN[SpheroidalS1[2, 1, 1 / 3, -11], 20]The precision of the output tracks the precision of the input:
SpheroidalS1[2, 2, 1 / 3, 0.211111111111111111]//ChopSpheroidalS1[2, 0.211111111111111000111, 5, 6]N[SpheroidalS1[23, 5 - I, 2, I]]Evaluate efficiently at high precision:
SpheroidalS1[31, 5, 1, 2`100]//TimingSpheroidalS1[15, 5, 1, 5`1000];//TimingCompute the elementwise values of an array using automatic threading:
SpheroidalS1[1, 0, .2, {{1.3, .5}, {0, 0.2}}]Or compute the matrix SpheroidalS1 function using MatrixFunction:
MatrixFunction[SpheroidalS1[1, 0, .2, #]&, {{1.3, .5}, {0, 0.2}}]Specific Values (4)
Simple exact values are generated automatically:
Table[SpheroidalS1[1, 1, π / 2, x], {x, -1 / 2, 3 / 2}]Find the first positive maximum of SpheroidalS1[2,0,5,x]:
xmax = x /. FindRoot[D[SpheroidalS1[2, 0, 5, x], x] == 0, {x, .5}]//ChopPlot[SpheroidalS1[2, 0, 5, x], {x, 0, 3}, Epilog -> Style[Point[{xmax, SpheroidalS1[2, 0, 5, xmax ]}], PointSize[Large], Red]]SpheroidalS1 functions become elementary if m=1 and γ=n π/2 :
Table[SpheroidalS1[n, 1, n π / 2, x], {n, 1, 3}]//FullSimplifyTraditionalForm typesetting:
SpheroidalS1[n, m, c, η]//TraditionalFormVisualization (3)
Plot the SpheroidalS1 function for integer orders:
Plot[{SpheroidalS1[1, 0, 2, x], SpheroidalS1[2, 0, 2, x], SpheroidalS1[3, 0, 2, x]}, {x, -10, 10}]Plot the SpheroidalS1 function for noninteger parameters:
Plot[{SpheroidalS1[1, 0, 1 / 2, x], SpheroidalS1[2, 0, 1 / 2, x], SpheroidalS1[3, 0, 1 / 2, x]}, {x, -10, 10}]ComplexContourPlot[Re[SpheroidalS1[2, 0, 1, z]], {z, -4 - 6 I, 4 + 6 I}, Contours -> 20]ComplexContourPlot[Im[SpheroidalS1[2, 0, 1, z]], {z, -4 - 6 I, 4 + 6 I}, Contours -> 20]Function Properties (5)
SpheroidalS1 is not an analytic function:
FunctionAnalytic[SpheroidalS1[1, 1, π / 2, x], x]
has both singularities and discontinuities for
:
FunctionSingularities[SpheroidalS1[1, 2, π / 2, x], x]//QuietFunctionDiscontinuities[SpheroidalS1[1, 2, π / 2, x], x]//Quiet
is neither non-decreasing nor non-increasing:
FunctionMonotonicity[SpheroidalS1[1, 1, π / 2, x], x]FunctionInjective[SpheroidalS1[1, 1, π / 2, x], x]Plot[{SpheroidalS1[1, 1, π / 2, x], .3}, {x, -6, 6}]SpheroidalS1 is neither non-negative nor non-positive:
FunctionSign[SpheroidalS1[1, 1, π / 2, x], x]SpheroidalS1 is neither convex nor concave:
FunctionConvexity[SpheroidalS1[1, 1, π / 2, x], x]Differentiation (2)
First derivative with respect to z:
D[SpheroidalS1[n, m, γ, z], z]Higher derivatives with respect to z:
Table[D[SpheroidalS1[n, m, γ, z], {z, k}], {k, 1, 3}]//FullSimplifyPlot the higher derivatives with respect to z when n=10, m=2 and γ=1/3:
Plot[Evaluate[% /. { n -> 10, m -> 2, γ -> 1 / 3}], {z, -3, 3}, PlotLegends -> {"First Derivative", "Second Derivative", "Third Derivative"}]Series Expansions (2)
Find the Taylor expansion using Series:
Series[SpheroidalS1[n, m, γ, x], {x, 0, 3}]//NormalPlots of the first three approximations around
:
terms = Normal@Table[Series[SpheroidalS1[2, 0, 1, x], {x, 0, m}], {m, 1, 5, 2}]//N;
Plot[{SpheroidalS1[2, 0, 1, x], terms}, {x, 0, 4}, MaxRecursion -> 1]//QuietTaylor expansion at a generic point:
Series[SpheroidalS1[n, m, γ, x], {x, x0, 2}]// Normal//FullSimplifyApplications (4)
Spheroidal angular harmonics are eigenfunctions of the Sinc transform on the interval
:
Integrate[Sinc[c (t - x)]SpheroidalPS[2, 0, c, x], {x, -1, 1}]Plot[SpheroidalS1[2, 0, c, 1] ^ 2, {c, 1, 20}]Find resonant frequencies for the Dirichlet problem in the prolate spheroidal cavity:
Plot[SpheroidalS1[0, 0, k / (2 3^1 / 3) , 3], {k, 0, 20}, PlotPoints -> 20]Determine the first few frequencies:
FindRoot[SpheroidalS1[0, 0, k / (2 3^1 / 3) , 3], {k, 2, 5}, WorkingPrecision -> 50]FindRoot[SpheroidalS1[0, 0, k / (2 3^1 / 3) , 3], {k, 5, 8}, WorkingPrecision -> 50]Plot the prolate and oblate functions:
Plot[SpheroidalS1[2, 0, 10, r], {r, 1, 4}]Plot[SpheroidalS1[2, 0, -10 I, I r], {r, 0, 4}]Build a near-spherical approximation to
:
nsLimitS1[n_, m_, z_, {c_, ord_Integer}] :=
(Normal[Series[SpheroidalPS[n1, -m1, c, x], {c, 0, ord}] /. LegendreP[k_, -m1, 2, x] :> SphericalBesselJ[k, z]](1 - c ^ 2 / z ^ 2)^m / 2 / SpheroidalRadialFactor[n, -m, c]) /. {n1 -> n, m1 -> m}First few terms of the approximation:
nsLimitS1[n, m, z, {c, 3}]//TraditionalForm(Table[nsLimitS1[3, 1, 1`30, {c, k}], {k, 2, 10, 2}] - SpheroidalS1[3, 1, c, 1 / c]) /. {c -> 1 / 100}Tech Notes
Related Links
History
Text
Wolfram Research (2007), SpheroidalS1, Wolfram Language function, https://reference.wolfram.com/language/ref/SpheroidalS1.html.
CMS
Wolfram Language. 2007. "SpheroidalS1." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/SpheroidalS1.html.
APA
Wolfram Language. (2007). SpheroidalS1. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/SpheroidalS1.html
BibTeX
@misc{reference.wolfram_2026_spheroidals1, author="Wolfram Research", title="{SpheroidalS1}", year="2007", howpublished="\url{https://reference.wolfram.com/language/ref/SpheroidalS1.html}", note=[Accessed: 13-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_spheroidals1, organization={Wolfram Research}, title={SpheroidalS1}, year={2007}, url={https://reference.wolfram.com/language/ref/SpheroidalS1.html}, note=[Accessed: 13-June-2026]}