皮尔西积分

皮尔西积分(Pearcey Integral)是一种在论述光的传播、光的衍射、分岔理論、突变理论以及关于特殊函数的渐进展开式的研究中常见的多鞍点积分，其定义为

P(x,y)=\int_{-\infty}^{\infty}\exp(I(t^4+xt^2+y*t))

={-(1/2)\sqrt(Pi)exp(-(1/4I)y^2/x)(limit(exp(It^4)erf(Ixt/\sqrt(-Ix)+(1/2I)y/\sqrt(-Ix)), t = infinity))/\sqrt(-Ix)+(1/2)\sqrt(Pi)exp(-(1/4I)y^2/x)(limit(exp(It^4)erf(Ixt/\sqrt(-Ix)+(1/2I)y/\sqrt(-Ix)), t = -infinity))/\sqrt(-Ix)+(2I)\sqrt(Pi)exp(-(1/4I)y^2/x)(int(exp(It^4)erf(Ixt/\sqrt(-Ix)+(1/2I)y/\sqrt(-Ix))t^3, t = -infinity .. infinity))/\sqrt(-Ix)}

被积分函数的鞍点
对于大的|x|,皮尔逊积分的被积分函数左右各有三个鞍点。

在右平面的鞍点是

rs1 := -(1/3)\sqrt(6)sinh((1/3)arcsinh((3/4)y\sqrt(2)\sqrt(3)/x^(3/2)))

rs2 := (1/3)\sqrt(6)sinh((1/3)arcsinh((3/4)y\sqrt(2)\sqrt(3)/x^(3/2))+(1/3I)Pi)

s3 := -(1/3I)\sqrt(6)cosh((1/3)arcsinh((3/4)y\sqrt(2)\sqrt(3)/x^(3/2))+(1/6I)*Pi)

左平面的鞍点是：

ls1 := -(1/3)\sqrt(6)sin((1/3)arcsin((3/4)y\sqrt(2)\sqrt(3)/x^(3/2))+(1/3)*Pi)

ls3 := (1/3)\sqrt(6)cos((1/3)arcsin((3/4)y\sqrt(2)\sqrt(3)/x^(3/2))+(1/6)*Pi)

ls2 := (1/3)\sqrt(6)sin((1/3)arcsin((3/4)y\sqrt(2)\sqrt(3)/x^(3/2)))

分岔
皮尔西积分的分岔曲线为

27*x^2=-8y^3
斯托克斯曲线
皮尔西积分的斯托克斯曲线为

y^3=\frac{27}{4}(\sqrt(27)-5)x^2

尖点突变
在（x,y)平面中，分岔曲线和斯托克斯曲线将平面分化为尖点突变区。

参考文献
*Frank Oliver, NIST Handbook of Mathematical Functions,2010,Cambridge University Press.
*R.B. Paris,D. Kaminski,Hyperasymptotic evaluation of the Pearcey integral via
Hadamard expansions.Journal of Computational and Applied Mathematics 190 (2006) 437–452