WhittakerM[k,m,z]
gives the Whittaker function
.
WhittakerM
WhittakerM[k,m,z]
gives the Whittaker function
.
Details
- Mathematical function, suitable for both symbolic and numerical manipulation.
- WhittakerM is related to the Kummer confluent hypergeometric function by
.
vanishes at
for
. - For certain special arguments, WhittakerM automatically evaluates to exact values.
- WhittakerM can be evaluated to arbitrary numerical precision.
- WhittakerM automatically threads over lists.
- WhittakerM[k,m,z] has a branch cut discontinuity in the complex
plane running from
to
. - WhittakerM can be used with Interval and CenteredInterval objects. »
Examples
open all close allBasic Examples (6)
WhittakerM[2, 3, 4.5]Use FunctionExpand to expand in terms of hypergeometric functions:
FunctionExpand[WhittakerM[k, m, x]]Plot over a subset of the reals
:
Plot[WhittakerM[2, 1 / 2, x], {x, 0, 15}]Plot over a subset of the complexes:
ComplexPlot3D[WhittakerM[2, 1 / 2, z], {z, -2 - 2I, 2 + 2I}, PlotLegends -> Automatic]Series expansion at the origin:
Series[WhittakerM[2, 1 / 2, x], {x, 0, 6}]Series expansion at Infinity:
Series[WhittakerM[2, 1 / 2, x], {x, ∞, 5}]//NormalScope (35)
Numerical Evaluation (6)
WhittakerM[2, 3, 1.7]WhittakerM[2, 0.5, 0]N[WhittakerM[2, 1 / 3, 1 / 7], 100]The precision of the output tracks the precision of the input:
WhittakerM[2, 1.33333333333333333333, 1]N[WhittakerM[4 / 5, 2 - I, 2]]Evaluate efficiently at high precision:
WhittakerM[2, 1 / 3, 1 / 7`100]//TimingWhittakerM[2, 7 / 3, 4 / 7`10000];//TimingWhittakerM can be used with Interval and CenteredInterval objects:
WhittakerM[0.3, 0.4, Interval[{0.5, 0.6}]]WhittakerM[1, 2, CenteredInterval[3, 1 / 100]]Compute the elementwise values of an array:
WhittakerM[1, 3.1, {{2, 3}, {1, 2}}]Or compute the matrix WhittakerM function using MatrixFunction:
MatrixFunction[WhittakerM[1.3, 3, #]&, {{2, 0}, {1.2, -2}}]Specific Values (7)
WhittakerM for symbolic parameters:
WhittakerM[k, m, 2]//FunctionExpandWhittakerM[1, 1, 0]WhittakerM[0, 0, 0]Find the first positive maximum of WhittakerM[5,1/2,x]:
xmax = x /. FindRoot[D[WhittakerM[5, 1 / 2, x ], x] == 0, {x, 0.5}]Plot[WhittakerM[5, 1 / 2, x ], {x, 0, 10}, Epilog -> Style[Point[{xmax, WhittakerM[5, 1 / 2, xmax ]}], PointSize[Large], Red]]Compute the associated WhittakerM[3,1/2,x] function:
WhittakerM[3, 1 / 2, x]//FunctionExpandCompute the associated WhittakerM function for half-integer parameters:
WhittakerM[1 / 2, 1 / 2, x]//FunctionExpandDifferent cases of WhittakerM give different symbolic forms:
Table[WhittakerM[k, m, x], {k, {-1 / 2, 1 / 2}}, {m, {-1, 0, 1}}]//FunctionExpandWhittakerM threads elementwise over lists:
WhittakerM[2, 3, {1.5, 2.5, 3.5}]Visualization (3)
Plot the WhittakerM function for various orders:
Plot[{WhittakerM[1, 1 / 2, x], WhittakerM[2, 1 / 2, x], WhittakerM[3, 1 / 2, x], WhittakerM[4, 1 / 2, x]}, {x, 0, 3}]ComplexContourPlot[Re[WhittakerM[10, 1 / 2, z]], {z, -2 - 2I, 2 + 2I}, Contours -> 24]ComplexContourPlot[Im[WhittakerM[10, 1 / 2, z]], {z, -2 - 2I, 2 + 2I}, Contours -> 24]Plot as real parts of two parameters vary:
Plot3D[Re[WhittakerM[k, 1 / 2, z]], {k, 0, 3}, {z, 0, 2}, PlotRange -> All]Function Properties (11)
FunctionDomain[WhittakerM[2, 0, x], x]Complex domain of WhittakerM:
FunctionDomain[WhittakerM[k, m, z], z, Complexes]FunctionRange[WhittakerM[2, 0, z], z, y]//QuietWhittakerM may reduce to simpler functions:
WhittakerM[-1 / 4, 1 / 4, x ^ 2]//FunctionExpandWhittakerM[n / 2 + 1 / 4, -1 / 4, 1 / 2z ^ 2]//FunctionExpand
is not an analytic function of
for integer values of
:
FunctionAnalytic[WhittakerM[2, 0, x], x]FunctionMeromorphic[WhittakerM[2, 0, x], x]It is analytic for other values of
:
FunctionAnalytic[WhittakerM[3 / 2, 1 / 2, x], x, ℂ]
is neither non-decreasing nor non-increasing:
FunctionMonotonicity[WhittakerM[2, 0, x], x]FunctionInjective[WhittakerM[2, 0, x], x]Plot[{WhittakerM[2, 0, x], .2}, {x, 0, 6}]FunctionSurjective[WhittakerM[2, 1 / 2, x], x]Plot[{WhittakerM[2, 1 / 2, x], 1}, {x, -3, 6}]
is neither non-negative nor non-positive on its real domain:
FunctionSign[{WhittakerM[2, 0, x], x >= 0}, x]WhittakerM has both singularity and discontinuity in (-∞,0]:
FunctionSingularities[WhittakerM[a, m, x], x]FunctionDiscontinuities[WhittakerM[a, m, x], x]
is neither convex nor concave on its real domain:
FunctionConvexity[{WhittakerM[2, 0, x], x >= 0}, x]TraditionalForm formatting:
WhittakerM[m, k, x]//TraditionalFormDifferentiation (3)
First derivative with respect to z:
D[WhittakerM[k, m, z] , z]Higher derivatives with respect to z when k=1/3 and m=1/2:
Table[D[WhittakerM[1 / 3, 1 / 2, z], {z, n}], {n, 1, 3}]//FullSimplifyPlot the higher derivatives with respect to z when k=1/3 and m=1/2:
Plot[%, {z, -10, 10}, PlotLegends -> {"First Derivative", "Second Derivative", "Third Derivative"}]Formula for the ![]()
derivative with respect to z:
D[WhittakerM[k, m, z], {z, k}]// FullSimplifySeries Expansions (5)
Find the Taylor expansion using Series:
Series[WhittakerM[k, m, x], {x, 0, 4}]Plots of the first three approximations around
:
terms = Normal@Table[Series[WhittakerM[2, 1 / 2, x], {x, 0, m}], {m, 1, 5, 2}];
Plot[{WhittakerM[2, 1 / 2, x], terms}, {x, 0, 10}, PlotRange -> {-2, 2}]General term in the series expansion using SeriesCoefficient:
SeriesCoefficient[WhittakerM[k, m, x], {x, 1, n}]Find the series expansion at Infinity:
Series[WhittakerM[k, m, x], {x, Infinity, 1}]Find series expansion for an arbitrary symbolic direction
:
Series[WhittakerM[k, m, x], {x, DirectedInfinity[z], 1}, Assumptions -> x > 0]// FullSimplifyTaylor expansion at a generic point:
Series[WhittakerM[k, m, x], {x, x0, 2}]// FullSimplifyApplications (2)
The bound-state Coulomb eigenfunction in parabolic coordinates:
Subscript[ψ, par][{n1_, n2_, m_}, ζ_, η_, ϕ_] := Module[{mm = Abs[m], n = n1 + n2 + Abs[m] + 1}, (Exp[I m ϕ]/Sqrt[π]Sqrt[ζ η]n)Sqrt[((n1 + mm)!/n1!)((n2 + mm)!/n2!)]WhittakerM[n1 + (mm + 1) / 2, mm / 2, (1/n) η]WhittakerM[n2 + (mm + 1) / 2, mm / 2, (1/n) ζ]]Decompose the eigenfunction in terms of spherical eigenfunctions:
Subscript[ψ, sph][{n_, l_, m_}, r_, θ_, ϕ_] := Sqrt[( (n - l - 1)!/(n + l)!)] E^-(r/n) ((2r/n))^l(2/n^2) LaguerreL[n - l - 1, 2l + 1, (2r/n)]SphericalHarmonicY[l, m, θ, ϕ]Parabolic coordinates relate to radial coordinates as
and
:
FunctionExpand[Subscript[ψ, par][{1, 1, -1}, r(1 + Cos[θ]), r(1 - Cos[θ]), ϕ]]//FullSimplify[#, 0 <= θ < π && r ≥ 0]&Simplify[% + (3/Sqrt[15]) Subscript[ψ, sph][{4, 3, -1}, r, θ, ϕ] - (2/Sqrt[10])Subscript[ψ, sph][{4, 1, -1}, r, θ, ϕ]]Green's function of the 3D Coulomb potential:
gf[r1_, r2_, e_] := Module[{ν, a1, a2, δ, 𝓋ℓ = Sqrt[Total[#1 ^ 2]]&}, ν = e ^ (-1 / 2);δ = 𝓋ℓ[r2 - r1];{a1, a2} = 𝓋ℓ[r2] + 𝓋ℓ[r1] + {δ, -δ};(-Gamma[1 - I ν]/4 π δ)Det[(| | |
| ------------------------------------------- | ------------------------------------------- |
| WhittakerW[I ν, (1/2), -(I a2/ν)] | WhittakerM[I ν, (1/2), -(I a1/ν)] |
| WhittakerW^(0, 0, 1)[I ν, (1/2), -(I a2/ν)] | WhittakerM^(0, 0, 1)[I ν, (1/2), -(I a1/ν)] |)]]Plot[gf[{1, 0, 0}, {0, 0, 1}, e], {e, -2, .01}]Properties & Relations (4)
Use FunctionExpand to expand WhittakerM into other functions:
FunctionExpand[WhittakerM[5, 1 / 2, x]]FunctionExpand[WhittakerM[-1 / 4, 1 / 4, x]]FunctionExpand[WhittakerM[0, m, x]]Integrate expressions involving Whittaker functions:
Integrate[Exp[-x]WhittakerM[3, 2, x]WhittakerW[1, -2, x], {x, 0, ∞}]WhittakerM can be represented as a DifferentialRoot:
DifferentialRootReduce[WhittakerM[n, m, x], x]WhittakerM can be represented as a DifferenceRoot:
DifferenceRootReduce[WhittakerM[k, y, z], k]DifferenceRootReduce[WhittakerM[y, k, z], k]Neat Examples (1)
With[{n = 3 / 5, m = 1 / 3, ε = 1*^-12},
Show[Table[ParametricPlot3D[{r Cos[φ], r Sin[φ], Im[(-1)^kExp[2π I m k]WhittakerM[n, m, r Exp[I φ]]]}, {r, 0, 3}, {φ, -π + ε, π - ε}, Mesh -> None, PlotPoints -> 21, PlotStyle -> Directive[Hue[0.85], Opacity[0.6]]], {k, -1, 1}], BoxRatios -> {1, 1, 3}, PlotRange -> {All, All, {-5, 5}}]]Tech Notes
Related Guides
Related Links
History
Text
Wolfram Research (2007), WhittakerM, Wolfram Language function, https://reference.wolfram.com/language/ref/WhittakerM.html.
CMS
Wolfram Language. 2007. "WhittakerM." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/WhittakerM.html.
APA
Wolfram Language. (2007). WhittakerM. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/WhittakerM.html
BibTeX
@misc{reference.wolfram_2026_whittakerm, author="Wolfram Research", title="{WhittakerM}", year="2007", howpublished="\url{https://reference.wolfram.com/language/ref/WhittakerM.html}", note=[Accessed: 13-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_whittakerm, organization={Wolfram Research}, title={WhittakerM}, year={2007}, url={https://reference.wolfram.com/language/ref/WhittakerM.html}, note=[Accessed: 13-June-2026]}