यहां {{math|''μ<sub>j</sub>''}} आव्यूह <math>S_{zz}''(z^0)</math> के [[आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स|आइगेनवैल्यूज़ एवं आइगेनसदिश्स]] हैं।
यहां {{math|''μ<sub>j</sub>''}} आव्यूह <math>S_{zz}''(z^0)</math> के [[आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स|आइगेनवैल्यूज़ एवं आइगेनसदिश्स]] हैं।
{{math proof|title=Proof of complex Morse lemma|proof=
'''एकल अपतित सैडल बिंदु के विषय में स्पर्शोन्मुख विस्तार'''
The following proof is a straightforward generalization of the proof of the real [[Morse theory#Morse lemma|Morse Lemma]], which can be found in.<ref>{{harvtxt|Poston|Stewart|1978}}, page 54; see also the comment on page 479 in {{harvtxt|Wong|1989}}.
</ref> We begin by demonstrating
:'''Auxiliary statement.''' Let {{math| ''f''  : '''C'''<sup>''n''</sup> → '''C'''}} be [[Holomorphic function|holomorphic]] in a neighborhood of the origin and {{math| ''f'' (0) {{=}} 0}}. Then in some neighborhood, there exist functions {{math|''g<sub>i</sub>'' : '''C'''<sup>''n''</sup> → '''C'''}} such that <math display="block">f(z) = \sum_{i=1}^n z_i g_i(z),</math> where each {{math|''g<sub>i</sub>''}} is [[Holomorphic function|holomorphic]] and <math display="block">g_i(0) = \left. \tfrac{\partial f(z)}{\partial z_i} \right|_{z=0}.</math>
From the identity
:<math>f(z) = \int_0^1 \frac{d}{dt} f \left (t z_1,\cdots, t z_n \right ) dt = \sum_{i=1}^n z_i \int_0^1 \left. \frac{\partial f(z)}{\partial z_i}\right|_{z=(t z_1, \ldots, t z_n)} dt,</math>
Without loss of generality, we translate the origin to {{math|''z''<sup>0</sup>}}, such that {{math|''z''<sup>0</sup> {{=}} 0}} and {{math|''S''(0) {{=}} 0}}. Using the Auxiliary Statement, we have
Recall that an arbitrary matrix {{mvar|A}} can be represented as a sum of symmetric {{math|''A''<sup>(''s'')</sup>}} and anti-symmetric {{math|''A''<sup>(''a'')</sup>}} matrices,
Thus, {{math|''h<sub>ij</sub>''(''z'')}} in equation (1) can be assumed to be symmetric with respect to the interchange of the indices {{mvar|i}} and {{mvar|j}}. Note that
hence, {{math|det(''h<sub>ij</sub>''(0)) ≠ 0}} because the origin is a non-degenerate saddle point.
Let us show by [[Mathematical induction|induction]] that there are local coordinates {{math|''u'' {{=}} (''u''<sub>1</sub>, ... ''u<sub>n</sub>''), ''z'' {{=}} '''''ψ'''''(''u''), 0 {{=}} '''''ψ'''''(0)}}, such that
First, assume that there exist local coordinates {{math|''y'' {{=}} (''y''<sub>1</sub>, ... ''y<sub>n</sub>''), ''z'' {{=}} '''''φ'''''(''y''), 0 {{=}} '''''φ'''''(0)}}, such that
where {{math|''H<sub>ij</sub>''}} is symmetric due to equation (2). By a linear change of the variables {{math|(''y<sub>r</sub>'', ... ''y<sub>n</sub>'')}}, we can assure that {{math|''H<sub>rr</sub>''(0) ≠ 0}}. From the [[chain rule]], we have
The matrix {{math|(''H<sub>ij</sub>''(0))}} can be recast in the [[Jordan normal form]]: {{math|(''H<sub>ij</sub>''(0)) {{=}} ''LJL''<sup>−1</sup>}}, were {{mvar|L}} gives the desired non-singular linear transformation and the diagonal of {{mvar|J}} contains non-zero [[Eigenvalues and eigenvectors|eigenvalues]] of {{math|(''H<sub>ij</sub>''(0))}}. If {{math|''H<sub>ij</sub>''(0) ≠ 0}} then, due to continuity of {{math|''H<sub>ij</sub>''(''y'')}}, it must be also non-vanishing in some neighborhood of the origin. Having introduced <math>\tilde{H}_{ij}(y) = H_{ij}(y)/H_{rr}(y)</math>, we write
The change of the variables {{math|''y'' ↔ ''x''}} is locally invertible since the corresponding [[Jacobian matrix and determinant|Jacobian]] is non-zero,
Comparing equations (4) and (5), we conclude that equation (3) is verified. Denoting the [[Eigenvalues and eigenvectors|eigenvalues]] of <math>S''_{zz}(0)</math> by {{math|''μ<sub>j</sub>''}}, equation (3) can be rewritten as
From equation (6), it follows that <math>\det S''_{ww} (\boldsymbol{\varphi}(0)) = \mu_1 \cdots \mu_n</math>. The [[Jordan normal form]] of <math>S''_{zz}(0)</math> reads <math>S''_{zz}(0) = P J_z P^{-1}</math>, where {{math|''J<sub>z</sub>''}} is an upper diagonal matrix containing the [[Eigenvalues and eigenvectors|eigenvalues]] and {{math|det ''P'' ≠ 0}}; hence, <math>\det S''_{zz} (0) = \mu_1 \cdots \mu_n</math>. We obtain from equation (7)
If <math>\det \boldsymbol{\varphi}'_w(0) = -1</math>, then interchanging two variables assures that <math>\det \boldsymbol{\varphi}'_w(0) = +1</math>.
}}'''एकल अपतित सैडल बिंदु के विषय में स्पर्शोन्मुख विस्तार'''
मान लीजिए
मान लीजिए
# {{math| ''f'' (''z'')}} और {{math|''S''(''z'')}} [[ खुला सेट |संवृत, परिबद्ध]],और साधारण रूप से [[ बस जुड़ा हुआ स्थान |जुड़े हुए]] समुच्चय Ωx ⊂ Cn में होलोमोर्फिक फलन हैं जैसे कि Ix = Ωx ∩ Rn जुड़ा हुआ है;
# {{math| ''f'' (''z'')}} और {{math|''S''(''z'')}} [[ खुला सेट |संवृत, परिबद्ध]],और साधारण रूप से [[ बस जुड़ा हुआ स्थान |जुड़े हुए]] समुच्चय Ωx ⊂ Cn में होलोमोर्फिक फलन हैं जैसे कि Ix = Ωx ∩ Rn जुड़ा हुआ है;
Line 162:
Line 79:
यह कथन फेडोर्युक (1987) में प्रस्तुत अधिक सामान्य परिणामों का विशेष विषय है।<ref>{{harvtxt|Fedoryuk|1987}}, pages 417-420.</ref>
यह कथन फेडोर्युक (1987) में प्रस्तुत अधिक सामान्य परिणामों का विशेष विषय है।<ref>{{harvtxt|Fedoryuk|1987}}, pages 417-420.</ref>
{{math proof|title=Derivation of equation (8)|proof=
[[File:Illustration To Derivation Of Asymptotic For Saddle Point Integration.pdf|thumb|center|An illustration to the derivation of equation ((8)]]
First, we deform the contour {{math|''I<sub>x</sub>''}}into a new contour<math>I'_x \subset \Omega_x</math> passing through the saddle point {{math|''x''<sup>0</sup>}} and sharing the boundary with {{math|''I<sub>x</sub>''}}. This deformation does not change the value of the integral {{math|''I''(''λ'')}}. We employ the [[Method of steepest descent#Complex Morse Lemma|Complex Morse Lemma]] to change the variables of integration. According to the lemma, the function {{math|'''''φ'''''(''w'')}} maps a neighborhood {{math|''x''<sup>0</sup> ∈ ''U'' ⊂ Ω''<sub>x</sub>''}} onto a neighborhood {{math|Ω''<sub>w</sub>''}} containing the origin. The integral {{math|''I''(''λ'')}} can be split into two: {{math|''I''(''λ'') {{=}} ''I''<sub>0</sub>(''λ'') + ''I''<sub>1</sub>(''λ'')}}, where {{math|''I''<sub>0</sub>(''λ'')}} is the integral over <math>U\cap I'_x</math>, while {{math|''I''<sub>1</sub>(''λ'')}} is over <math>I'_x \setminus (U\cap I'_x)</math> (i.e., the remaining part of the contour {{math|''I′<sub>x</sub>''}}). Since the latter region does not contain the saddle point {{math|''x''<sup>0</sup>}}, the value of {{math|''I''<sub>1</sub>(''λ'')}} is exponentially smaller than {{math|''I''<sub>0</sub>(''λ'')}} as {{math|''λ'' → ∞}};<ref>This conclusion follows from a comparison between the final asymptotic for {{math|''I''<sub>0</sub>(''λ'')}}, given by equation (8), and [[Method of steepest descent#A simple estimate .5B1.5D|a simple estimate]] for the discarded integral {{math|''I''<sub>1</sub>(''λ'')}}.</ref> thus, {{math|''I''<sub>1</sub>(''λ'')}} is ignored. Introducing the contour {{math|''I<sub>w</sub>''}} such that <math>U\cap I'_x = \boldsymbol{\varphi}(I_w)</math>, we have
Recalling that {{math|''x''<sup>0</sup> {{=}} '''''φ'''''(0)}} as well as <math>\det \boldsymbol{\varphi}_w'(0) = 1</math>, we expand the pre-exponential function <math>f[\boldsymbol{\varphi}(w)]</math> into a Taylor series and keep just the leading zero-order term
Here, we have substituted the integration region {{math|''I<sub>w</sub>''}} by {{math|'''R'''<sup>''n''</sup>}} because both contain the origin, which is a saddle point, hence they are equal up to an exponentially small term.<ref>This is justified by comparing the integral asymptotic over {{math|'''R'''<sup>''n''</sup>}} [see equation (8)] with [[Method of steepest descent#A simple estimate .5B1.5D|a simple estimate]] for the altered part.</ref> The integrals in the r.h.s. of equation (11) can be expressed as
From this representation, we conclude that condition (9) must be satisfied in order for the r.h.s. and l.h.s. of equation (12) to coincide. According to assumption 2, <math>\Re \left( S_{xx}''(x^0) \right)</math> is a [[Definite bilinear form|negatively defined quadratic form]] (viz., <math>\Re(\mu_j)<0</math>) implying the existence of the integral <math>\mathcal{I}_j</math>, which is readily calculated
समीकरण (8) को इस प्रकार भी लिखा जा सकता है{{NumBlk|:|<math>I(\lambda) = \left( \frac{2\pi}{\lambda}\right)^{\frac{n}{2}} e^{\lambda S(x^0)} \left ( \det (-S_{xx}''(x^0)) \right )^{-\frac{1}{2}} \left (f(x^0) + O\left(\lambda^{-1}\right) \right),</math>|{{EquationRef|13}}}}
समीकरण (8) को इस प्रकार भी लिखा जा सकता है{{NumBlk|:|<math>I(\lambda) = \left( \frac{2\pi}{\lambda}\right)^{\frac{n}{2}} e^{\lambda S(x^0)} \left ( \det (-S_{xx}''(x^0)) \right )^{-\frac{1}{2}} \left (f(x^0) + O\left(\lambda^{-1}\right) \right),</math>|{{EquationRef|13}}}}
गणित में, तीव्रतम अवतरण की विधि या सैडल-बिंदु की विधि इंटीग्रल का अनुमान लगाने के लिए लाप्लास की विधि का विस्तार है, जहां स्थिर बिंदु (सैडल बिंदु) के समीप से निकलने के लिए कठोर समतल में समोच्च इंटीग्रल को तीव्रतम अवतरण या स्थिर चरण की दिशा में विकृत किया जाता है। सैडल-पॉइंट सन्निकटन का उपयोग कठोर समतल में इंटीग्रल्स के साथ किया जाता है, जबकि लाप्लास की विधि का उपयोग वास्तविक इंटीग्रल्स के साथ किया जाता है।
अनुमान लगाया जाने वाला इंटीग्रल प्रायः निम्नलिखित रूप का होता है
जहां C समोच्च है, और λ बड़ा है। तीव्रतम अवतरण की विधि का संस्करण एकीकरण C के समोच्च को नवीन पथ एकीकरण C' में विकृत कर देता है जिससे निम्नलिखित स्थितियाँ बनी रहें:
C′ व्युत्पन्न g′(z) के एक या अधिक शून्य से होकर निकलता है,
g(z) का काल्पनिक भाग C′ पर स्थिर है।
तीव्रतम अवतरण की विधि सर्वप्रथम किसके द्वारा प्रकाशित की गई थी? डेबी (1909) harvtxt error: no target: CITEREFडेबी1909 (help), जिन्होंने बेसेल फलन का अनुमान लगाने के लिए इसका उपयोग किया और बताया कि यह हाइपरज्यामितीय फलन के विषय में रीमैन (1863) harvtxt error: no target: CITEREFरीमैन1863 (help) अप्रकाशित नोट में हुआ था। तीव्रतम अवतरण के समोच्च में न्यूनतम गुण होता है, देखें फेडोर्युक (2001) harvtxt error: no target: CITEREFफेडोर्युक2001 (help) देखें। सीगल (1932) harvtxt error: no target: CITEREFसीगल1932 (help) रीमैन के कुछ अन्य अप्रकाशित नोट्स का वर्णन किया, जहां उन्होंने रीमैन-सीगल सूत्र प्राप्त करने के लिए इस विधि का उपयोग किया था।
तीव्रतम अवतरण की विधि प्रपत्र के कठोर इंटीग्रल का अनुमान लगाने की विधि है
बड़े के लिए, जहाँ और , के विश्लेषणात्मक फलन हैं। क्योंकि इंटीग्रैंड विश्लेषणात्मक है, रूपरेखा नये स्वरूप में इंटीग्रल को परिवर्तित किए बिना विकृत किया जा सकता है। विशेष रूप से, कोई नई रूपरेखा की शोध करता है जिस पर काल्पनिक भाग हो स्थिर है। तब
और शेष इंटीग्रल का अनुमान लाप्लास की विधि जैसी अन्य विधियों से लगाया जा सकता है।[1]
व्युत्पत्ति
विश्लेषणात्मक होने के कारण इस विधि को तीव्रतम अवतरण की विधि कहा जाता है, स्थिर चरण समोच्च तीव्रतम अवरोही समोच्चों के समतुल्य हैं।
यदि का विश्लेषणात्मक फलन है, यह कॉची-रीमैन समीकरण को संतुष्ट करता है,
तब
इसलिए स्थिर चरण की आकृतियाँ भी तीव्रतम अवतरण की आकृतियाँ हैं।
साधारण अनुमान
मान लीजिए f, S : Cn → C और C ⊂ Cn, यदि
जहाँ वास्तविक भाग को दर्शाता है, और धनात्मक वास्तविक संख्या λ0 सम्मिलित है जो इस प्रकार है,
अपतित सैडल बिंदु, z0 ∈ Cn, होलोमोर्फिक फलन S(z) का महत्वपूर्ण बिंदु है (अर्थात, ∇S(z0) = 0) जहां फलन के हेसियन आव्यूह में अलुप्त होने वाला निर्धारक (अर्थात, ) है।
गैर-अपक्षयी सैडल बिंदु के विषय में इंटीग्रल के एसिम्प्टोटिक्स के निर्माण के लिए निम्नलिखित मुख्य उपकरण है:
कॉम्प्लेक्स मोर्स लेम्मा
वास्तविक-मूल्यवान फलन के लिए मोर्स लेम्मा होलोमोर्फिक फलन के लिए निम्नानुसार सामान्यीकृत करता है।[3] होलोमोर्फिक फलन S(z) के अपतित सैडल बिंदु z0 के पास, ऐसे निर्देशांक होते हैं जिनके संदर्भ में S(z) − S(z0) सम्पूर्ण द्विघात है। इसे त्रुटिहीन बनाने के लिए S डोमेन W ⊂ Cn के साथ होलोमोर्फिक फलन मान लीजिए, और W में z0 को S का अपतित सैडल बिंदु मान लीजिए, अर्थात, ∇S(z0) = 0 और , फिर z0 के नेबर U ⊂ W और w = 0 के V ⊂ Cn और φ(0) के साथ विशेषण होलोमोर्फिक फलन सम्मिलित है, φ: V → U φ : V → U साथ φ(0) = z0 इस प्रकार है कि
एकल अपतित सैडल बिंदु के विषय में स्पर्शोन्मुख विस्तार
मान लीजिए
f (z) और S(z)संवृत, परिबद्ध,और साधारण रूप से जुड़े हुए समुच्चय Ωx ⊂ Cn में होलोमोर्फिक फलन हैं जैसे कि Ix = Ωx ∩ Rn जुड़ा हुआ है;
के x0 ∈ Ix सम्पूर्ण बिंदु के लिए एकल अधिकतम है;
x0 अपतित सैडल बिंदु (अर्थात, ∇S(x0) = 0 और ) है,
फिर, निम्नलिखित स्पर्शोन्मुख धारण करता है,
(8)
जहाँ μj हेस्सियन आव्यूह और के आइगेनवैल्यू हैं जो नियमों से परिभाषित किये गये हैं,
(9)
यह कथन फेडोर्युक (1987) में प्रस्तुत अधिक सामान्य परिणामों का विशेष विषय है।[4]
समीकरण (8) को इस प्रकार भी लिखा जा सकता है
(13)
निम्नानुसार चयन किया गया है,
महत्वपूर्ण विषयों पर विचार करें:
यदि S(x), Rn (अर्थात, बहुआयामी लाप्लास विधि) में वास्तविक x और x0 के लिए वास्तविक मूल्य है, फिर[5]
यदि S(x)x के लिए वास्तव में पूर्णतया काल्पनिक है (अर्थात, सभी के लिए x में Rn) और x0 में Rn (अर्थात, बहुआयामी स्थिर चरण विधि),[6] तब[7]
जहाँ आव्यूह के जड़त्व के नियम को दर्शाता है। प्रमेय का कथन , जो ऋणात्मक आइगेनवैल्यू की संख्या घटाकर धनात्मक आइगेनवैल्यू की संख्या के समान है। यह उल्लेखनीय है कि क्वांटम यांत्रिकी (साथ ही प्रकाशिकी में) में बहुआयामी WKB सन्पासन के लिए स्थिर चरण विधि के अनुप्रयोगों में, Ind मास्लोव सूचकांक से संबंधित है, उदाहरण के लिए, चाइचियन & डेमीचेव (2001) harvtxt error: no target: CITEREFचाइचियनडेमीचेव2001 (help) और शुलमैन (2005) harvtxt error: no target: CITEREFशुलमैन2005 (help) है।
एकाधिक अक्षतिग्रस्त सैडल बिंदुओं का विषय
यदि फलन S(x) में कई भिन्न-भिन्न अपतित सैडल बिंदु हैं, अर्थात,
जहाँ
Ωx का संवृत आवरण है, तो एकता के विभाजन को नियोजित करके इंटीग्रल एसिम्प्टोटिक की गणना को एकल सैडल बिंदु के विषय में कम कर दिया जाता है। एकता का विभाजन हमें निरंतर फलन ρk(x) : Ωx → [0, 1], 1 ≤ k ≤ K, का समुच्चय बनाने की अनुमति देता है जो इस प्रकार है,
जहाँ से,
इसलिए जैसे λ → ∞ हमारे पास है:
जहां अंतिम चरण में समीकरण (13) और पूर्व-घातीय फलन का उपयोग किया गया था f (x) कम से कम निरंतर होना चाहिए।
अन्य विषय
जब ∇S(z0) = 0 और , बिंदु z0 ∈ Cn को किसी फलन S(z) का अपतित सैडल पॉइंट कहा जाता है।
स्पर्शोन्मुख की गणना
जब λ → ∞, f (x) सतत है, और S(z) में पतित सैडल बिंदु है, यह अधिक समृद्ध समस्या है, जिसका समाधान अधिकतम सीमा तक कैटास्ट्रोफ सिद्धांत पर निर्भर करता है। यहां, कैटास्ट्रोफ सिद्धांत मोर्स लेम्मा की विधि को, S(z) विहित अभ्यावेदन से एक में परिवर्तित करने के लिए प्रतिस्थापित करता है जो केवल अपतित विषय में मान्य है। अधिक जानकारी के लिए, उदाहरण, पोस्टन & स्टीवर्ट (1978) harvtxt error: no target: CITEREFपोस्टनस्टीवर्ट1978 (help) और फेडोर्युक (1987) harvtxt error: no target: CITEREFफेडोर्युक1987 (help) है।
विकृत सैडल बिंदुओं वाले इंटीग्रल स्वाभाविक रूप से कास्टिक (प्रकाशिकी) और क्वांटम यांत्रिकी में बहुआयामी डब्ल्यूकेबी सहित कई अनुप्रयोगों में दिखाई देते हैं।
अन्य विषय जैसे f (x) और S(x) असंतत हैं या जब S(x) का शीर्ष एकीकरण क्षेत्र की सीमा पर स्थित है,तो विशेष सुरक्षा की आवश्यकता है, उदाहरण के लिए, फेडोर्युक (1987) harvtxt error: no target: CITEREFफेडोर्युक1987 (help) और वोंग (1989) harvtxt error: no target: CITEREFवोंग1989 (help) है।
विस्तार और सामान्यीकरण
सबसे तीव्र अवतरण विधि का विस्तार तथाकथित अरेखीय स्थिर चरण/सबसे तीव्र अवतरण विधि है। यहां, इंटीग्रल के अतिरिक्त, किसी को रीमैन-हिल्बर्ट फ़ैक्टराइज़ेशन समस्याओं के स्पर्शोन्मुख समाधानों का मूल्यांकन करने की आवश्यकता है।
कठोर क्षेत्र में समोच्च C को देखते हुए, उस समोच्च पर परिभाषित फलन f और विशेष बिंदु, मान लीजिए अनंत, समोच्च C से दूर फलन f C में निर्धारित जंप के साथ, और अनंत पर दिए गए सामान्यीकरण के साथ होलोमोर्फिक की शोध करता है। यदि f और इसलिए M अदिश के अतिरिक्त आव्यूह हैं तो यह ऐसी समस्या है जो सामान्य रूप से स्पष्ट समाधान स्वीकार नहीं करती है।
तब रैखिक स्थिर चरण/तीव्रतम अवतरण विधि की के विषय पर स्पर्शोन्मुख मूल्यांकन संभव है। विचार यह है कि दी गई रीमैन-हिल्बर्ट समस्या के समाधान को असम्बद्ध रूप से कम करके सरल, स्पष्ट रूप से समाधान करने योग्य, रीमैन-हिल्बर्ट समस्या बना दिया जाए। कॉची के प्रमेय का उपयोग जम्प समोच्च की विकृतियों को उचित बताने के लिए किया जाता है।
रूसी गणितज्ञ अलेक्जेंडर इट्स के पूर्व फलन के आधार पर, 1993 में डेफ्ट और झोउ द्वारा नॉनलाइनियर स्थिर चरण का प्रारम्भ किया गया था। लैक्स, लीवरमोर, डेफ्ट, वेनाकिड्स और झोउ के पूर्व फलन के आधार पर, 2003 में फलनविसिस, के. मैकलॉघलिन और पी. मिलर द्वारा नॉनलाइनियर स्टीपेस्ट डीसेंट विधि प्रस्तुत की गई थी। जैसा कि रैखिक विषय में होता है, सबसे तीव्र अवरोही आकृतियाँ न्यूनतम-अधिकतम समस्या का समाधान करती हैं। अरैखिक विषय में वे S-वक्र बन जाते हैं (80 के दशक में स्टाल, गोन्चर और राखमनोव द्वारा भिन्न संदर्भ में परिभाषित)।
↑Rigorously speaking, this case cannot be inferred from equation (8) because the second assumption, utilized in the derivation, is violated. To include the discussed case of a purely imaginary phase function, condition (9) should be replaced by
Chaichian, M.; Demichev, A. (2001), Path Integrals in Physics Volume 1: Stochastic Process and Quantum Mechanics, Taylor & Francis, p. 174, ISBN075030801X
Deift, P.; Zhou, X. (1993), "A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics for the MKdV equation", Ann. of Math., The Annals of Mathematics, Vol. 137, No. 2, vol. 137, no. 2, pp. 295–368, arXiv:math/9201261, doi:10.2307/2946540, JSTOR2946540, S2CID12699956.
Fedoryuk, M. V. (1987), Asymptotic: Integrals and Series, Nauka, Moscow [in Russian].
Kamvissis, S.; McLaughlin, K. T.-R.; Miller, P. (2003), "Semiclassical Soliton Ensembles for the Focusing Nonlinear Schrödinger Equation", Annals of Mathematics Studies, Princeton University Press, vol. 154.
Riemann, B. (1863), Sullo svolgimento del quoziente di due serie ipergeometriche in frazione continua infinita (Unpublished note, reproduced in Riemann's collected papers.)
Siegel, C. L. (1932), "Über Riemanns Nachlaß zur analytischen Zahlentheorie", Quellen und Studien zur Geschichte der Mathematik, Astronomie und Physik, 2: 45–80 Reprinted in Gesammelte Abhandlungen, Vol. 1. Berlin: Springer-Verlag, 1966.
Translated in Deift, Percy; Zhou, Xin (2018), "On Riemanns Nachlass for Analytic Number Theory: A translation of Siegel's Uber", arXiv:1810.05198 [math.HO].
Poston, T.; Stewart, I. (1978), Catastrophe Theory and Its Applications, Pitman.
Schulman, L. S. (2005), "Ch. 17: The Phase of the Semiclassical Amplitude", Techniques and Applications of Path Integration, Dover, ISBN0486445283
Wong, R. (1989), Asymptotic approximations of integrals, Academic Press.