Mathematical Functions
Interesting and relevant developments in the MathematicalFunctions and FunctionAdvisor projects happened for Maple 2016, regarding both the user-interface and the mathematics, making this release a more complete and user-friendly environment to work with mathematical functions.
Gaps were filled regarding mathematical formulas, with more identities for all of BesselI, BesselK, BesselY, ChebyshevT, ChebyshevU, Chi, Ci, FresnelC, FresnelS, GAMMA(z), HankelH1, HankelH2, InverseJacobiAM, the twelve InverseJacobiPQ for P, Q in [C,D,N,S], KelvinBei, KelvinBer, KelvinKei, KelvinKer, LerchPhi, arcsin, arcsinh, arctan, ln;
Developments happened in the Mathematical function package, to both compute with symbolic sequences and symbolic nth order derivatives of algebraic expressions and functions;
The input FunctionAdvisor⁡differentiate_rule,mathematical_function now returns both the first derivative (old behavior) and the nth symbolic derivative (new behavior) of a mathematical function;
A new topic, plot, used as FunctionAdvisor⁡plot,mathematical_function, now returns 2-D and 3-D plots for each mathematical function, following the NIST Digital Library of Mathematical Functions;
The FunctionAdvisor⁡display,mathematical_function was redesigned, so that the display keyword is not necessary anymore. The command now displays more information about any mathematical function, and organized into a Section with subsections for each of the different topics, making it simpler to find the information one needs without getting distracted by a myriad of formulas that are not related to what one is looking for.
To display special functions and sequences using textbook notation as shown in this page, use extended typesetting:
interfacetypesetting=extended: Typesetting:-EnableTypesetRuleTypesetting:-SpecialFunctionRules:
More mathematics
More powerful symbolic differentiation (nth order derivative)
Mathematical handling of symbolic sequences
Visualization of mathematical functions
Section and subsections displaying properties of mathematical functions
More mathematical knowledge is in place, more identities, differentiation rules of special functions with respect to their parameters, differentiation of functions whose arguments involve symbolic sequences with an indeterminate number of operands, and sum representations for special functions under different conditions on the functions' parameters.
Examples
More identities for all of BesselI, BesselK, BesselY, ChebyshevT, ChebyshevU, Chi, Ci, FresnelC, FresnelS, GAMMA(z), HankelH1, HankelH2, InverseJacobiAM, the twelve InverseJacobiPQ for P, Q in [C,D,N,S], KelvinBei, KelvinBer, KelvinKei, KelvinKer, LerchPhi, arcsin, arcsinh, arctan, ln
FunctionAdvisoridentities, ln
ln⁡z=I⁢arg⁡z+ln⁡z,ln⁡ⅇz=z,z::real,ln⁡y+z=ln⁡y+ln⁡12⁢y+z+ln⁡2⁢y+z+2⁢arctanh⁡z2⁢y+z,0<y,ln⁡y⁢z=ln⁡y+ln⁡z,0≤y+z,ln⁡zy=y⁢ln⁡z+2⁢I⁢π⁢π−ℑ⁡y⁢ln⁡z2⁢π,ln⁡za⁢yb=a⁢ln⁡z+b⁢ln⁡y+2⁢I⁢π⁢π−ℑ⁡a⁢ln⁡z−ℑ⁡b⁢ln⁡y2⁢π
FunctionAdvisoridentities, BesselK
Ka⁡I⁢z=−π⁢Ya⁡z2⁢Ia+Ja⁡z⁢ln⁡z−ln⁡I⁢zIa,a::ℤ,Ka⁡I⁢z=−π⁢za⁢Ya⁡z2⁢I⁢za+π⁢Ja⁡z⁢−I⁢zaza+za⁢cos⁡π⁢aI⁢za⁢csc⁡π⁢a2,a::¬ℤ,Ka⁡−z=−1a⁢Ka⁡z+Ia⁡z⁢ln⁡z−ln⁡−z,a::ℤ,Ka⁡−z=za⁢Ka⁡z−za+π⁢za−za−−zaza⁢Ia⁡z⁢csc⁡π⁢a2,a::¬ℤ,Ka⁡b⁢c⁢zqp=b⁢cp⁢zp⁢qa⁢Ka⁡b⁢cp⁢zp⁢qb⁢c⁢zqpa−π⁢csc⁡π⁢a⁢Ia⁡b⁢cp⁢zp⁢q⁢b⁢c⁢zqpab⁢cp⁢zp⁢qa−b⁢cp⁢zp⁢qab⁢c⁢zqpa2,a::¬ℤ∧2⁢p::ℤ,Ka⁡b⁢c⁢zqp=c⁢zqpcp⁢zp⁢qa⁢Ka⁡b⁢cp⁢zp⁢q−−1a⁢Ia⁡b⁢cp⁢zp⁢q⁢ln⁡b⁢c⁢zqp−ln⁡b⁢cp⁢zp⁢q,a::ℤ∧2⁢p::ℤ,Ka⁡z=2⁢a−1⁢Ka−1⁡zz+Ka−2⁡z,Ka⁡z=−2⁢a+1⁢Ka+1⁡zz+Ka+2⁡z
FunctionAdvisoridentities, Ci
Ci⁡−z=Ci⁡z+ln⁡−z−ln⁡z,Ci⁡I⁢z=Chi⁡z−ln⁡z+ln⁡I⁢z,Ci⁡z=−Ei1⁡I⁢z2−Ei1⁡−I⁢z2+I⁢csgn⁡z−1⁢csgn⁡I⁢z⁢π2
FunctionAdvisoridentities,InverseJacobiSN
sn⁡sn−1⁡z|k|k=z,sn−1⁡z|k=I⁢cs−1⁡Iz|−k2+1,sn−1⁡z|k=am−1⁡arcsin⁡z|k,sn−1⁡z|k=K⁡k−cd−1⁡z|k,sn−1⁡z|k=K⁡k−dc−1⁡1z|k,sn−1⁡z|k=I⁢sc−1⁡−I⁢z|−k2+1,sn−1⁡z|k=F⁡z,k,z<1∧k2<1
More differentiation rules of special functions with respect to their parameters
Equating the inert derivative (on hold) to the active derivative (computed)
%diff = diffLaguerreLa, b, z, a
ⅆⅆaLab⁡z=−Γ⁡1+b+a⁢∑k=0∞⁡zk⁢−ak⁢Ψ⁡k−ak!⁢Γ⁡1+k+ba!+π⁢cot⁡π⁢a+Ψ⁡1+b+a⁢Lab⁡z
%diff = diffLaguerreLa, b, z, b
ⅆⅆbLab⁡z=Ψ⁡1+b+a⁢Lab⁡z−Γ⁡1+b+a⁢∑k=0∞⁡zk⁢−ak⁢Ψ⁡1+k+bk!⁢Γ⁡1+k+bΓ⁡1+a
Differentiation rules of the hypergeometric pFq and MeijerG functions for an indeterminate (symbolic sequence) number of parameters:
%diff = diffhypergeomai$i=1..p, bi$i=1..q, z,z
ⅆⅆzFqp⁡a1,...,ap;b1,...,bq;z=∏i=1p⁡ai⁢Fqp⁡a1+1,...,ap+1;b1+1,...,bq+1;z∏i=1q⁡bi
The system can now also compute the nth symbolic order derivative of these hypergeometric functions of an indeterminate number of parameters:
%diff = diffhypergeomai$i=1..p, bi$i=1..q, z,z$n
ⅆnⅆznFqp⁡a1,...,ap;b1,...,bq;z=∏i=1p⁡ain⁢Fqp⁡n+a1,...,n+ap;n+b1,...,n+bq;z∏i=1q⁡bin
In this development converge a number of lower level developments a) the Maple system now operates mathematically with symbolic sequences, addition, multiplication and differentiation, b) there is new typesetting for displaying of symbolic sequences, c) there is more mathematical knowledge in the differentiation rules, taking advantage of a) and b).
The kth order derivative of the more general MeijerG function with an indeterminate number (symbolic sequence) of parameters:
%diff = diffMeijerGai $ i=1..n,bi $ i=n+1..p, bi $ i=1..m,bi $ i=m+1..q, z ,z$k
ⅆkⅆzkGp,qm,n⁡z|a1,...,an,bn+1,...,bpb1,...,bm,bm+1,...,bq=G1+p,q+1m,n+1⁡z|−k,a1−k,...,an−k,bn+1−k,...,bp−kb1−k,...,bm−k,0,bm+1−k,...,bq−k
The first order derivative of this function:
eval, k=1
ⅆⅆzGp,qm,n⁡z|a1,...,an,bn+1,...,bpb1,...,bm,bm+1,...,bq=G1+p,q+1m,n+1⁡z|−1,a1−1,...,an−1,bn+1−1,...,bp−1b1−1,...,bm−1,0,bm+1−1,...,bq−1
These formulas involve a rather high level of abstraction and required a number of underlying supporting routines to do all the mathematics correctly.
More sum representations of mathematical functions under different conditions on their parameters
FunctionAdvisorsum,polylog
* Partial match of "sum" against topic "sum_form".
Lia⁡z=∑_k1=1∞⁡z_k1_k1a,z<1,Lia⁡z=z2+∑_k1=−∞∞⁡Γ⁡1−a,2⁢I⁢_k1⁢π−ln⁡z⁢2⁢I⁢_k1⁢π−ln⁡z−1+a,a::ℤ+,Lia⁡z=z2+∑_k1=−∞∞⁡∑_k2=0∞⁡−2⁢I⁢_k1⁢π+ln⁡z_k22⁢I⁢_k1⁢π−ln⁡z−1+a⁢Γ⁡1+_k2⁢−_k2−1+a+Γ⁡1−a⁢2⁢I⁢_k1⁢π−ln⁡z−1+a,a::¬ℤ+
FunctionAdvisorsum,Zeta
ζ⁡s=∑_k1=1∞⁡1_k1s,1<ℜ⁡s,ζ⁡s=π−12+s⁢−12−s2!⁢∑_k1=1∞⁡_k1−1+ss2−1!,ℜ⁡s<0,ζn⁡s=ⅆnⅆsn∑_k1=1∞⁡1_k1s,n::ℤ0,+∧1<ℜ⁡s,ζn⁡s=∂n∂snπ−12+s⁢−12−s2!⁢∑_k1=1∞⁡_k1−1+ss2−1!,n::ℤ0,+∧ℜ⁡s<0,ζn⁡s,a=ⅆnⅆsn∑_k1=0∞⁡1a+_k1s,n::ℤ0,+∧1<ℜ⁡s,ζn⁡s,a=∂n∂sn∑_k1=0∞⁡−1_k1⁢s_k1⁢∑_k2=1∞⁡1_k2_k1+s1<ℜ⁡_k1+sπ−12+_k1+s⁢−12−_k12−s2!⁢∑_k2=1∞⁡_k2−1+_k1+s_k12+s2−1!ℜ⁡_k1+s<0⁢a_k1_k1!+1as,n::ℤ0,+∧ℜ⁡s<0
Significant developments happened in the computation of the nth order derivative of mathematical functions and algebraic expressions involving them.
%diff=difffα⁢z+β,z$n
∂n∂znfα⁢z+β=αn⁢fα⁢z+β⁢ln⁡fn
The symbolic differentiation of binomial(z, m)
%diff=diffbinomialz,n,z$n
ⅆnⅆznzn=∑_k1=1n⁡−1_k1+n⁢Sn_k1⁢_k1−n+1n⁢z−n+1_k1−nn!
And for the first time in computer algebra systems, we now have the Faà di Bruno formula for the nth derivative of a composite function working, using the IncompleteBell polynomials and taking advantage of the new developments in the mathematical handling and display of symbolic sequences:
%diff=difffg⁡z,z$n
ⅆnⅆznf⁡g⁡z=∑_k2=0n⁡D_k2⁡f⁡g⁡z⁢IncompleteBellB⁡n,_k2,ⅆⅆzg⁡z,...,ⅆn−_k2+1ⅆzn−_k2+1g⁡z
All these results can also be verified with ease, for instance, the third derivative of the composite function f⁡g⁡z is given by
eval,n=3
ⅆ3ⅆz3f⁡g⁡z=∑_k2=03⁡D_k2⁡f⁡g⁡z⁢IncompleteBellB⁡3,_k2,ⅆⅆzg⁡z,...,ⅆ4−_k2ⅆz4−_k2g⁡z
value
D3⁡f⁡g⁡z⁢ⅆⅆzg⁡z3+3⁢D2⁡f⁡g⁡z⁢ⅆⅆzg⁡z⁢ⅆ2ⅆz2g⁡z+D⁡f⁡g⁡z⁢ⅆ3ⅆz3g⁡z=D3⁡f⁡g⁡z⁢ⅆⅆzg⁡z3+3⁢D2⁡f⁡g⁡z⁢ⅆⅆzg⁡z⁢ⅆ2ⅆz2g⁡z+D⁡f⁡g⁡z⁢ⅆ3ⅆz3g⁡z
Check that the left-hand side in (2.1.5) is actually equal to the right-hand side
evalb
true
These developments regarding nth order symbolic differentiation are now displayed by the FunctionAdvisor when the differentiation rule of a mathematical function is requested
FunctionAdvisordiff, ln
* Partial match of "diff" against topic "differentiation_rule".
ⅆⅆzln⁡z=1z,ⅆnⅆznln⁡z=ln⁡zn=0−1n−1⁢n−1!znotherwise
Symbolic sequences enter various formulations in mathematics. Their computerized mathematical handling, however, was never implemented - only a representation for them existed in the Maple system. In connection with this, a new subpackage, Sequences, within the MathematicalFunctions package, has been developed.
The most typical cases of symbolic sequences are:
1) A sequence of numbers - say from n to m - frequently displayed as
n,...,m
2) A sequence of one object, say a, repeated say p times, frequently displayed as
a,...,a⏟p times
3) A more general sequence, as in 1), but of different objects and not necessarily numbers, frequently displayed as
an,...,am
or likewise a sequence of functions
f⁡n,...,f⁡m
In all these cases, of course, none of n, m, or p are known: they are just symbols, or algebraic expressions, representing integer values.
A representation for these typical cases of symbolic sequences has been present in Maple for a long time using the `$` operator. Cases 1), 2) and 3) above are respectively entered as `$`n .. m, `$`a,p, and `$`ai, i = n .. m or `$`fi, i = n .. m. The typesetting of these symbolic sequences, however, did not exist. More relevant: too little could be done with these objects; the rest of Maple did not know how to add, multiply, differentiate, or map an operation over the elements of a symbolic sequence, nor for instance count the sequence's number of elements.
All these operations on symbolic sequences are now implemented and functional.
First of all, now these three types of sequences have textbook-like typesetting:
`$`n .. m
a$p
ai $ i = n .. m
Moreover, this now permits textbook display of mathematical functions that depend on sequences of parameters, for example:
hypergeomai $ i = 1 .. p,bi $ i = 1 .. q,z
Fqp⁡a1,...,ap;b1,...,bq;z
More interestingly, these new developments now permit differentiating these functions even when their arguments are symbolic sequences, and displaying the result as in textbooks, with copy and paste working properly, for instance
%diff = diff,z
This enhances the representation capabilities in different relevant ways; to mention but one, this made possible the implementation of the Faà di Bruno formula for the nth symbolic derivative of composite functions for the first time in computer algebra systems.
To access the mathematics of symbolic sequences in Maple, first load the corresponding package:
withMathematicalFunctions:-Sequences
Add,Differentiate,Map,Multiply,Nops
With these commands, it is now possible to add, multiply, differentiate, or map an operation over the elements of a symbolic sequence, as well as count the sequence's number of elements.
For example, here are the three types of symbolic sequences mentioned, with textbook-like typesetting, and on which the operations Add, Multiply, Differentiate, Map, and Nops can now be performed:
`$`a,p
`$`ai, i = n .. m
Nops
m−n+1
Add
m−n+1⁢n+m2
Multiply
m!n−1!
MapInt,,x
∫nⅆx,...,∫mⅆx
p
a⁢p
ap
Mapf,
f⁡a,...,f⁡a⏟p times
Differentiate,a
ⅆⅆaf⁡a,...,ⅆⅆaf⁡a⏟p times
∑i=nm⁡ai
∏i=nm⁡ai
Differentiate,ak
1k=n0otherwise,...,1k=m0otherwise
When working with mathematical functions, it is frequently desired to have a rapid glimpse of the shape of the function for some sample values of their parameters. Following the NIST Digital Library of Mathematical Functions, a new option, plot, has now been implemented.
The Jacobi elliptic sn and Weierstrass P functions,
FunctionAdvisorplot, JacobiSN
In the first of the 3-D plots, for real values of the parameters n and z, sn⁡z|1−ⅇ−n is real and its value is on the vertical axis, while in last three 3-D complex plots, the coloring of the surface follows the value of the argument θ in z=z⁢ⅇI⁢θ, while on the vertical axis the absolute value of the function is plotted.
For the Weierstrass P function,
FunctionAdvisorplot, WeierstrassP
Each of these plots can be rotated, selected with the mouse, or copied and pasted elsewhere in the worksheet for further analysis.
Until recently, the display of a whole set of mathematical information regarding a function was somehow cumbersome, appearing all together on the screen. That display was and is still available via entering, for instance for the sin function, FunctionAdvisor⁡table, sin . That returns a table of information that can be used programmatically.
With time however, the FunctionAdvisor evolved into a consultation tool, where a better organization of the information being displayed is required, making it simpler to find the information we need without being distracted by a screen full of complicated formulas.
To address this requirement, the FunctionAdvisor now returns the information organized into a Section with subsections, built using the DocumentTools package. This enhances the presentation significantly.
For example, for the Ei and GAMMA functions
FunctionAdvisorEi
Ei
describe
Ei=exponential integral
definition
Ei⁡z=PV⁢∫−∞zⅇ_k1_k1ⅆ_k1
Andz::real
Eia⁡z=∫1∞1ⅇ_k1⁢z⁢_k1aⅆ_k1
And0<ℜ⁡z
analytic extension
Ei⁡z=γ−ln⁡1z2+ln⁡z2+z⁢F22⁡1,1;2,2;z
classify function
Ei_related
1F1
periodicity
Ei⁡z
No periodicity
Eia⁡z
plot
singularities
z=∞+∞⁢I
a=∞+∞⁢I
branch points
z∈0,∞+∞⁢I
branch cuts
z<0
special values
Ei⁡∞⁢I=I⁢π
Ei⁡−∞⁢I=−I⁢π
Ei⁡0=−∞
Ei⁡∞=∞
Ei⁡−∞=0
Ei0⁡z=ⅇ−zz
Ei−1⁡z=ⅇ−z⁢1+zz2
Ei0⁡0=∞+∞⁢I
Ei1⁡0=∞+undefined⁢I
Eia⁡0=1a−1
And1<ℜ⁡a
Ei⁡z=−U⁡1,1,−z⁢ⅇz
Andz<0
Ei⁡z=Li⁡ⅇz
−π≤ℑ⁡z≤π
Ei12⁡z=1−erf⁡z⁢πz
Andz≠0
Ei−12⁡z=−π⁢erf⁡z2⁢z32+π2⁢z32+1z⁢ⅇz
Eia⁡z=−Li⁡1ⅇz−ln⁡−1z2+ln⁡−z2−ln⁡z
a=1And−π≤ℑ⁡zAndℑ⁡z≤π
Eia⁡0=∞+∞⁢I
Andℜ⁡a<1
identities
Ei⁡z2=Ei⁡z+z2z−1⁢Shi⁡z+ln⁡1z2+ln⁡−I⁢z2+ln⁡I⁢z2−ln⁡z2
Ei1⁡z=−Ei⁡−z+ln⁡−z2−ln⁡−1z2−ln⁡z
Ei−a⁡z=a!⁢ⅇ−z⁢∑_k1=0a⁡z_k1−a−1_k1!
Anda::nonnegint
Eia⁡z=z⁢Ei−2+a⁡z+−2+a−z⁢Ei−1+a⁡z−1+a
Eia⁡z=−a+z⁢Ei1+a⁡z+1+a⁢Ei2+a⁡zz
sum form
Ei⁡z=γ−ln⁡1z2+ln⁡z2+∑_k1=1∞⁡z_k1_k1⁢_k1!
with no restrictions on ⁡z
Eia⁡z=∑_k1=0∞⁡−z_k1Γ⁡_k1+1⁢−1+a−_k1+za−1⁢Γ⁡1−a
Anda::Not⁡posint
Eia⁡z=∑_k1=0∞⁡−1a+2⁢_k1⁢−Ψ⁡_k1+1+ln⁡z⁢za−1+_k1Γ⁡_k1+1⁢Γ⁡a⁢ⅇz+∑_k1=0−2+a⁡−1_k1⁢Γ⁡−1+a−_k1⁢z_k1Γ⁡a⁢ⅇz
a::posintAndz≠0
Eia⁡z=∑_k1=0−a⁡−a!⁢za−1+_k1ⅇz⁢_k1!
a::negintAndz≠0
series
series⁡Ei⁡z,z,4=γ+ln⁡z+z+14⁢z2+118⁢z3+O⁡z4
series⁡Eia⁡z,a,4=ⅇ−zz−G1,22,0⁡z|0−1,−1⁢a+G2,33,0⁡z|0,0−1,−1,−1⁢a2−G3,44,0⁡z|0,0,0−1,−1,−1,−1⁢a3+O⁡a4
asymptotic expansion
asympt⁡Ei⁡z,z,4=1z+1z2+2z3+O⁡1z4⁢ⅇz
asympt⁡Eia⁡z,z,4=1z−a!a−1!⁢z2+a+1!a−1!⁢z3+O⁡1z4ⅇz
integral form
Ei⁡z=∫0zⅇ_k1−1_k1ⅆ_k1+ln⁡z+γ
differentiation rule
∂∂aEia⁡z=−za⁢G2,33,0⁡z|0,0−1,−1,−a
∂∂zEia⁡z=−Eia−1⁡z
∂n∂znEia⁡z=−−1n⁢G2,31,2⁡z|0,a0,a−1,n−zn+π⁢za−1−nΓ⁡a−n⁢sin⁡π⁢a
DE
f⁡z=Ei⁡z
ⅆ2ⅆz2f⁡z=ⅆⅆzf⁡z⁢z−1z
f⁡z=Eia⁡z
ⅆ2ⅆz2f⁡z=a−z−2⁢ⅆⅆzf⁡zz+a−1⁢f⁡zz
FunctionAdvisorGAMMA
GAMMA
Γ=Gamma and incomplete Gamma functions
Γ⁡z=∫0∞_k1z−1ⅇ_k1ⅆ_k1
Γ⁡a,z=Γ⁡a−za⁢∫01_t1a−1ⅇ_t1⁢zⅆ_t1
And0<ℜ⁡a
Γ⁡z=πsin⁡π⁢z⁢Γ⁡1−z
Andℜ⁡z<0
GAMMA_related
Γ⁡z
Γ⁡a,z
z::nonposint
No branch points
a::Not⁡posintAndz∈0,∞+∞⁢I
No branch cuts
a::Not⁡posintAndz<0
Γ⁡−12=−2⁢π
Γ⁡12=π
Γ⁡1=1
Γ⁡−1=∞+∞⁢I
Γ⁡0=∞+∞⁢I
Γ⁡∞=∞
Γ⁡−∞=undefined
Γ⁡∞⁢I=0
Γ⁡−∞⁢I=0
Γ⁡−1,z=Ei2⁡zz
Γ⁡−12,z=−2⁢π⁢erfc⁡z+2⁢ⅇ−zz
Γ⁡0,z=Ei1⁡z
Γ⁡12,z=π⁢erfc⁡z
Γ⁡1,z=ⅇ−z
Γ⁡a,0=Γ⁡a
Γ⁡a,∞=0
Γ⁡a,−∞=∞+∞⁢I
Γ⁡−z=−π⁢csc⁡π⁢zΓ⁡z+1
Γ⁡z+y=zy⁢Γ⁡z
Γ⁡z=22⁢z⁢∏_k1=1∞⁡Γ⁡12+z2_k1πz
Γ⁡z=234−z+cos⁡2⁢π⁢z4⁢πsin⁡π⁢z22⁢2⁢z−2!!
Γ⁡a+n,z=an⁢Γ⁡a,z+za+n−1⁢ⅇ−z⁢∑_k1=0n−1⁡1−a−n_k1⁢−z−_k1
Andn::nonnegint
Γ⁡a−n,z=−1n⁢Γ⁡a,z1−an−ⅇ−z⁢za−n⁢∑_k1=0n−1⁡z_k1a−n_k1+1
Γ⁡a,z=∑_k1=0∞⁡−za⁢−z_k1Γ⁡_k1+1⁢_k1+a+Γ⁡a
Anda::Not⁡nonposint
Γ⁡a,z=∑_k1=0−a−1⁡Γ⁡−a−_k1⁢z_k1+a⁢−1_k1Γ⁡1−a⁢ⅇz+∑_k1=0∞⁡Ψ⁡_k1+1−ln⁡z⁢z_k1⁢−1−a+2⁢_k1ⅇz⁢Γ⁡1−a⁢Γ⁡_k1+1
a::nonposintAndz≠0
Γ⁡a,z=∑_k1=0a−1⁡a−1!⁢z_k1ⅇz⁢_k1!
series⁡Γ⁡z,z,4=z−1−γ+π212+γ22⁢z+−ζ⁡33−π2⁢γ12−γ36⁢z2+π4160+ζ⁡3⁢γ3+π2⁢γ224+γ424⁢z3+O⁡z4
series⁡Γ⁡a,z,a,4=Ei1⁡z+Ei1⁡z⁢ln⁡z+G2,33,0⁡z|1,10,0,0⁢a+Ei1⁡z⁢ln⁡z22+ln⁡z⁢G2,33,0⁡z|1,10,0,0+G3,44,0⁡z|1,1,10,0,0,0⁢a2+Ei1⁡z⁢ln⁡z36+ln⁡z2⁢G2,33,0⁡z|1,10,0,02+ln⁡z⁢G3,44,0⁡z|1,1,10,0,0,0+G4,55,0⁡z|1,1,1,10,0,0,0,0⁢a3+O⁡a4
asympt⁡Γ⁡z,z,4=2⁢π⁢1z+2⁢π⁢1z3212+2⁢π⁢1z52288−139⁢2⁢π⁢1z7251840+O⁡1