The Series function computes a truncated series expansion of the function call math_function, with respect to the variable z indicated in eqn, about the point z0 (default is z = 0), up to order n (default is n = 6). So, the call Series(math_func, z) is equivalent to Series(math_func, z = 0, 6).

The series returned is computed (outside singularities/branch cuts) using a unified approach, based on the knowledge of the series expansion around the origin (herein called the Sum form).

The Series function is not a replacement for the Maple series command, nor does it have the generality of series. As a function of the MathematicalFunctions package, the goal of Series is to return valuable mathematical information not returned by series (see the examples). As such, Series can be seen as a complement to the series command.

Unlike the Maple series command, the first argument in a call to Series must be a known mathematical function call or a Sum. When the expansion is around z = 0, Series computes the requested series expansion directly from the Sum form of the mathematical function. In this case, the output may consist of a sequence of two objects: first the series, then any restrictions on the function parameters known to the FunctionAdvisor such that the Sum representation and so the series expansion is valid. If the series cannot be computed directly from this sum representation, Series returns FAIL.

If the expansion point is specified as z = z0 with $\mathrm{z0}\ne 0$, Series first attempts to compute a Taylor series expansion. If that fails, it shifts the expansion variable $z\mapsto z+\mathrm{z0}$, then attempts to compute the series around the origin from the Sum form.

In cases that the series is computed from the expansion around the origin, when a call of the form Series( f(g(z)), z ) is processed (where $f$ is a mathematical function and $g$ an arbitrary mathematical expression), the series computation proceeds directly from the Sum form whenever $g\left(0\right)=0$, and otherwise through a call to series with the Sum form.

For some mathematical functions, the series expansion can be computed only by first factoring out abstract powers of the expansion variable, $z$. In this case, the Series command attempts the factorization, and, if successful, returns a structure of the form ${\mathrm{p1}\left(z\right)}^{\mathrm{k1}}\mathrm{series\_1}+{\mathrm{p2}\left(z\right)}^{\mathrm{k2}}\mathrm{series\_2}+...$, where $\mathrm{series\_1}$, $\mathrm{series\_2}$, ... are of type/series and ${\mathrm{p1}\left(z\right)}^{\mathrm{k1}},{\mathrm{p2}\left(z\right)}^{\mathrm{k2}}$, ... are abstract powers of polynomials $\mathrm{p1}$, $\mathrm{p2}$, ... in $z$.

$\mathrm{with}\left(\mathrm{MathematicalFunctions}\,\mathrm{Series}\right)$

$\left[\mathrm{Series}\right]$

$\mathrm{Series}\left(\sin \left(z\right)\,z\right)$

$z-\frac{1}{6}{z}^3+\frac{1}{120}{z}^5+O\left(z^7\right)$

$\mathrm{Series}\left(\arcsin \left(z\right)\,z\right)$

$z+\frac{1}{6}{z}^3+\frac{3}{40}{z}^5+O\left(z^7\right),\wedge \left({\left|z\right|}^2<1\right)$

$\mathrm{series}\left(\arcsin \left(z\right)\&comma;z\right)$

$z+\frac{1}{6}{z}^3+\frac{3}{40}{z}^5+O\left(z^7\right)$

$\mathrm{Series}\left(\mathrm{BesselJ}\left(n\&comma;z\right)\&comma;z\right)$

$z^n\left(\frac{1}{2^n\mathrm{\Gamma }\left(1+n\right)}-\frac{1}{4}\frac{1}{2^n\mathrm{\Gamma }\left(n+2\right)}{z}^2+\frac{1}{32}\frac{1}{2^n\mathrm{\Gamma }\left(3+n\right)}{z}^4+O\left(z^6\right)\right),\wedge \left(n::\left(\neg {\mathbb{Z}}^{-}\right)\right)$

$\mathrm{Series}\left(\mathrm{BesselK}\left(n\&comma;z\right)\&comma;z\&comma;3\right)$

$z^n\left(-\frac{\mathrm{\pi }}{2\mathrm{\Gamma }\left(1+n\right)\sin \left(n\mathrm{\pi }\right)2^n}-\frac{1}{8}\frac{\mathrm{\pi }}{\mathrm{\Gamma }\left(n+2\right)\sin \left(n\mathrm{\pi }\right)2^n}{z}^2+O\left(z^4\right)\right)+z^{-n}\left(\frac{\mathrm{\pi }{2}^n}{2\mathrm{\Gamma }\left(1-n\right)\sin \left(n\mathrm{\pi }\right)}+\frac{1}{8}\frac{\mathrm{\pi }{2}^n}{\mathrm{\Gamma }\left(2-n\right)\sin \left(n\mathrm{\pi }\right)}{z}^2+O\left(z^4\right)\right),\wedge \left(n::\left(\neg \mathbb{Z}\right)\right)$

