SinIntegral[z]
gives the sine integral function
.
SinIntegral
SinIntegral[z]
gives the sine integral function
.
Details
- Mathematical function, suitable for both symbolic and numerical manipulation.
. - SinIntegral[z] is an entire function of
with no branch cut discontinuities. - For certain special arguments, SinIntegral automatically evaluates to exact values.
- SinIntegral can be evaluated to arbitrary numerical precision.
- SinIntegral automatically threads over lists.
- SinIntegral can be used with Interval and CenteredInterval objects. »
Examples
open all close allBasic Examples (6)
SinIntegral[2.8]Plot[SinIntegral[x], {x, -20, 20}]Plot over a subset of the complexes:
ComplexPlot3D[SinIntegral[z], {z, -2 - 2I, 2 + 2I}, PlotLegends -> Automatic]SinIntegral'[x]Series expansion at the origin:
Series[SinIntegral[x], {x, 0, 10}]Asymptotic expansion at Infinity:
Series[SinIntegral[z], {z, ∞, 3}]//NormalScope (37)
Numerical Evaluation (5)
Evaluate numerically to high precision:
N[SinIntegral[2], 50]The precision of the output tracks the precision of the input:
SinIntegral[2.0000000000000000000000]Evaluate for complex arguments:
SinIntegral[2.5 + I]Evaluate SinIntegral efficiently at high precision:
SinIntegral[2`500]//TimingSinIntegral[2`10000];//TimingCompute worst-case guaranteed intervals using Interval and CenteredInterval objects:
SinIntegral[Interval[{1.1, 1.2}]]SinIntegral[CenteredInterval[2, 1 / 100]]Or compute average-case statistical intervals using Around:
SinIntegral[Around[2, 0.01]]Compute the elementwise values of an array:
SinIntegral[{{-.2, 0.2}, {0, -.1}}]Or compute the matrix SinIntegral function using MatrixFunction:
MatrixFunction[SinIntegral, {{-.2, 0.2}, {0, -.1}}]Specific Values (3)
SinIntegral[0]SinIntegral[{-Infinity, Infinity, -I Infinity, I Infinity}]SinIntegral[ComplexInfinity]Find a local maximum as a root of
:
xmax = Solve[D[SinIntegral[x], x] == 0 && 0 < x < 4, x][[1, 1]]Plot[SinIntegral[x], {x, -3π, 3π}, Epilog -> Style[Point[{xmax[[2]], SinIntegral[xmax[[2]]]}], PointSize[Large], Red]]Visualization (2)
Plot the SinIntegral function:
Plot[SinIntegral[x], {x, -4π, 4π}]//RasterizeComplexContourPlot[Re[SinIntegral[z]], {z, -2π - π I, 2π + π I}, Contours -> 24]ComplexContourPlot[Im[SinIntegral[z]], {z, -2π - π I, 2π + π I}, Contours -> 24]Function Properties (10)
SinIntegral is defined for all real and complex values:
FunctionDomain[SinIntegral[x], x]FunctionDomain[SinIntegral[z], z, Complexes]Approximate function range of SinIntegral:
FunctionRange[SinIntegral[x], x, y]//NSinIntegral is an odd function:
SinIntegral[-x]SinIntegral is an analytic function of x:
FunctionAnalytic[SinIntegral[x], x]SinIntegral is neither non-decreasing nor non-increasing:
FunctionMonotonicity[SinIntegral[x], x]SinIntegral is not injective:
FunctionInjective[SinIntegral[x], x]Plot[{SinIntegral[x], 1.6}, {x, -15, 15}]SinIntegral is not surjective:
FunctionSurjective[SinIntegral[x], x]Plot[{SinIntegral[x], 3}, {x, -15, 15}]SinIntegral is neither non-negative nor non-positive:
FunctionSign[SinIntegral[x], x]SinIntegral has no singularities or discontinuities:
FunctionSingularities[SinIntegral[x], x]FunctionDiscontinuities[SinIntegral[x], x]SinIntegral is neither convex nor concave:
FunctionConvexity[SinIntegral[x], x]Differentiation (3)
D[SinIntegral[x], x]Table[D[SinIntegral[x], {x, n}], {n, 1, 4}]Plot[Evaluate[%], {x, -2π, 2π}, PlotLegends -> {"First Derivative", "Second Derivative", "Third Derivative", "Fourth Derivative"}]//RasterizeD[SinIntegral[x], {x, n}]Integration (3)
Indefinite integral of SinIntegral:
Integrate[SinIntegral[x], x]Definite integral of an odd integrand over an interval centered at the origin is 0:
Integrate[SinIntegral[x], {x, -2π, 2 π}]Integrate[x SinIntegral[Sqrt[x]], x]Integrate[SinIntegral[x]Sin[x], x]Integrate[(SinIntegral[x] - π / 2)^2, {x, 0, Infinity}]Series Expansions (4)
Taylor expansion for SinIntegral:
Series[SinIntegral[x], {x, 0, 7}]Plot the first three approximations for SinIntegral around
:
terms = Normal@Table[Series[SinIntegral[x], {x, 0, m}], {m, 1, 5, 2}];
Plot[{SinIntegral[x], terms}, {x, -1.5 π, 1.5π}]//RasterizeGeneral term in the series expansion of SinIntegral:
SeriesCoefficient[SinIntegral[x], {x, 0, n}]Find series expansion at infinity:
Series[SinIntegral[x], {x, Infinity, 3}]Give the result for an arbitrary symbolic direction
:
Series[SinIntegral[x], {x, DirectedInfinity[z], 1}]SinIntegral can be applied to power series:
SinIntegral[x + (x^2/2) + (x^3/9) + O[x]^4]Function Identities and Simplifications (3)
Use FullSimplify to simplify expressions containing sine integrals:
FullSimplify[2SinIntegral[z] + I ExpIntegralEi[I z] - I ExpIntegralEi[-I z], z > 0]Simplify expressions to SinIntegral:
(I/2)(Gamma[0, -I z] - Gamma[0, I z] + Log[-I z] - Log[I z])//FullSimplifySinIntegral[I z]SinIntegral[Sqrt[z^2]]//FunctionExpandFunction Representations (4)
Series representation of SinIntegral:
z Sum[((-1)^k z^2k/(1 + 2k)^2 (2k)!), {k, 0, Infinity}]//FullSimplifySinIntegral can be represented in terms of MeijerG:
MeijerGReduce[SinIntegral[x], x]Activate[%]SinIntegral can be represented as a DifferentialRoot:
DifferentialRootReduce[SinIntegral[x], x]TraditionalForm formatting:
SinIntegral[x]//TraditionalFormGeneralizations & Extensions (1)
Applications (6)
Plot the absolute value in the complex plane:
Plot3D[Abs[SinIntegral[x + I y]], {x, -20, 20}, {y, -3, 3}]Real part of the Euler–Heisenberg effective action:
action[k_, {m_, e_}, {a_, b_}] := With[{x = k π (m^2/e a), y = k π (m^2/e b)}, -(e^2a b/4 π^3k)(Coth[k π (b/a)](CosIntegral[x] Cos[x] + (SinIntegral[x] - π / 2)Sin[x]) -
Coth[k π (a/b)](ExpIntegralEi[y] Exp[-y] + ExpIntegralEi[-y] Exp[y]) / 2)]Assuming[{q > 0 && a > 0 && b > 0 && m > 0 && k > 0},
Series[action[k, {m, q}, {a, b}], {q, 0, 4}]//Normal//FullSimplify]Gibbs phenomenon for a square wave:
Plot[Evaluate[Table[1 / 2 - 4x Sum[Sinc[2Pi (2k - 1)x], {k, n}], {n, 20}]], {x, 0, 1}]Show[%, PlotRange -> {{0.45, 0.55}, {0.8, 1.2}}]Compute the asymptotic overshoot:
N[1 / 2 + SinIntegral[Pi] / Pi]Solve a differential equation:
DSolve[x (y[x] + Derivative[2][y][x]) == Cos[x], y[x], x]Integrate a composition of trigonometric functions:
Integrate[Sin[Tan[x]], x]nielsen[t_] := {CosIntegral[t], SinIntegral[t]}ParametricPlot[nielsen[t], {t, 1 / 25, 25}, PlotRange -> All]The curvature is a simple function of the parameter:
FullSimplify[FrenetSerretSystem[nielsen[t], t][[1, 1]], t∈Reals]Properties & Relations (7)
Parity transformation is automatically applied:
SinIntegral[-z]SinIntegral[I z]Use FullSimplify to simplify expressions containing sine integrals:
FullSimplify[2SinIntegral[z] + I ExpIntegralEi[I z] - I ExpIntegralEi[-I z], z > 0]FindRoot[SinIntegral[z] == 1 / Pi, {z, 1}]Obtain SinIntegral from integrals and sums:
Integrate[(Sin[t]/t), {t, 0, z}]Underoverscript[∑, k = 0, ∞]((-1)^k z^2 k/(1 + 2 k)^2 (2 k)!)Obtain SinIntegral from a differential equation:
w[z] /. DSolve[z w'''[z] + 2 w''[z] + z w'[z] == 0, w[z], z]//FunctionExpandDet[Outer[D[#1, {x, #2}]&, {1, SinIntegral[x], CosIntegral[x]}, {0, 1, 2}]]//FullSimplifyCompare with Wronskian:
Wronskian[{1, SinIntegral[x], CosIntegral[x]}, x]Integrate[z SinIntegral[Sqrt[z]], z]LaplaceTransform[SinIntegral[t], t, s]Possible Issues (2)
SinIntegral can take large values for moderate‐size arguments:
SinIntegral[10000. (1 + I)]A larger setting for $MaxExtraPrecision can be needed:
N[SinIntegral[10 ^ 60] - Pi / 2, 20]Block[{$MaxExtraPrecision = 70}, N[SinIntegral[10 ^ 60] - Pi / 2, 20]]Tech Notes
Related Guides
Related Links
History
Introduced in 1991 (2.0) | Updated in 2021 (13.0) ▪ 2022 (13.1)
Text
Wolfram Research (1991), SinIntegral, Wolfram Language function, https://reference.wolfram.com/language/ref/SinIntegral.html (updated 2022).
CMS
Wolfram Language. 1991. "SinIntegral." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2022. https://reference.wolfram.com/language/ref/SinIntegral.html.
APA
Wolfram Language. (1991). SinIntegral. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/SinIntegral.html
BibTeX
@misc{reference.wolfram_2026_sinintegral, author="Wolfram Research", title="{SinIntegral}", year="2022", howpublished="\url{https://reference.wolfram.com/language/ref/SinIntegral.html}", note=[Accessed: 13-June-2026]}
BibLaTeX
@online{reference.wolfram_2026_sinintegral, organization={Wolfram Research}, title={SinIntegral}, year={2022}, url={https://reference.wolfram.com/language/ref/SinIntegral.html}, note=[Accessed: 13-June-2026]}