$\mathrm{Series}\left(\mathrm{KummerU}\left(a\&comma;b\&comma;z\right)\&comma;z\&comma;3\right)$

$\left(\frac{\mathrm{\pi }\csc \left(\mathrm{\pi }b\right)}{\mathrm{\Gamma }\left(b\right)\mathrm{\Gamma }\left(a-b+1\right)}+\frac{\mathrm{\pi }\csc \left(\mathrm{\pi }b\right)a}{\mathrm{\Gamma }\left(b+1\right)\mathrm{\Gamma }\left(a-b+1\right)}z+\frac{1}{2}\frac{\mathrm{\pi }\csc \left(\mathrm{\pi }b\right)a\left(a+1\right)}{\mathrm{\Gamma }\left(b+2\right)\mathrm{\Gamma }\left(a-b+1\right)}{z}^2+O\left(z^3\right)\right)+z^{-b}\left(-\frac{\mathrm{\pi }\csc \left(\mathrm{\pi }b\right)}{\mathrm{\Gamma }\left(2-b\right)\mathrm{\Gamma }\left(a\right)}z-\frac{\mathrm{\pi }\csc \left(\mathrm{\pi }b\right)\left(a-b+1\right)}{\mathrm{\Gamma }\left(3-b\right)\mathrm{\Gamma }\left(a\right)}{z}^2+O\left(z^3\right)\right),a::\left(\neg {\mathbb{Z}}^{\left(0\&comma;-\right)}\right)\wedge \left(b-1\right)::\left(\neg {\mathbb{Z}}^{\left(0\&comma;-\right)}\right)\wedge \left(1-b\right)::\left(\neg {\mathbb{Z}}^{\left(0\&comma;-\right)}\right)$

$\mathrm{Series}\left(\mathrm{KummerU}\left(a\&comma;b\&comma;z\right)\&comma;z\&comma;3\right)\mathbf{assuming}a::\mathrm{nonposint}$

$\mathrm{pochhammer}\left(a-b+1\&comma;-a\right)+\frac{\mathrm{pochhammer}\left(a-b+1\&comma;-a\right)a}{b}z+\frac{1}{2}\frac{\mathrm{pochhammer}\left(a-b+1\&comma;-a\right)\mathrm{pochhammer}\left(a\&comma;2\right)}{\mathrm{pochhammer}\left(b\&comma;2\right)}{z}^2+O\left(z^3\right)$

$\mathrm{convert}\left(\mathrm{InverseJacobiAM}\left(z\&comma;k\right)\&comma;\mathrm{Sum}\right)$

$\sum _{\mathrm{\_m2}=0}^{\mathrm{\infty }}\sum _{\mathrm{\_m1}=0}^{\mathrm{\infty }}\frac{{\left(k^2\right)}^{\mathrm{\_m1}}{\sin \left(z\right)}^{2\mathrm{\_m2}+2\mathrm{\_m1}+1}\mathrm{pochhammer}\left(\frac{1}{2}\&comma;\mathrm{\_m2}\right)\mathrm{pochhammer}\left(\frac{1}{2}\&comma;\mathrm{\_m1}\right)}{\left(2\mathrm{\_m2}+2\mathrm{\_m1}+1\right)\mathrm{\_m2}!\mathrm{\_m1}!}$

$\mathrm{Series}\left(\&comma;k\right)\mathbf{assuming}z::\mathrm{real}$

$\arcsin \left(\sin \left(z\right)\right)+\left(-\frac{\left|\cos \left(z\right)\right|\sin \left(z\right)}{4}+\frac{\arcsin \left(\sin \left(z\right)\right)}{4}\right)k^2+\left(\frac{9\arcsin \left(\sin \left(z\right)\right)}{64}+\frac{3\left|\cos \left(z\right)\right|\sin \left(z\right){\cos \left(z\right)}^2}{32}-\frac{15\left|\cos \left(z\right)\right|\sin \left(z\right)}{64}\right)k^4+O\left(k^6\right)$