text:
# Asymptotic behaviour of a solution for Kadomtsev-Petviashvili-2 equation This work was partially supported by RFBR (97-01-459), grant for scientific schools (001596038) and INTAS (No1068).

## 1 Introduction

Properties of solutions and process of integration for the Kadomtsev-Petviashvili equation (KP) :

$$_x(_tu+6u_xu+_x^3u)=3\sigma ^2_y^2u$$
(1)
depend on sign of $`\sigma ^2`$. This equation is called KP-2 if $`\sigma ^2=1`$.

In this work the asymptotic behaviour of decreasing solution of equation KP-2 is obtained as $`t\mathrm{}`$. The main term of the asymptotics has an order by $`O(t^1)`$ and fast oscillates. An envelope of these oscillations depends on $`\xi =x/t`$ and $`\eta =y/t`$. Results of the work are formulated in terms of scattering data for auxiliary linear problem, which is associated with equation KP-2 in the inverse scattering transform (IST) .

The equation KP plays important role in modern mathematical physics. Many applications of these equation are known in plasma physics, water waves and other fields of waves process . Therefore questions about solvability of this equation in different functional classes were studied in detail. In particular, the existence of global solution of Cauchy problem in class of distribution functions was proved in . A norm of solution in Sobolev space was estimated and an asymptotics of a solution were obtained as $`t\mathrm{}`$ for KP-like equation but with more high nonlinearity, which called generalized equation KP, in .

The IST formalism allows to reduce a constructing of the solution of nonlinear integrable equation into solving of linear problems. One of most important achievements of IST is the constructing of asymptotic behaviours of solutions as $`t\mathrm{}`$. These results are well-known for 1+1-dimensional equations (one spatial and one temporal variable) -. Asymptotic behaviours of solutions for (2+1)-dimensional equations were studied not so in detail. Rigorous results about temporal asymptotics of solutions, which is nonuniform with respect to spatial variables, for a special class of nondecreasing solutions were obtained in , . In work  an formal asymptotics of decreasing solution of equation KP-1, which is nonuniform with respect to spatial variables, was constructed. Results about asymptotic behaviour of decreasing solutions of equation KP-2 as $`t\mathrm{}`$ and about a remainder of this asymptotics, which are uniform with respect to spatial variables, are obtained in present work.

The formalism IST for solving of Cauchy problem for equation KP-2 was presented in . This formalism will be used in this work. Below we remind basic steps of solving of Cauchy problem for equation KP-2 by IST.

Let us denote an initial condition for equation KP-2 as:

$$u|_{t=0}=u_0(x,y).$$
(2)
The process of solving of the Cauchy problem (1), (2) consists of several steps.

First step is solving direct scattering problem. On this step a boundary problem is solved for function $`\phi `$:

$$_y\phi +_x^2\phi +2ik_x\phi +u\phi =0,\phi |_{|k|\mathrm{}}=1,$$
(3)
and scattering data are constructed by a formula:

$`F(k)=`$ $`{\displaystyle \frac{\text{sgn}(\text{Re}(k))}{2\pi }}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)\phi (x,y,k,0)\times `$ (4)
$`\times \mathrm{exp}(i(k+\overline{k})x(k^2\overline{k}^2)y).`$
It is useful to note, if $`u`$ is real, then the function $`F(k)`$ has property $`F(\overline{k})=\overline{F}(k)`$. This follows from formulas (3) and (4).

On the next step an evolution of scattering data is determined. This evolution is very simple:

$$(k;t)=F(k)\mathrm{exp}(4it(k^3+\overline{k}^3)).$$

Third step is solving of inverse scattering problem. This is reduced to so-called $`\overline{D}`$-problem:

$$\begin{array}{c}_{\overline{k}}\phi =\psi F(\overline{k})\mathrm{exp}(itS),\\ _k\psi =\phi F(k)\mathrm{exp}(itS);\end{array}\left(\begin{array}{c}\phi \\ \psi \end{array}\right)|_{|k|\mathrm{}}=\left(\begin{array}{c}1\\ 1\end{array}\right),$$
(5)
where $`S=4(k^3+\overline{k}^3)+(k+\overline{k})\xi i(k^2\overline{k}^2)\eta ,\xi =x/t,\eta =y/t`$, the function $`F(k)`$ is nonanalytic with respect to complex variable $`k\text{}`$. Solving of this problem allows to obtain the functions $`\phi `$ and $`\psi `$ at any time, if we know the evolution of scattering data.

At last we can obtain the solution of the Cauchy problem by using the formula :

$$u(x,y,t)=_x_{\text{}}𝑑kd\overline{k}F(k)\psi (k,x,y,t)\mathrm{exp}(itS).$$
(6)
However, this successive method of solving for Cauchy problems for general initial conditions allows to obtain very implicit answer with respect to initial data. Therefore the problem (5) and integral (6) are used usually for solving of integrable nonlinear equations and the scattering data is taking from some suitable functional class (see, for example, -). On one-hand side this gives for studying a functional class of solutions, but on the other hand-side this leads to implicit restrictions on class of initial conditions.

The auxiliary linear problem (3) and some related problems were studied in -. In the work  it was proved that if the function $`u_0(x,y)`$ decreases exponentially, then the direct scattering problem (3) is solvable for $`k\text{}`$.

To construct the asymptotics of $`u`$ as $`t\mathrm{}`$ we use the asymptotics of $`\psi `$ and evaluate an asymptotics of the integral (6) by stationary phase method . As it turned out, the main term of asymptotics of the function $`u`$ may be defined by using only the main term of an asymptotic expansion of $`\psi `$ as $`t\mathrm{}`$. Taking into account the fast oscillating coefficients of the system from (5) one can guess, that $`\psi =1+o(1)`$ as $`t\mathrm{}`$. But we cannot say something definitely about the asymptotic behaviour of the function $`u`$, until we do not know analytic properties and an order of the remainder of asymptotics of $`\psi `$ more precisely. Thus we come to studying an asymptotic behaviour of a solution of the $`\overline{D}`$-problem (5).

The asymptotic behaviour of the solution of the $`\overline{D}`$-problem as $`t\mathrm{}`$ with continuous and fast oscillated coefficients was obtained in . Here we study the $`\overline{D}`$-problem with discontinuous coefficients on imaginary axis of complex parameter $`k`$. Constructing of the asymptotics of solution for such problem is more complicated. Not only stationary points of phase function of an oscillated exponent define a structure of the asymptotic solution for (5) as in  (see also ), but the location of these stationary points with respect to line of discontinuity of coefficients of the equations (5) as well. It contributes additional difficulties into evaluating and leads to changes in results. The uniform asymptotics of the solution of the problem (5) is constructed by matching method .

## 2 Main result

###### Theorem 1

Let $`(1+|k|)FL_1C`$, $`^\alpha FL_1C`$ as $`\text{Re}(k)0`$, $`|a|2`$ and:

$$\underset{z\text{}}{sup}_\text{}^2𝑑\kappa 𝑑\lambda \left|\frac{F(\kappa +i\lambda )}{\kappa +i\lambda z}\right|<\mathrm{\hspace{0.17em}2}\pi ,$$
(7)
then the solution of the Cauchy problem of equation KP-2 for corresponding initial condition exists as $`t>0`$. The asymptotic behaviour of the solution as $`t\mathrm{}`$ differs in different domains of variables $`(x,y,t)`$:
as $`(12\xi +\eta ^2)t^{1/3}1`$:

$`u(x,y,t)=4t^1{\displaystyle \frac{\pi }{12i\sqrt{\eta ^212\xi }}}f({\displaystyle \frac{1}{2}}\sqrt{\eta ^212\xi }+{\displaystyle \frac{i\eta }{12}})\times `$
$`\times \mathrm{exp}(11it\sqrt{{\displaystyle \frac{y^2}{t^2}}12{\displaystyle \frac{x}{t}}})+c.c.+o(1).`$
as $`(12\xi +\eta ^2)t^{1/3}1`$:

$$u=o(t^1);$$

as $`|12\xi +12\eta ^2|1`$:

$`u(x,y,t)=8it^1\sqrt{\pi }f(i\eta /12)({\displaystyle _0^{\mathrm{}}}dp_1\sqrt{p_1}\mathrm{cos}(8p_1^3zp_1)+`$
$`+{\displaystyle _0^{\mathrm{}}}dp_1\sqrt{p_1}\mathrm{sin}(8p_1^3zp_1))+o(t^1).`$
Here $`\xi =x/t`$, $`\eta =y/t`$,

$$z=8\left(\frac{y^2}{12t^{4/3}}+\frac{x}{t^{1/3}}\right);$$

$`f(k)={\displaystyle \frac{1}{2\pi }}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)\phi (x,y,k,0)\times `$
$`\times \mathrm{exp}(i(k+\overline{k})x(k^2\overline{k}^2)y).`$
The domains of validity for the asymptotics of the solution of equation KP-2 are intersected and therefore they give combined asymptotics of the solution uniformly on plane of $`x,y`$.

## 3 An analytic behaviours of the scattering data

In this section we demonstrate analytic behaviour of the scattering data which corresponds to sufficiently smooth and decreasing initial condition (2).

First of all we show, that the solution of (5) exists.

###### Theorem 2

Let $`F(k)`$ be such, that the condition (7) is fulfilled, then the solution of (5) exists in a space of continuous vector-functions bounded when $`k\text{}`$.

Proof. Let us consider a system of integral equations which is equivalent to (5):

$$\left(\begin{array}{c}\phi \\ \psi \end{array}\right)=\left(\begin{array}{c}1\\ 1\end{array}\right)+G[F]\left(\begin{array}{c}\phi \\ \psi \end{array}\right),$$

where

$`G[F]V={\displaystyle }{\displaystyle _m\text{}}dmd\overline{m}\times `$
$`\left(\begin{array}{cc}0& \frac{F(\overline{m})}{km}\mathrm{exp}(itS)\\ \frac{F(m)}{\overline{km}}\mathrm{exp}(itS)& 0\end{array}\right)V(m,\xi ,\eta ,t).`$
Using well-known results about integral operators (see, for example ) one can show that the operator $`G[F]`$ transforms the space of continuous vector-functions into itself.

The operator $`G[F]`$ is contracting operator. It is follows from inequality (7). Hence we obtain the theorem statement.

To construct the asymptotics we must study smoothness of scattering data $`F(k)`$ in neighborhood of line of discontinuity $`\text{Re}(k)=0`$.

###### Lemma 1

Let $`u_0(x,y)`$ be a finite function, then in neighborhood of some point $`k^{}\text{}`$ the scattering data have to be represented in the form:

$$F(k)=\text{sgn}[\text{Re}(k)]\underset{|\alpha |=0}{\overset{2}{}}f_{\alpha _1\alpha _2}(k^{})(kk^{})^{\alpha _1}\overline{(kk^{})^{\alpha _2}}+O(|kk^{}|^3),$$

where $`f_{\alpha _1,\alpha _2}(k^{})`$ are continuous functions, when $`\text{Re}(k^{})0`$ and

$$f_{\alpha _1,\alpha _2}(k^{})=f_{\alpha _1\alpha _2}^{(1)}(k^{})+\text{sgn}[\text{Re}(k)]f_{\alpha _1\alpha _2}^{(2)}(k^{}),$$

when $`\text{Re}(k^{})=0`$.

Proof of this lemma consists of successive evaluating of partial derivatives of scattering data with respect to $`k`$ and $`\overline{k}`$. For example let us evaluate $`f_{10}^{(j)}(k^{})`$ and $`f_{10}^{(j)}(k^{})`$ as $`\text{Re}(k^{})=0`$.

Denote

$$f(k)=\frac{1}{2\pi }_\text{}^2𝑑x𝑑yu_0(x,y)\phi (x,y,k,0)\mathrm{exp}(i(k+\overline{k})x(k^2\overline{k}^2)y).$$

Then $`f_{00}^{(1)}(k^{})=f(k^{})`$, $`f_{00}^{(0)}(k^{})0`$.

Evaluate $`f_{10}(k^{})=f_{10}^{(1)}(k^{})+\text{sgn}[\text{Re}(k)]f_{10}^{(0)}(k^{})`$.

$`_kf(k)={\displaystyle \frac{1}{2\pi }}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)\mathrm{exp}(i(k+\overline{k})x(k^2\overline{k}^2)y)\times `$
$`\left[(ix2ky)\phi (x,y,k,0)+_k\phi (x,y,k,0)\right].`$
To evaluate the derivative $`_k\phi `$ we use the integral equation for the function $`\phi `$:

$`_k\phi ={\displaystyle \frac{1}{2i\pi }}V.P.{\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dmd\overline{m}}{(km)^2}}F(\overline{m})\times `$
$`\mathrm{exp}(i(m+\overline{m})(m^2\overline{m}^2)y)\psi (x,y,m).`$
The functions $`F(\overline{m})`$ and $`\psi (x,y,m)`$ are smooth on the left and right-hand sides of complex plane of $`m`$ and one can show that the integral exists. Let us represent it into more convenient form. For this we rewrite the integral into sum of two integrals over left-hand side and right-hand side of complex plane. Let us integrate by parts these integrals. As a result we obtain

$`_k\phi ={\displaystyle \frac{1}{2i\pi }}{\displaystyle \underset{\pm }{}}{\displaystyle _{\mathrm{\Omega }^\pm }}{\displaystyle \frac{(\pm )d\overline{m}}{km}}[f(\overline{m})\times `$
$`\mathrm{exp}(i(m+\overline{m})x(m^2\overline{m}^2)y)\psi (x,y,m)]`$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle \underset{\pm }{}}{\displaystyle }{\displaystyle _{\mathrm{\Omega }^\pm }}{\displaystyle \frac{dmd\overline{m}}{km}}[(ix2my)(\pm )f(\overline{m})\psi (x,y,m)+`$
$`(\pm )_mf(\overline{m})\psi (x,y,m)\phi (x,y,m)f(\overline{m})f(m)]\times `$
$`\times \mathrm{exp}(i(m+\overline{m})x(m^2\overline{m}^2)y).`$
Here $`\mathrm{\Omega }^\pm `$ is right-hand side (+) and left-hand side (-) of complex plane. We mean the integral over $`\mathrm{\Omega }^+`$ as sum of integral over right-hand side and left-hand side of the circle with center at origin of coordinates and with large radius (which tends into infinity) and of improper integral over imaginary line and of an integral over the circle with small radius $`\epsilon `$ with center at $`k=m`$. The direction of the integration over $`\mathrm{\Omega }^+`$ is determined by standard way.

Consider the sum of the integrals over $`\mathrm{\Omega }^\pm `$. Each of the integrals over large half-circles is equal to zero as $`R\mathrm{}`$ because the integrand decreases. The sum of the integrals over the imaginary axis gives doubled integral over the imaginary axis in positive direction. The integral over small circle equals zero as $`\epsilon 0`$.

Let us evaluate the derivative $`_mf(\overline{m})`$ of the integrand of the double integrals. For this we note that $`\phi (x,y,k,0)=\psi (x,y,\overline{k},0)`$, hence $`f(\overline{m})`$ has to be represented by the function $`\psi (x,y,m,0)`$ and then:

$$_mf(\overline{m})=\frac{1}{2\pi }_\text{}^2dxdyu_0(x,y)[f(m)\phi (x,y,m,0)+$$

$$+(ix2my)\psi (x,y,m,0)]\mathrm{exp}(i(m+\overline{m})x(m^2\overline{m}^2)y).$$

So we obtain a formula $`f_{10}(k^{})`$:

$$f_{10}(k^{})=f_{10}^{(1)}(k^{})+\text{sgn}[\text{Re}(k)]f_{10}^{(2)}(k^{}),$$

where

$`f_{10}^{(1)}(k^{})=f(\overline{k}^{})\psi (x,y,k^{},0){\displaystyle \frac{1}{2\pi i}}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)\times `$
$`[(ix2ky)\varphi (x,y,k,0)+`$
$`({\displaystyle \frac{1}{2i\pi }}{\displaystyle \underset{\pm }{}}V.P.{\displaystyle _i\mathrm{}^i\mathrm{}}{\displaystyle \frac{d\overline{m}}{k^{}m}}(\pm )f(\overline{m})\psi (x,y,m,0)`$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle \underset{\pm }{}}{\displaystyle }{\displaystyle _{\mathrm{\Omega }^\pm }}{\displaystyle \frac{dmd\overline{m}}{km}}[(ix2my)(\pm )f(\overline{m})\psi (x,y,m)+`$
$`(\pm )_mf(\overline{m})\psi (x,y,m)\phi (x,y,m)f(\overline{m})f(m)]\times `$
$`\mathrm{exp}(i(m+\overline{m})x(m^2\overline{m}^2)y))\times `$
$`\mathrm{exp}(i(k^{}+\overline{k}^{})x(k^2\overline{k}^2)y)];`$
$$f_{10}^{(2)}=\frac{1}{\pi }_\text{}^2𝑑x𝑑yu_0(x,y)f(\overline{k}^{})\psi (x,y,k^{},0).$$

An expression of $`f_{01}(k^{})=f_{01}^{(1)}(k^{})+\text{sgn}[\text{Re}(k)]f_{01}^{(2)}(k^{})`$ has the form:

$`f_{01}^{(1)}(k^{})`$ $`={\displaystyle \frac{1}{2i\pi }}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)(ix+2\overline{k}^{}y)\phi (x,y,k^{},0)\times `$
$`\mathrm{exp}(i(k^{}+\overline{k}^{})x(k^2\overline{k}^2)y),`$
$`f_{01}^{(2)}(k^{})`$ $`={\displaystyle \frac{1}{2i\pi }}{\displaystyle }{\displaystyle _\text{}^2}dxdyu_0(x,y)\psi (x,y,k^{},0)\times `$
$`\times \mathrm{exp}(2i(k^{}+\overline{k}^{})x2(k^2\overline{k}^2)y).`$
The expressions of other coefficients of the expansion have to be evaluated by the same way. The lemma is proved.

## 4 Asymptotic solution of $`\overline{D}`$-problem

In this section we construct an asymptotic solution as $`t\mathrm{}`$ of the $`\overline{D}`$-problem:

$$\begin{array}{c}_{\overline{k}}\mu =\nu F(k)\mathrm{exp}(itS),\\ _k\nu =\mu F(\overline{k})\mathrm{exp}(itS);\end{array}\left(\begin{array}{c}\mu \\ \nu \end{array}\right)|_{|k|\mathrm{}}=\left(\begin{array}{c}1\\ 0\end{array}\right).$$
(9)
The solution of the problem (5) has to be obtained using the solution of (9) and the formula:

$$\left(\begin{array}{c}\phi \\ \psi \end{array}\right)=\left(\begin{array}{c}\mu (k,\xi ,\eta ,t)\\ \nu (k,\xi ,\eta ,t)\end{array}\right)+\left(\begin{array}{c}\nu (\overline{k},\xi ,\eta ,t)\\ \mu (\overline{k},\xi ,\eta ,t)\end{array}\right).$$
(10)
The asymptotic solution of the problem (9) is combined. Here we construct the asymptotics as $`t\mathrm{}`$ uniformly with respect to all parameters. The stationary points of the functions $`S(k,\overline{k},\xi ,\eta )`$ with respect to parameters $`k,\overline{k}`$ plays the important role in these constructions.

The asymptotic expansion of the solutions is constructed using the inverse powers of large parameter $`t`$ as a asymptotic sequence of asymptotic expansion outside of small neighborhoods of the stationary points of the function. Near the stationary points the asymptotic sequence has the form $`t^{n/2}`$, where $`n=0,1,2,\mathrm{}`$. Following the terminology of matching method , we call the asymptotics outside of small neighborhoods of stationary points by outer expansion and we call by interior expansion the asymptotics near the stationary points . The right-hand side in the system of equations (9) is discontinuous on the line $`\text{Re}(k)=0`$, therefore the asymptotic solution is differentiable out of the line $`\text{Re}(k)=0`$.

The domains of validity for interior and exterior asymptotic expansions are intersected. This fact is used by matching method in order to construct unique combined asymptotic expansions.

The phase function $`S`$ depends on two parameters $`(\xi ,\eta )`$ $`\text{}^2`$. On the curve $`12\xi +\eta ^2=0`$ the confluence of two stationary points of the function $`S`$ occurs. In this case we have one confluent stationary point. The structure of the asymptotic expansion of solution of (9) is changed here. As $`|12\xi +\eta ^2|1`$ the expansion is constructed on the powers of $`t^{n/3},`$ $`n=0,1,2,\mathrm{}`$ as a asymptotic sequence.

The formulas for the uniform asymptotic expansions of the solution of (9) are large and it seems convenient to formulate the results about general case in section 4.1 and about confluence case in section 4.2 for convenience .

### 4.1 Asymptotics in a general case

In this section the combined asymptotic solution of (9) is constructed when the phase function $`S`$ has nondegenerate stationary points $`k=k_1`$ and $`k=k_2`$. Here we suppose that $`k\text{},t^{1/3}\left|_k^2S|_{k=k_{1,2}}\right|1`$ as $`t\mathrm{}`$. This leads to restriction on values of the parameters $`\xi `$ and $`\eta `$, namely, $`t^{1/3}|12\xi +\eta ^2|1`$. The asymptotic expansion which is uniform with respect to $`k\text{}`$ is formulated in the end of this section.

To construct the combined asymptotic solution which is valuable as $`k\text{}`$ we obtain exterior and interior asymptotic expansions. These expansions are valid outside of small neighborhoods of $`k_j`$, $`j=1,2`$ and in the small neighborhoods of $`k_j`$ respectively. Denote: $`\theta =\sqrt{12\xi \eta ^2}`$. Then one can obtain an expressions for the stationary points of $`S`$ by using the parameter $`\theta `$: $`k_1=\frac{1}{12}(i\eta +\theta ),k_2=\frac{1}{12}(i\eta \theta )`$.

#### 4.1.1 The stationary points outside of the break line

Consider the case when $`\text{Re}(\theta )0`$, i.e. when the stationary points $`k_{1,2}`$ are outside of the discontinuity line $`\text{Re}(k)=0`$. Let us formulate a result of this section about the combined asymptotic solution:

###### Lemma 2

Let the system of the equations (9) have not the homogeneous solutions, $`F(k)C^2L_1`$ as $`\text{Re}(k)0`$ and the parameters $`\xi `$ and $`\eta `$ satisfy the inequality $`t^{1/3}|\theta |^21`$, then:
when $`\sqrt{t|\theta |}|kk_{1,2}|1`$ the formal asymptotic solution of the system (9) with respect to $`mod(O(t^2|_kS|^3))`$ has the form:

$$\stackrel{~}{\mu }=1+t^1\stackrel{1}{\mu }(k,\xi ,\eta ),$$

$$\stackrel{~}{\nu }=(t^1\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )+t^2\underset{1}{\overset{2}{\nu }}(k,\xi ,\eta ))\mathrm{exp}(itS)+t^1\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta );$$

the functions $`\stackrel{1}{\mu },\underset{1}{\overset{1}{\nu }},\underset{1}{\overset{2}{\nu }},\underset{0}{\overset{1}{\nu }}`$ are defined by (19), (13), (17), (30);
when $`|\theta |^1|kk_j|1`$ the formal asymptotic solution of the system (9) with respect to $`mod(O(t^1|\theta |^1))`$ has the form:

$$\stackrel{~}{\mu }=1+t^1\stackrel{1}{M}(l_j,\xi ,\eta ),$$

$$\stackrel{~}{\nu }=\left(t^{1/2}\stackrel{1}{N}(l_j,\xi ,\eta )+t^1\stackrel{2}{N}(l_j,\xi ,\eta )\right)\mathrm{exp}(itS),$$

where $`l_j,j=1,2`$, are defined by formula:

$$l_j=\sqrt{t}\left(kk_j\right)\sqrt{\frac{_k^2S_j}{2}+\mathrm{\hspace{0.17em}4}(kk_j)},$$

the functions $`\stackrel{1}{M}(l_j,\xi ,\eta ),\stackrel{1}{N}(l_j,\xi ,\eta ),\stackrel{2}{N}(l_j,\xi ,\eta )`$ are defined by (32), (27), (31).

Proof. Let us construct the external asymptotic solution in the form:

$$\mu ^{ex}=1+t^1\stackrel{1}{\mu }(k,\xi ,\eta )+t^2\underset{1}{\overset{2}{\mu }}(k,\xi ,\eta )\mathrm{exp}(itS)+\mathrm{},$$
(11)
$`\nu ^{ex}=\left(t^1\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )+t^2\stackrel{2}{\nu }(k,\xi ,\eta )+\mathrm{}\right)\mathrm{exp}(itS)+`$
$`+t^1\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta )+\mathrm{}.`$ (12)
Let us substitute the formulas (11) and (12) into (9), equate coefficients with equal power of $`t`$. As a result we obtain:

$`i_kS\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )=\text{sgn}(\text{Re}(k))f(\overline{k}),`$
$`_{\overline{k}}\stackrel{1}{\mu }(k,\xi ,\eta )+i_{\overline{k}}S\underset{1}{\overset{2}{\mu }}(k,\xi ,\eta )\mathrm{exp}(itS)=`$
$`\text{sgn}(\text{Re}(k))f(k)\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )+`$
$`\text{sgn}(\text{Re}(k))f(k)\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta )\mathrm{exp}(itS),`$
$`i_kS\mathrm{exp}(itS)\underset{1}{\overset{2}{\nu }}(k,\xi ,\eta )+_k\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta )=`$
$`\left(\text{sgn}(\text{Re}(k))f(\overline{k})\stackrel{1}{\mu }(k,\xi ,\eta )_k\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )\right)\mathrm{exp}(itS).`$
If we equate to zero the coefficients of oscillated terms and nonoscillated terms of the expansions (11) and (12) respectively , then we obtain formulas:

$$\underset{1}{\overset{1}{\nu }}=\frac{\text{sgn}(\text{Re}(k))f(\overline{k})}{i_kS};$$
(13)
$$_k\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta )=0;$$
(14)
$$_{\overline{k}}\stackrel{1}{\mu }=\frac{f(\overline{k})f(k)}{i_kS};$$
(15)
$$\underset{1}{\overset{2}{\mu }}(k,\xi ,\eta )=\frac{\text{sgn}(\text{Re}(k))f(k)\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta )}{i_{\overline{k}}S};$$
(16)
$$\stackrel{2}{\nu }=\frac{1}{i_kS}(\text{sgn}(\text{Re}(k))f(\overline{k})\stackrel{1}{\mu }_k\underset{1}{\overset{1}{\nu }}).$$
(17)
The function $`\underset{1}{\overset{1}{\nu }}`$ has a jump on the imaginary axis of $`k`$. This jump is $`\frac{2f(k)}{i_kS}`$. In order for the coefficients of the asymptotics of function $`\nu `$ as $`t^{}1`$ to be continuous, we add an analytic function of $`\overline{k}`$, which has the same jump on the imaginary axis in inverse direction, to $`\underset{1}{\overset{1}{\nu }}\mathrm{exp}(iSt)`$:

$$\underset{0}{\overset{1}{\nu ^{}}}(\overline{k})=\frac{1}{i\pi }_i\mathrm{}^i\mathrm{}\frac{dnf(n)}{(\overline{k}n)_nS}.$$
(18)
The function $`\underset{1}{\overset{1}{\nu ^{}}}`$ define an analytic function $`\underset{0}{\overset{1}{\nu }}`$ of $`\overline{k}`$ uncompletely. The rest terms will be defined when we will match the exterior and interior expansion.

A Cauchy-Green formula gives solution of (15) which is decreasing as $`|k|\mathrm{}`$ and bounded when $`k\text{}`$:

$$\stackrel{1}{\mu }=\frac{1}{2i\pi }_{\text{}}\frac{dpd\overline{p}}{kp}\frac{f(\overline{p})f(p)}{i_pS}.$$
(19)
We can obtain the domain of values of $`k`$, where the external expansion is valid, using the condition $`|t^1\underset{1}{\overset{1}{\nu }}/(t^2\stackrel{2}{\nu })|1`$. As a result of calculations we obtain:

$$\sqrt{t|\theta |}|kk_j|1.$$

Let us construct the interior asymptotic solution which is valid in the neighborhood of point $`k_j`$ as $`j=1,2`$. Denote new scaling variable by $`l_j`$:

$$l_j^2=t(kk_j)^2\frac{_k^2S_j}{2}+\mathrm{\hspace{0.17em}4}t(kk_j)^3.$$
(20)
When $`|l_j|`$ is not so large ($`t^{1/2}|_k^2S_j|^{3/2}|l_j|`$) an asymptotic formula is valid as $`t\mathrm{}`$

$$(kk_j)=\sqrt{\frac{2}{t_k^2S_j}}l_j\frac{8}{t(_k^2S_j)^2}l_j^2+\mathrm{}.$$
(21)
Rewrite the system (9) into terms of new variables $`l_j`$ and $`\overline{l}_j`$. Substitute the asymptotic expansion

$$\mu ^{in}=1+t^1\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta )+\mathrm{},$$
(22)
$$\nu ^{in}=(t^{1/2}\underset{j}{\overset{1}{N}}(l_j,\xi ,\eta )+t^1\underset{j}{\overset{2}{N}}(l_j,\xi ,\eta )+\mathrm{})\mathrm{exp}(itS),$$
(23)
into system (9).

As a result one can obtain equations for the coefficients of expansions (22) and (23):

$$_{l_j}\underset{j}{\overset{1}{N}}2il_j\underset{j}{\overset{1}{N}}=\sqrt{\frac{2}{_k^2S_j}}\text{sgn}(\text{Re}(k_j))f(\overline{k}_j);$$
(24)
$$_{l_j}\underset{j}{\overset{2}{N}}2il_j\underset{j}{\overset{2}{N}}=\text{sgn}(\text{Re}(k_j))[(\frac{16}{(_k^2S_j)^2}f(\overline{k}_j)\frac{2}{_k^2S_j}f_{10}(\overline{k}_j))l_j$$

$$\frac{2\overline{l}_j}{|_k^2S_j|}f_{01}(\overline{k}_j)];$$
(25)
$$_{\overline{l}_j}\underset{j}{\overset{1}{M}}=\text{sgn}(\text{Re}(k_j))f(k_j)\underset{j}{\overset{1}{N}}\sqrt{\frac{2}{_{\overline{k}}^2S_j}}.$$
(26)
Construct the solutions of the equations (24)–(26). Since the external expansion has not the terms of order $`t^{1/2}`$, then the boundary condition for function $`\underset{1}{\overset{1}{N}}(l_j,\xi ,\eta )`$ has the form:

$$\underset{j}{\overset{1}{N}}(l_j,\xi ,\eta )|_{|l_j|\mathrm{}}=0.$$

The solution of the boundary condition for $`\underset{j}{\overset{1}{N}}(l_j,x,y)`$ is evaluated by formula:

$`\underset{j}{\overset{1}{N}}(l_j,\xi ,\eta )=\text{sgn}(\text{Re}(k_j)){\displaystyle \frac{\sqrt{2}f(\overline{k}_j)}{\sqrt{_k^2S_j}}}{\displaystyle \frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2i\pi }}\times `$
$`\times {\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dnd\overline{n}}{\overline{l_jn}}}\mathrm{exp}(i(n^2+\overline{n}^2));`$ (27)
The solution of nonhomogeneous equation for the function $`\underset{j}{\overset{2}{N}}(l_j,\xi ,\eta )`$ has the form:

$`\underset{j}{\overset{s}{\stackrel{2}{N}}}(l_j,\xi ,\eta )=\text{sgn}(\text{Re}(k_j))[{\displaystyle \frac{8f(\overline{k}_j)}{(_k^2S_j)^2}}{\displaystyle \frac{f_{10}(\overline{k}_j)}{_k^2S_j}}`$
$`{\displaystyle \frac{2f_{01}(\overline{k}_j)}{|_k^2S_j|}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2))`$
$`\overline{l}_j{\displaystyle \frac{2f_{01}(\overline{k}_j)}{|_k^2S_j|}}{\displaystyle \frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2i\pi }}{\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dnd\overline{n}}{\overline{l_jn}}}\mathrm{exp}(i(n^2+\overline{n}^2))].`$ (28)
The partial solution of the equation for the function $`\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta )`$ has to be written as four-multiple integral:

$$\underset{j}{\overset{s}{\stackrel{1}{M}}}=\frac{2f(k_j)f(\overline{k}_j)}{|_{\overline{k}}^2S_j|}J,$$
(29)
where

$`J={\displaystyle \frac{1}{2i\pi }}{\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dnd\overline{n}}{l_jn}}{\displaystyle \frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\text{}}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}.`$
It is possible to reduce this four-multiple integral into two-multiple integral (see Appendix):

$`J=\overline{l}_j{\displaystyle \frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2i\pi }}{\displaystyle _{\text{}}\frac{dnd\overline{n}}{\overline{l_jn}}\mathrm{exp}(i(n^2+\overline{n}^2))}`$
$`\mathrm{exp}(i(l_j^2+\overline{l}_j^2)).`$
There exists an domain of values of complex parameter $`k`$, in which the external and internal asymptotic expansions of solution for the problem (9) are valid. In this domain these asymptotic expansions are equal up to the terms of order $`o(t^1)`$. In the domain, where $`t^{1/2}|kk_j|1`$ and $`|kk_j|t^{1/4}`$, the external and internal expansion are valid. We compute asymptotics of external expansion as $`kk_j`$ and internal expansion as $`|l_j|\mathrm{}`$.

Let us present the asymptotics of the functions $`\underset{j}{\overset{1}{N}}`$ and $`\underset{j}{\overset{s}{\stackrel{2}{N}}}`$ as $`|l_j|\mathrm{}`$:

$`\underset{j}{\overset{1}{N}}|_{|l_j|\mathrm{}}=\text{sgn}(\text{Re}(k_j))\sqrt{{\displaystyle \frac{2}{_k^2S_j}}}f(\overline{k}_j)\times `$
$`\left({\displaystyle \frac{1}{2il_j}}+{\displaystyle \frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{\overline{l}_j}}+O(|l_j|^3)\right).`$
$`\underset{j}{\overset{2}{N}}(l_j,\xi ,\eta )|_{|l_j|\mathrm{}}=\text{sgn}(\text{Re}(k_j))[({\displaystyle \frac{8if(\overline{k}_j)}{(_k^2S_j)^2}}{\displaystyle \frac{if_{10}(\overline{k}_j)}{_k^2S_j}})`$
$`+{\displaystyle \frac{\overline{l}_j}{2il_j}}{\displaystyle \frac{2f_{01}(\overline{k}_j)}{|_{\overline{k}}^2S_j|}}+O(|l_j|^1)].`$
The matching condition for $`\stackrel{~}{\nu }`$ means that in the domain $`t^{1/3}|kk_j|1`$ and $`|kk_j|=o(1)`$:

$$\left(t^{1/2}\underset{j}{\overset{1}{N}}(l_j,\xi ,\eta )+t^1\underset{j}{\overset{2}{N}}(l_j,\xi ,\eta )\right)\mathrm{exp}(itS)$$

$$t^1(\underset{1}{\overset{1}{\nu }}(k,\xi ,\eta )\mathrm{exp}(itS)+\underset{0}{\overset{1}{\nu }}(k,\xi ,\eta ))=o(t^1).$$

Using the asymptotics of interior expansions of $`\underset{j}{\overset{2}{N}}`$, which are rewritten in the terms of external variable $`k`$ as $`kk_{1,2}`$, and external expansion of $`\underset{0}{\overset{1}{\nu }}`$, we can obtain:

$`\underset{0}{\overset{1}{\nu }}(k,x,y)`$ $`=\underset{0}{\overset{1}{\nu ^{}}}(\overline{k})+{\displaystyle \frac{2f(\overline{k}_1)\mathrm{exp}(itS_1)}{|_k^2S_1|\overline{(kk_1)}}}+`$ (30)
$`+{\displaystyle \frac{2f(\overline{k}_2)\mathrm{exp}(itS_2)}{|_k^2S_2|\overline{(kk_2)}}},`$
where the function $`\underset{0}{\overset{}{\stackrel{1}{\nu }}}(\overline{k})`$ is defined by (18);

$`\underset{j}{\overset{2}{N}}=\underset{j}{\overset{s}{\stackrel{2}{N}}}+\underset{0}{\overset{1}{\nu ^{}}}(k_j)\mathrm{exp}(i(l_j^2+\overline{l}_j^2))\mathrm{exp}(itS_j)+`$
$`+{\displaystyle \frac{2\text{sgn}(\text{Re}(k_m))f(\overline{k}_m)\mathrm{exp}(itS_m)}{|_k^2S_m|\overline{(k_jk_m)}}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2)),`$ (31)
where $`mj`$, the function $`\underset{j}{\overset{s}{\stackrel{2}{N}}}`$ is defined by (28).

Let us construct $`\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta )`$. The asymptotics of (29) as $`|l_j|\mathrm{}`$ has the form:

$$\underset{j}{\overset{s}{\stackrel{1}{M}}}|_{|l_j|\mathrm{}}=\frac{\overline{l}_j}{l_j}\sqrt{\frac{_k^2S_j}{_{\overline{k}}^2S_j}}\frac{i}{12\pi }\frac{f(\overline{k}_j)f(k_j)}{_k^2S}.$$

Evaluate the asymptotics of $`\underset{0}{\overset{1}{\mu }}`$ as $`kk_j`$:

$$\underset{0}{\overset{1}{\mu }}|_{kk_j}=\frac{\overline{(kk_j)}}{kk_j}\frac{f(\overline{k}_j)f(k_j)}{12i(k_jk_n)}\frac{f(k_j)f(\overline{k}_j)}{12i(k_jk_n)}+$$

$$+\frac{1}{2i\pi }_{\text{}}𝑑pd\overline{p}\frac{f(\overline{p})f(p)f(k_j)f(\overline{k}_j)}{12i(pk_j)^2(pk_n)},$$

here $`n1,2,nj`$.

The matching condition for $`\stackrel{~}{\mu }`$ when $`t^{1/3}|kk_j|1`$ and $`|kk_j|=o(t^{1/4})`$ has the form:

$$(1+t^1\underset{0}{\overset{1}{\mu }}(k,\xi ,\eta ))(1+t^1\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta ))=o(t^1).$$

These condition allows to define $`\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta )`$:

$$\underset{j}{\overset{1}{M}}(l_j,\xi ,\eta )=\underset{j}{\overset{s}{\stackrel{1}{M}}}(l_j,\xi ,\eta )+C_j(\xi ,\eta ).$$
(32)
Here the function $`\underset{j}{\overset{s}{\stackrel{1}{M}}}(l_j,\xi ,\eta )`$ is defined by (29), $`C_j(\xi ,\eta )`$ has the form:

$`C_j(\xi ,\eta )={\displaystyle \frac{f(k_j)f(\overline{k}_j)}{12i(k_jk_n)}}+`$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle _{\text{}}𝑑p}d\overline{p}{\displaystyle \frac{f(\overline{p})f(p)f(k_j)f(\overline{k}_j)}{12i(pk_j)^2(pk_n)}},`$ (33)
where $`n1,2,nj`$.

Thus we have constructed the interior expansion in neighborhood of nondegenerate stationary point of $`S`$ as $`|12\xi +\eta ^2|t^{1/4}`$.

Lemma is proved.

Constructed asymptotic solutions are nonuniform with respect to $`k`$. But one can obtain an uniform asymptotic solution by using their combination. This uniform solution has the form (see, for example, ):

$`\left(\begin{array}{c}\widehat{\mu }\\ \widehat{\nu }\end{array}\right)=\left(\begin{array}{c}\stackrel{~}{\mu }\\ \stackrel{~}{\nu }\end{array}\right)+\left(\begin{array}{c}\stackrel{~}{M}_1\\ \stackrel{~}{N}_1\end{array}\right)+\left(\begin{array}{c}\stackrel{~}{M}_2\\ \stackrel{~}{N}_2\end{array}\right)`$ (42)
$`A_{1,k}\left(\begin{array}{c}\stackrel{~}{M}_1\\ \stackrel{~}{N}_1\end{array}\right)A_{1,k}\left(\begin{array}{c}\stackrel{~}{M}_2\\ \stackrel{~}{N}_2\end{array}\right).`$ (47)
The result of action of operator $`A_{n,k}`$ on the function $`\stackrel{~}{M}`$ is defined as follows (). Take the formula for $`\stackrel{~}{M}(l_j,\xi ,\eta ,t)`$, and change the dependence on variable $`l_j`$ into the dependence on variable $`k`$ using the formula (20). Rewrite the sum of all terms of asymptotic expansion up to $`t`$ with powers equal to $`m`$, where $`0mn`$. For example for the functions $`\stackrel{~}{M}(l_1,\xi ,\eta ,t)`$ and $`\stackrel{~}{N}(l_1,\xi ,\eta ,t)`$ this process leads to the formulas:

$$A_{1,k}(\stackrel{~}{M}_1)=1+t^1\frac{\overline{\sqrt{\theta (kk_1)^2+4(kk_1)^3}}}{\sqrt{\theta (kk_1)^2+4(kk_1)^3}}\frac{f(\overline{k}_1)f(k_1)}{2|_k^2S_1|};$$

$`A_{1,k}(\stackrel{~}{N}_1)=t^1[{\displaystyle \frac{f(\overline{k}_1)}{2i\theta \sqrt{\theta (kk_1)^2+4(kk_1)^3}}}`$
$`{\displaystyle \frac{f(\overline{k}_1)\mathrm{exp}(it(SS_1))}{|\theta |\overline{\sqrt{\theta (kk_1)^2+4(kk_1)^3}}}}+`$
$`+\underset{0}{\overset{1}{\nu ^{}}}{\displaystyle \frac{f_{01}(\overline{k}_1)}{|\theta |}}{\displaystyle \frac{\overline{\sqrt{\theta (kk_1)^2+4(kk_1)^3}}}{\sqrt{\theta (kk_1)^2+4(kk_1)^3}}}].`$
If we substitute (47) into (9) and evaluate a remainder, then we obtain:

###### Theorem 3

The formal asymptotic solution of the problem (9) with respect to $`mod\left(O\left((t|\theta |)^1\right)\right)`$, which is uniformly valuable when $`k\text{}`$ and $`\theta ^2t1`$, has the form (47).

#### 4.1.2 The stationary points on the imaginary axis

If $`\text{Re}(\theta )=0`$, then the stationary point of the phase function $`S`$ belongs to the line, where the coefficients of the equation (9) are discontinuous. In this case constructing of the formal asymptotic solution of problem (9) differs from the asymptotic solution which was constructed before. The main result about the combined asymptotic solution is:

###### Lemma 3

Let the system of equations (9) have no homogeneous solutions, $`F(k)C^2L_1`$ as $`\text{Re}(k)0`$; the parameters $`\xi `$ and $`\eta `$ are $`t^{2/3}(12\xi +\eta ^2)1`$, then:
when $`\sqrt{t|\theta |}|kk_{1,2}|1`$the formal asymptotic solution of (9) with respect to $`mod(O(t^2|_kS|^3))`$ has the form:

$$\stackrel{~}{\mu }=1+t^1\stackrel{1}{\alpha }(k,\xi ,\eta ),$$

$$\stackrel{~}{\nu }=(t^1\underset{1}{\overset{1}{\beta }}(k,\xi ,\eta )+t^2\underset{1}{\overset{2}{\beta }}(k,\xi ,\eta ))\mathrm{exp}(itS)+t^1\underset{0}{\overset{1}{\beta }}(k,\xi ,\eta );$$

the functions $`\stackrel{1}{\alpha },\underset{1}{\overset{1}{\beta }},\underset{1}{\overset{2}{\beta }},\underset{0}{\overset{1}{\beta }}`$ are defined by (50), (51), (54), (60);
when $`|\theta ||kk_j|1`$ the formal asymptotic solution of (9) with respect to $`mod(O(t^1|\theta |^1))`$ has the form:

$$\stackrel{~}{Y}=1+t^1\stackrel{1}{Y}(l_j,\xi ,\eta ),$$

$$\stackrel{~}{Z}=\left(t^{1/2}\stackrel{1}{Z}(l_j,\xi ,\eta )+t^1\stackrel{2}{Z}(l_j,\xi ,\eta )\right)\mathrm{exp}(itS),$$

where the functions $`\stackrel{1}{Y}(l_j,\xi ,\eta ),\stackrel{1}{Z}(l_j,\xi ,\eta ),\stackrel{2}{Z}(l_j,\xi ,\eta )`$ are defined by (62), (61).

The proof. Let us construct the asymptotics. The external expansion is constructed similarly as in 4.1.1. The main difference is that the stationary points of the function $`S`$ are on the line of discontinuity of the coefficients of system (9). It leads to sufficient modifications in the formulas for asymptotics. Let us find the external expansion in the form:

$$\mu ^{ex}=1+t^1\stackrel{1}{\alpha }(k,\xi ,\eta )+t^2\underset{1}{\overset{2}{\alpha }}(t,\xi ,\eta )\mathrm{exp}(itS)+\mathrm{},$$
(48)
$$\nu ^{ex}=t^1(\stackrel{1}{\beta }(k,\xi ,\eta )\mathrm{exp}(itS)+\underset{0}{\overset{1}{\beta }}(k,\xi ,\eta ))+t^2\stackrel{2}{\beta }\mathrm{exp}(itS)+\mathrm{}.$$
(49)
After substituting of (48) and (49) into the system (9) we equate the coefficients with identical powers of $`t`$ and with oscillated and nonoscillated terms correspondingly. As a result we obtain equations for the coefficients of asymptotic expansions (48) and (49). The obvious formulas are

$$\stackrel{1}{\alpha }=\frac{1}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{kn}\frac{f(\overline{n})f(n)}{i_nS(n,\xi ,\eta )},$$
(50)
$$\underset{1}{\overset{1}{\beta }}=\frac{\text{sgn}(\text{Re}(k))f(\overline{k})}{i_kS}.$$
(51)
The function $`\underset{1}{\overset{1}{\beta }}`$ is discontinuous on the imaginary axis. We add an analytic function with the same jump on the imaginary axis in back direction then the coefficient of the expansion (49) at $`t^1`$ is continuous:

$$\underset{0}{\overset{1}{\beta ^{}}}=\frac{1}{\pi i_{\overline{k}}S}_i\mathrm{}^i\mathrm{}\frac{dnf(n)_{\overline{n}}S}{(\overline{k}n)_nS}=\frac{1}{\pi i_{\overline{k}}S}_i\mathrm{}^i\mathrm{}\frac{dnf(n)}{\overline{k}n}.$$
(52)
Unlike $`\underset{0}{\overset{1}{\nu ^{}}}`$, this function have singularities on the discontinuous line.

The formulas for the other coefficients of the expansions (48) and (49) have the forms:

$$\underset{1}{\overset{2}{\alpha }}=\frac{\text{sgn}(\text{Re}(k))f(k)\underset{0}{\overset{1}{\beta }}(k,\xi ,\eta )}{i_{\overline{k}}S},$$
(53)
$$\stackrel{2}{\beta }=\frac{1}{i_kS}(\text{sgn}(\text{Re}(k)f(\overline{k})\stackrel{1}{\alpha }_k\underset{1}{\overset{1}{\beta }}).$$
(54)
Here the analytic function $`\underset{0}{\overset{1}{\beta }}`$ is still undefined. We will finally define this function after matching of external and internal expansions.

The interior expansion in the neighborhood of the point $`k_j`$ depend on the scaling variable $`l_j`$. The asymptotics of the functions $`\mu `$ and $`\nu `$ has the same asymptotic sequence, but the equations for $`\underset{j}{\overset{1}{Z}}`$, $`\underset{j}{\overset{2}{Z}}`$ and $`\underset{j}{\overset{1}{Y}}`$ have discontinuous right-hand sides.

$`_{l_j}\underset{j}{\overset{1}{Z}}2il_j\underset{j}{\overset{1}{Z}}=\sqrt{{\displaystyle \frac{2}{_k^2S_j}}}f(\overline{k}_j)\text{sgn}(\text{Re}(\overline{l}\mathrm{exp}(i\pi /4)));`$ (55)
$`_{l_j}\underset{j}{\overset{2}{Z}}2il_j\underset{j}{\overset{2}{Z}}=\text{sgn}(\text{Re}[\overline{l}\mathrm{exp}(i\pi /4)])\times `$
$`\times \left({\displaystyle \frac{16}{(_k^2S_j)^2}}f_{00}(\overline{k}_j)+{\displaystyle \frac{2l_j}{_k^2S_j}}f_{01}^{(1)}(\overline{k}_j)+{\displaystyle \frac{2i\overline{l}_j}{|_k^2S_j|}}f_{10}^{(1)}(\overline{k}_j)\right)+`$
$`+{\displaystyle \frac{2l_j}{_k^2S_j}}f_{01}^{(2)}(\overline{k}_j)+{\displaystyle \frac{2i\overline{l}_j}{|_k^2S_j|}}f_{10}^{(2)}(\overline{k}_j).`$ (56)
$`_{\overline{l}_j}\underset{j}{\overset{1}{Y}}=\text{sgn}(\text{Re}(l\mathrm{exp}(i\pi /4)))f(k_j)\underset{j}{\overset{1}{Z}}\sqrt{{\displaystyle \frac{2}{_{\overline{k}}^2S_j}}}.`$ (57)
Continuous partial solutions of the equations (55), (56) and (57) are obtained by Cauchy-Green formula. The formulas for the partial solutions of the equations (55) and (56) are differing from the formulas obtained in section 4.1.1 by function $`\text{sgn}(\text{Re}(\overline{n}\mathrm{exp}(i\pi )))`$ under the sign of double integral:

$$\underset{j}{\overset{s}{\stackrel{1}{Z}}}=𝒥[\stackrel{1}{H_j}](l_j),\underset{j}{\overset{s}{\stackrel{2}{Z}}}=𝒥[\stackrel{2}{H_j}](l_j),$$
(58)
where $`\underset{j}{\overset{1}{H}}`$ and $`\underset{j}{\overset{2}{H}}`$ are the right-hand side of the equations (55) and (56). An operator $`𝒥`$ has the form:

$$𝒥[h](l_j)=\frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{\overline{l_jn}}h(n)\mathrm{exp}(i(l_j^2+\overline{l}_j^2)).$$

The partial solution of (57) contains of the double integral the same as in (29) and the term, which is defined by integral over the discontinuous line. The function $`\underset{j}{\overset{s}{\stackrel{1}{Y}}}`$ has the form:

$$\underset{j}{\overset{s}{\stackrel{1}{Y}}}=\frac{2f(k_j)f(\overline{k}_j)}{|_{\overline{k}}^2S_j|}(J2J_+2J_+),$$
(59)
where:

$`J_\pm ={\displaystyle }{\displaystyle _\mathrm{\Omega }^{}}{\displaystyle \frac{dmd\overline{m}}{lm}}{\displaystyle \frac{\mathrm{exp}(i(m^2+\overline{m}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^\pm }\frac{dnd\overline{n}}{\overline{nm}}\mathrm{exp}(i(n^2+\overline{n}^2))}.`$
Here $`\mathrm{\Omega }^\pm =\{\pm \text{Re}[l\mathrm{exp}(i\pi /4)>0\}`$. Note that $`J_+(l)=J_+(l)`$. In the Appendix we show, that we can reduce the integral $`J_+`$ into double integral. As a result we obtain:

$`J_+={\displaystyle \frac{3}{4}}i\pi +`$
$`il{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{il\overline{m}}\mathrm{exp}(i(m^2+\overline{m}^2))},\text{when}l\mathrm{\Omega }^+;`$
$`J_+=\overline{l}\mathrm{exp}(i(l^2+\overline{l}^2)){\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{il\overline{m}}\mathrm{exp}(i(m^2+\overline{m}^2))}`$
$`{\displaystyle \frac{1}{4}}i\pi {\displaystyle \frac{1}{2}}i\pi \mathrm{exp}(i(l^2+\overline{l}^2)),\text{when}l\mathrm{\Omega }^{};`$
Thus, the function $`\underset{j}{\overset{s}{\stackrel{1}{Y}}}`$ is represented by sum of double integrals.

Let us consider matching of the external and internal asymptotic expansions. For that we need the asymptotic behaviour of the partial solutions of (55) and (56) as $`|l|\mathrm{}`$. The asymptotic behaviour of the double integrals is evaluated in the Appendix. By using it we obtain:

$`\underset{j}{\overset{s}{\stackrel{1}{Z}}}|_{|l|\mathrm{}}=\text{sgn}\left[\text{Re}[\overline{l}_j\mathrm{exp}(i\pi /4)]\right]f(\overline{k}_j)\sqrt{{\displaystyle \frac{2}{_k^2S_j}}}({\displaystyle \frac{1}{2il_j}}+`$
$`+{\displaystyle \frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2\overline{l}_j}}+O(|l_j|^2));`$
$`\underset{j}{\overset{s}{\stackrel{2}{Z}}}|_{|l_j|\mathrm{}}=\text{sgn}\left[\text{Re}[\overline{l}_j\mathrm{exp}(i\pi /4)]\right]\times `$
$`[({\displaystyle \frac{\overline{l}_j}{2il_j}}{\displaystyle \frac{1}{2}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2))){\displaystyle \frac{2if_{10}^{(1)}}{|_k^2S_j|^2}}+`$
$`+({\displaystyle \frac{1}{4}}{\displaystyle \frac{1}{4}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2)))2f_{01}^{(1)}]+`$
$`{\displaystyle \frac{2if_{10}^{(2)}}{|_k^2S_j|}}\left({\displaystyle \frac{\overline{l}_j}{2il}}+{\displaystyle \frac{1}{2}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2))\right)+if_{01}^{(2)}+O(|l_j|^1).`$
The asymptotic behaviour of the function $`\underset{0}{\overset{\prime \prime }{\stackrel{1}{\beta }}}`$ is:

$`\underset{0}{\overset{\prime \prime }{\stackrel{1}{\beta }}}|_{kk_j}={\displaystyle \frac{1}{12i\overline{(kk_j)(k_jk_m)}}}[f(k_j)\text{sgn}[\text{Re}(\overline{k})]+`$
$`+{\displaystyle \frac{1}{\pi i}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dn{\displaystyle \frac{f(in)f(k_j)}{n\overline{ik_j}}}]+`$
$`{\displaystyle \frac{1}{12i(k_jk_m)^2}}[f(k_j)\text{sgn}[\text{Re}(\overline{k})]+`$
$`{\displaystyle \frac{1}{\pi i}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dn{\displaystyle \frac{f(in)f(k_j)}{n\overline{ik_j}}}]+{\displaystyle \frac{1}{12i\overline{(k_jk_m)}}}[f_{01}^{(1)}(k_j)+`$
$`+f_{01}^{(2)}(k_j)\text{sgn}[\text{Im}(\overline{ik})]+f_{10}^{(1)}(k_j)+f_{10}^{(2)}(k_j)\text{sgn}[\text{Im}(\overline{ik})]+`$
$`+{\displaystyle \frac{1}{\pi i}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dn{\displaystyle \frac{f(in)f(k_j)\overline{(inik_j)}f_{01}^{(2)}(k_j)(inik_j)f_{10}^{(2)}(k_j)}{(n\overline{ik_j})^2}}]`$
$`+o(1).`$
Using the matching conditions for the asymptotics of the external $`\nu ^{ex}`$ and internal $`\nu ^{in}`$ expansions we obtain:

$$\underset{0}{\overset{1}{\beta }}=\underset{0}{\overset{1}{\beta ^{}}}\frac{C_1^1}{12i\overline{(kk_1)(k_1k_2)}}\frac{C_2^1}{12i\overline{(kk_2)(k_2k_1)}},$$
(60)
where $`\underset{0}{\overset{}{\stackrel{1}{\beta }}}`$ is defined by formula (52),

$$C_j^1=\frac{1}{\pi i}_{\mathrm{}}^{\mathrm{}}𝑑n\frac{f(n)f(k_j)}{n\overline{ik_j}},j=1,2;$$

$$\underset{j}{\overset{1}{Z}}=\underset{j}{\overset{s}{\stackrel{1}{Z}}}\underset{j}{\overset{2}{Z}}=\underset{j}{\overset{s}{\stackrel{2}{Z}}}\mathrm{exp}(i(l_j^2+\overline{l}_j^2))C_j^2,$$
(61)
where functions $`\underset{j}{\overset{s}{\stackrel{1}{Z}}}`$ and $`\underset{j}{\overset{s}{\stackrel{2}{Z}}}`$ are defined by (58),

$`C_j^1={\displaystyle \frac{i}{_k^2S_j}}f_{10}^{(2)}(\overline{k}_j)`$
$`{\displaystyle \frac{1}{\pi i}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑n{\displaystyle \frac{f(n)f(\overline{k}_j)f_{01}^{(2)}(\overline{k}_j)\overline{(nik_j)}(nik_j)f_{10}^{(2)}(k_j)}{(n\overline{ik_j})^2}}.`$
The asymptotic behaviour of the function $`\underset{j}{\overset{s}{\stackrel{1}{Y}}}`$ is :

$$\underset{j}{\overset{s}{\stackrel{1}{Y}}}=(i\pi +\frac{i\overline{l}_j}{12\pi l_j})\frac{f(\overline{k}_j)f(k_j)}{|_k^2S_j|}.$$

The matching condition for the external expansion of $`\mu ^{ex}`$ and for the internal expansion of $`\mu ^{in}`$ gives:

$$\underset{j}{\overset{1}{Y}}=\underset{j}{\overset{s}{\stackrel{1}{Y}}}+C_j(\xi ,\eta )+i\pi \frac{f(\overline{k}_j)f(k_j)}{|_k^2S_j|},$$
(62)
where $`\underset{j}{\overset{s}{\stackrel{1}{Y}}}`$ is defined by (59), The function $`C_j(\xi ,\eta )`$ is defined by (33).

The lemma is proved.

Constructed external and internal expansions are irregular with respect to parameters $`k,\xi ,\eta `$. The uniform expansion with respect to $`k`$ when$`(12\xi +\eta ^2)t^{1/3}1`$ has to be constructed by the same way as in preceding section:

$$\left(\begin{array}{c}\widehat{\mu }\\ \widehat{\nu }\end{array}\right)=\left(\begin{array}{c}\stackrel{~}{\mu }\\ \stackrel{~}{\nu }\end{array}\right)+\left(\begin{array}{c}\stackrel{~}{Y}_1\\ \stackrel{~}{Z}_1\end{array}\right)+\left(\begin{array}{c}\stackrel{~}{Y}_2\\ \stackrel{~}{Z}_2\end{array}\right)A_{1,k}\left(\begin{array}{c}\stackrel{~}{Y}_1\\ \stackrel{~}{Z}_1\end{array}\right)A_{1,k}\left(\begin{array}{c}\stackrel{~}{Y}_2\\ \stackrel{~}{Z}_2\end{array}\right).$$
(63)
The operator $`A_{n,k}`$ was defined above. Let us substitute (63) into (9) and evaluate a remainder. As a result we obtain

###### Theorem 4

The formal asymptotic solution of the problem (9) with respect to $`mod\left(O\left((t|\theta |)^1\right)\right)`$, which is uniform when $`k\text{}`$ and $`\theta ^2t^{2/3}1`$, has the form (63).

### 4.2 Asymptotics in neighborhood of confluent stationary point

The system of equations in the problem (9) depends on two control parameters $`\xi `$ and $`\eta `$. On the parabola $`12\xi +\eta ^2=0`$ the degeneracy of stationary points occurs: $`k_1=k_2=k_0=i\eta `$. For this reason the asymptotics constructed in section 4.1 is invalid when parameter $`\theta =\sqrt{\eta ^2+12\xi }`$ is close to zero. For example, the asymptotics of $`\stackrel{1}{\mu }`$ as $`kk_0`$ and $`\theta 0`$ is discontinuous:

$`\left[\stackrel{1}{\mu }|_{kk_0}\right]|_{\theta 0}=({\displaystyle \frac{|kk_0|}{12i((kk_0)^2\theta ^2)}}{\displaystyle \frac{\overline{\theta }}{\theta }}{\displaystyle \frac{kk_0}{12i((kk_0)^2\theta ^2)}})\times `$
$`f(\overline{k}_0)f(k_0)+o(1).`$
This shows that we need new scaling of the parameter $`\theta `$. The scaled control parameter is:

$$v=t^{1/3}\frac{\theta }{\sqrt{12}}.$$

The external expansion constructed in sec. 4.1 becomes discontinuous at $`\theta =0`$ and the internal expansions constructed in sec. 4.1 become singular at the point $`\theta =0`$ and lose their asymptotic properties. Therefore here we must change the internal variable for the internal asymptotic expansion:

$$p=t^{1/3}(kk_0).$$
(64)
In this section we construct a formal asymptotic solution of the problem (9) with respect to $`mod(O(t^1))`$ when $`|\theta |1`$ uniform with respect to $`k\text{}`$. The result is formulated in the end of this section.

To construct the uniform asymptotic solution we need the external and internal asymptotics outside and inside a small neighborhood of the point $`k_0`$ respectively.

###### Lemma 4

Let the system of the equations (9)be have no the homogeneous nontrivial solutions, $`F(k)C^2L_1`$ and the parameters $`\xi `$ and $`\eta `$ satisfy the inequality $`|12\xi +\eta ^2|1`$, then:
when $`|kk_0|t^{1/3}1`$ the formal asymptotic solution of the system (9) with respect to $`mod(O(t^{2/3}/|kk_0|)+O(t^1))`$ has the form:

$$\stackrel{~}{\mu }=1+t^1\stackrel{1}{m}(k,\xi ,\eta ),$$

$$\stackrel{~}{\nu }=t^{2/3}\underset{0}{\overset{1}{n}}+t^1(\underset{1}{\overset{1}{n}}\mathrm{exp}(itS)+\underset{0}{\overset{2}{n}})$$

the function $`\stackrel{1}{m}`$ is defined by (90), the functions $`\underset{0}{\overset{1}{n}}`$ and $`\underset{2}{\overset{2}{n}}`$ are defined by (86) and (87) respectively, the function $`\underset{1}{\overset{1}{n}}`$ is defined by (77);
when $`|kk_0|1`$ the asymptotic solution of the system (9)with respect to $`mod(O(t^{2/3}|kk_0|+O(t^1))`$ has the form:

$$\stackrel{~}{\mu }=1+t^{2/3}\stackrel{1}{}+t^1\stackrel{2}{},$$

$$\stackrel{~}{\nu }=(t^{1/3}\stackrel{1}{𝒩}+t^{2/3}\stackrel{2}{𝒩}+t^1\stackrel{3}{𝒩})\mathrm{exp}(itS);$$

the functions $`\stackrel{j}{}`$, $`j=1,2`$ and $`\stackrel{j}{𝒩}`$, $`j=1,2,3`$ are defined by (89), and (88) respectively.

Proof. Let us find the internal formal asymptotic expansion for the solution of the system (9) in the form:

$$^{in}=1+t^{2/3}\stackrel{1}{}+t^1\stackrel{2}{}+\mathrm{},$$
(65)
$$𝒩^{in}=(t^{1/3}\stackrel{1}{𝒩}+t^{2/3}\stackrel{2}{𝒩}+t^1\stackrel{3}{𝒩}+\mathrm{})\mathrm{exp}(itS).$$
(66)
Change the variable $`k`$ into the variable $`p`$ in the system (9). The phase function of the exponent is depend on the variable $`p`$ as:

$$tS\omega (p)4(p^3+\overline{p}^3)v^2(p+\overline{p}),$$

where $`v=t^{1/3}\theta /\sqrt{12}`$. Substitute (65) and (66) into the system (9), equate the coefficients with equal power of $`t`$. As a result we obtain a sequence of equation for the coefficients of the expansions (65) and (66).

$`_p\stackrel{1}{𝒩}i(12p^2v^2)\stackrel{1}{𝒩}=\text{sgn}[\text{Re}(\overline{p})]f(\overline{k}_0),`$ (67)
$`_p\stackrel{2}{𝒩}i(12p^2v^2)\stackrel{2}{𝒩}=\text{sgn}[\text{Re}(\overline{p})](f_{01}^{(1)}(\overline{k}_0)\overline{p}+`$
$`+f_{10}^{(1)}(\overline{k})p)f^{(2)}_{01}(\overline{k}_0)\overline{p}f^{(2)}_{10}(\overline{k}_0)p.`$ (68)
$`_{\overline{p}}\stackrel{1}{}=\text{sgn}[\text{Re}(\overline{p})]f(k_0)\stackrel{1}{𝒩},`$ (69)
$`_p\stackrel{3}{𝒩}i(12p^2v^2)\stackrel{3}{𝒩}=`$
$`{\displaystyle \frac{1}{2}}\left[\text{sgn}\left[\text{Re}(\overline{p})\right]f_{20}^{(1)}(\overline{k}_0)+f_{20}^{(2)}(\overline{k}_0)\right]p^2`$
$`{\displaystyle \frac{1}{2}}\left[\text{sgn}\left[\text{Re}(\overline{p})\right]f_{02}^{(1)}(\overline{k}_0)+f_{02}^{(2)}(\overline{k}_0)\right]\overline{p}^2`$
$`\left[\text{sgn}\left[\text{Re}(\overline{p})\right]f_{11}^{(1)}(\overline{k}_0)+f_{11}^{(2)}(\overline{k}_0)\right]|p|^2`$
$`\text{sgn}\left[\text{Re}(\overline{p})\right]f(\overline{k}_0)\stackrel{1}{}.`$ (70)
$`_{\overline{p}}\stackrel{2}{}=\text{sgn}\left[\text{Re}(p)\right]f(k_0)\stackrel{2}{𝒩}+`$
$`\left[\text{sgn}\left[\text{Re}(\overline{p})\right]f_{10}^{(1)}(k_0)+f_{10}^{(2)}(k_0)\right]\stackrel{1}{𝒩}p+`$
$`+\left[\text{sgn}[\text{Re}(k)]f_{01}^{(1)}(k_0)+f_{01}^{(2)}(k_0)\right]\stackrel{1}{𝒩}\overline{p}.`$ (71)
These equations are obtained in a supposition that $`|p|t^{1/3}1`$. The uniform bounded partial solutions of the equations for the coefficients of the expansion (66) are obtained by using an integral operator:

$$𝒫[g]=\frac{\mathrm{exp}(i\omega (p))}{2i\pi }_{\text{}}\frac{drd\overline{r}}{\overline{pr}}\mathrm{exp}(i\omega (r))g(r),$$
(72)
where $`g(r)`$ is the right-hand side of corresponding equation.

Thus, the bounded partial solutions of the equations (67) and (68) (the functions $`\stackrel{s}{\stackrel{1}{𝒩}}`$ and $`\stackrel{s}{\stackrel{2}{𝒩}}`$) are represented by double integrals.

Bounded partial solutions of equations for the coefficients of the expansion (65) are constructed by using the Cauchy-Green formula.

The formula for the partial solution of the equation (69) has the form:

$$\stackrel{s}{\stackrel{1}{}}(p,\xi ,\eta )=f(k_0)f(\overline{k}_0)\left(J_1(p,v^2)2J_1^+2J_1^+\right).$$
(73)
Here

$$J_1=\frac{1}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{pn}\frac{\mathrm{exp}(i\omega (n))}{2i\pi }_{\text{}}\frac{drd\overline{r}}{\overline{nr}}\mathrm{exp}(i\omega (r)).$$

$$J_1^\pm =\frac{1}{2i\pi }_\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\frac{\mathrm{exp}(i\omega (n))}{2i\pi }_{\mathrm{\Omega }^\pm }\frac{drd\overline{r}}{\overline{nr}}\mathrm{exp}(i\omega (r)).$$

In the Appendix the integrals $`J_1`$ and $`J_1^+`$ are reduced into double integrals. By the similar way we may represent the integral $`J_1^+`$ as the sum of double integrals. Therefore we may represent $`\stackrel{s}{\stackrel{1}{}}`$ as sum of double integrals also. Corresponded formula is very large and we don’t write it here. But we will use this formula to evaluate an asymptotic behaviour of $`\stackrel{s}{\stackrel{1}{}}`$ as $`|p|\mathrm{}`$.

The bounded partial solution of the equation (71) when $`p\text{}`$ has to be built similarly . We give following statement about $`\stackrel{2}{}`$.

###### Lemma 5

Continuous partial solution of the equation for $`\stackrel{2}{}`$ exists. This solution is uniformly bounded when $`p\text{}`$.

Sketch of proof. The partial solution of the equation for $`\stackrel{2}{}`$ has to be obtained as a result of using of the operator:

$$_{\overline{p}}^1\left[g\right]=\stackrel{s}{\stackrel{2}{}}\frac{1}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{pn}g(n)$$
(74)
to the right-hand side of the equation. The integrals, which are obtained, are continuous by virtue of continuity of the integral operator with respect to the parameter $`p`$. The boundedness of these integrals with respect to $`p`$ when $`p\text{}`$ may be obtained by using the asymptotic behaviour of the right-hand side terms of the equation for $`\stackrel{2}{}`$ as $`|p|\mathrm{}`$.

Let us to construct the external expansion in the form:

$$m=1+t^1\underset{0}{\overset{1}{m}}(k,v)+t^{5/3}(\underset{1}{\overset{2}{m}}\mathrm{exp}(itS)+\underset{0}{\overset{2}{m}})+\mathrm{};$$
(75)
$`n`$ $`=(t^1\underset{1}{\overset{1}{n}}(k,v)+t^{5/3}\underset{1}{\overset{2}{n}}(k,v)+t^2\underset{1}{\overset{3}{n}}\mathrm{})\mathrm{exp}(itS)`$ (76)
$`+t^{2/3}\underset{0}{\overset{1}{n}}(k,v)+t^1\underset{0}{\overset{2}{n}}(k,v)+\mathrm{}.`$
Substitute (75) and (76) into (9). As a result we get the equations for coefficients of the expansions:

$$\underset{1}{\overset{1}{n}}=\text{sgn}\left[\text{Re}(\overline{k})\right]\frac{f(\overline{k})}{12i(kk_0)^2};\underset{1}{\overset{2}{n}}=\text{sgn}\left[\text{Re}(\overline{k})\right]\frac{v^2f(\overline{k})}{144i(kk_0)^4};$$
(77)
$$_{\overline{k}}\stackrel{1}{m}=\frac{f(k)f(\overline{k})}{12i(kk_0)^2};$$
(78)
$$_k\underset{0}{\overset{1}{n}}(k,\xi ,\eta )=0;_k\underset{0}{\overset{2}{n}}(k,\xi ,\eta )=0;$$
(79)
$`_{\overline{k}}\underset{0}{\overset{2}{m}}={\displaystyle \frac{f(k)f(\overline{k})v^2}{144i(kk_0)^4}};`$ (80)
$`\underset{1}{\overset{2}{m}}(k,\xi ,\eta )=\text{sgn}\left[\text{Re}(\overline{k})\right]f(k){\displaystyle \frac{\underset{0}{\overset{1}{n}}(k,\xi ,\eta )}{12i\overline{(}kk_))^2}};`$ (81)
$`\underset{1}{\overset{3}{n}}(k,\xi ,\eta )={\displaystyle \frac{1}{12i(kk_0)^2}}(\text{sgn}\left[\text{Re}(\overline{k})\right]f(\overline{k})\underset{0}{\overset{1}{m}}`$
$`+{\displaystyle \frac{2f(\overline{k})}{12i(kk_0)^3}}\text{sgn}[\text{Re}(\overline{k})]{\displaystyle \frac{f_{10}^{(1)}(\overline{k})\text{sgn}[\text{Re}(\overline{k})]+f_{10}^{(2)}(\overline{k})}{12i(kk_0)^2}}).`$ (82)
The formulas (77) and (82) define the functions $`\underset{1}{\overset{1}{n}}`$, $`\underset{1}{\overset{2}{n}}`$ and $`\underset{1}{\overset{3}{n}}`$. Using the formulas (79) we can see, that the functions $`\underset{0}{\overset{1}{n}}`$ $`\underset{0}{\overset{2}{n}}`$ are analytic of variable $`\overline{k}`$. The obvious form of this dependence is defined by two conditions. First one is the continuity of the asymptotics and second one is matching condition for the external and internal asymptotic expansions.

One can see, that the sufficient condition for the continuity of the coefficient of the asymptotics (76) as $`t^1`$ with respect to $`k`$ is an addition into $`\underset{0}{\overset{2}{n}}`$ of the term:

$$\underset{0}{\overset{2}{n^{}}}=\frac{1}{\pi i}\frac{1}{12i\overline{(kk_0)}^2}_i\mathrm{}^i\mathrm{}\frac{d\lambda }{\overline{k\lambda }}f(\overline{l}).$$
(83)
It follows from the formulas (81), that the coefficient of the asymptotics (75) as $`t^{5/3}`$ is defined after evaluating of $`\underset{0}{\overset{1}{n}}`$, i.e. after matching of the coefficients as $`t^{2/3}`$ of the external and internal expansions.

Consider the problem for defining of $`\stackrel{1}{m}`$ in detail. This function satisfies the boundary condition:

$$\underset{0}{\overset{1}{m}}|_{|k|\mathrm{}}=0.$$

The external solution is valid outside of small neighborhoods of the points $`k_j`$, $`j=0,1,2.`$ Therefore solutions of the equation (78), which are smooth, bounded and decreased as $`|k|\mathrm{}`$:

$`\underset{0}{\overset{s}{\stackrel{1}{m}}}(k,\theta )={\displaystyle \frac{1}{2i\pi }}{\displaystyle _{\text{}}\frac{drd\overline{r}}{kr}\frac{f(\overline{r})f(r)f(\overline{k}_0)f(k_0)}{12i(rk_0)^2}}+`$
$`+{\displaystyle \frac{\overline{(kk_0)}f(k_0)f(\overline{k}_0)}{12i(kk_0)^2}}`$ (84)
are defined to within an analytic function with respect to variable $`k`$, which has poles when $`t^{1/3}|kk_0|1`$. The full definition of $`\underset{0}{\overset{1}{m}}`$ will be done by matching of the external and internal expansions.

For matching process we need asymptotic behaviour of the partial solutions of equations (67)-(70) as $`|p|\mathrm{}`$. Constructing of these asymptotics is reduced to evaluating of an integrals with weak singularity in the integrand. Evaluating of these integrals is done in the Appendix. Here we write the asymptotic behaviour of partial solutions of equations (67)-(70) as $`|p|\mathrm{}`$.

$`\stackrel{s}{\stackrel{1}{𝒩}}(p,\xi ,\eta )={\displaystyle \frac{1}{\overline{p}}}\mathrm{exp}(i\omega (p))f(\overline{k}_0)\left[\varphi _{00}^+(v^2)\varphi _{00}^{}(v^2)\right]+`$
$`+{\displaystyle \frac{1}{\overline{p}^2}}\mathrm{exp}(i\omega (p))f(\overline{k})\left[\varphi _{01}^+(v^2)\varphi _{01}^{}(v^2)\right]+`$
$`+{\displaystyle \frac{1}{12i}}f(\overline{k}_0)\text{sgn}\left[\text{Re}(\overline{p})\right]\left({\displaystyle \frac{1}{p^2}}{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{\overline{p}^2}}\right)+O(|v^2||p|^3).`$
$`\stackrel{s}{\stackrel{2}{𝒩}}(p,\xi ,\eta )=[{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{\overline{p}}}(\varphi _{10}^+(v^2)\varphi _{10}^{}(v^2))+`$
$`({\displaystyle \frac{1}{p}}{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{\overline{p}}}){\displaystyle \frac{\text{sgn}[\text{Re}(\overline{p})]}{12i}}]f_{10}^{(1)}(\overline{k}_0)+`$
$`\left[{\displaystyle \frac{1}{\overline{p}}}\mathrm{exp}(i\omega (p))\varphi _{10}(v^2)+{\displaystyle \frac{1}{12ip}}\right]f_{10}^{(2)}(\overline{k}_0)+`$
$`[{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{\overline{p}}}(\varphi _{01}^+(v^2)\varphi _{01}^{}(v^2))+({\displaystyle \frac{\overline{p}}{p^2}}{\displaystyle \frac{\mathrm{exp}(\omega (p))}{\overline{p}}})\times `$
$`{\displaystyle \frac{\text{sgn}\left[\text{Re}(\overline{p})\right]}{12i}}]f_{01}^{(1)}(\overline{k}_0)+[{\displaystyle \frac{1}{\overline{p}}}\mathrm{exp}(i\omega (p))\varphi _{01}(v^2)+`$
$`{\displaystyle \frac{\overline{p}}{12ip^2}}]f^{(2)}_{01}(\overline{k}_0)+O(|p|^2)+O(|v|^2|p|^2),`$
$`\stackrel{s}{\stackrel{1}{}}(p,\xi ,\eta )=f(k_0)f(\overline{k}_0)[{\displaystyle \frac{\overline{p}}{12ip^2}}+`$
$`{\displaystyle \frac{1}{p}}\left(12i\varphi _{01}(v^2)\overline{\varphi _{10}(v^2)}+iv^2\varphi _{00}(v^2)\overline{\varphi _{00}(v^2)}\right)`$
$`{\displaystyle \frac{2}{p}}[iv^2(\varphi _{00}^+(v^2)\psi _{00}^{}(v^2)+\varphi _{00}^{}(v^2)\psi _{00}^+(v^2))+`$
$`+12i(\varphi _{01}^+(v^2)\psi _{01}^{}(v^2)+\varphi _{01}^{}(v^2)\psi _{01}^+(v^2)){\displaystyle \frac{1}{2}}\psi _{01}]]+`$
$`O(|v|^2|p|^2);`$
$`\stackrel{s}{\stackrel{2}{}}(p,\xi ,\eta )=f(k_0)f_{10}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{p}}{12ip}}+`$
$`f(k_0)f_{10}^{(2)}(\overline{k}_0)\left({\displaystyle \frac{\overline{p}}{12ip}}+{\displaystyle \frac{1}{12i}}\right)\text{sgn}[\text{Re}(p)]+`$
$`f(k_0)f_{01}^{(2)}(\overline{k}_0)\left({\displaystyle \frac{\overline{p}^2}{24ip^2}}{\displaystyle \frac{1}{24i}}\right)\text{sgn}[\text{Re}(p)]+`$
$`f(k_0)f_{01}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{p}^2}{24ip^2}}+\mathrm{\Phi }(v^2)+O(|v|^2|p|^1),`$
$`\stackrel{s}{\stackrel{3}{𝒩}}(p,\xi ,\eta )={\displaystyle \frac{1}{24i}}\left[\text{sgn}[\text{Re}(\overline{p})]\left(1\mathrm{exp}(i\omega (p))\right)f_{20}^{(1)}(\overline{k}_0)\right]+`$
$`\left[\text{sgn}[\text{Re}(\overline{p})]\left({\displaystyle \frac{\overline{p}^2}{12ip^2}}{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{24i}}\right)f_{02}^{(1)}(\overline{k}_0)+{\displaystyle \frac{\overline{p}^2f_{02}^{(2)}(\overline{k}_0)}{24ip^2}}\right]+`$
$`\left[\text{sgn}[\text{Re}(\overline{p})]\left({\displaystyle \frac{\overline{p}}{12ip}}+{\displaystyle \frac{\mathrm{exp}(i\omega (p))}{12i}}\right)f_{11}^{(1)}(\overline{k}_0)+{\displaystyle \frac{\overline{p}}{12ip}}f_{11}^{(2)}(\overline{k}_0)\right]+`$
$`f_0^2(\overline{k}_0)f(k_0)\mathrm{exp}(i\omega (p))\mathrm{\Psi }(v^2)+O(|v||p|^1).`$
The functions $`\mathrm{\Phi }(v^2)`$ and $`\mathrm{\Psi }(v^2)`$ are smooth and uniformly bounded when $`v^2\text{}`$. Here we use notations:

$$\varphi _{mn}^\pm =_{\text{Re}(r)>0}𝑑rd\overline{r}r^m\overline{r}^n\mathrm{exp}(i\omega (r)),$$

$$\psi _{mn}^\pm =_{\text{Re}(r)>0}𝑑rd\overline{r}r^m\overline{r}^n\mathrm{exp}(i\omega (r)),$$

$$\varphi _{mn}=_{\text{}}𝑑rd\overline{r}r^m\overline{r}^n\mathrm{exp}(i\omega (r)),$$

$$\psi _{mn}=_{\text{}}𝑑rd\overline{r}r^m\overline{r}^n\mathrm{exp}(i\omega (r)).$$
(85)
Evaluate an asymptotic behaviours of the coefficients of the external expansion as $`kk_0`$ and $`_k^2S=o(1)`$.

$`\underset{1}{\overset{1}{n}}|_{kk_0}={\displaystyle \frac{f(\overline{k}_0)}{12i(kk_0)^2}}\text{sgn}[\text{Re}(\overline{k})]+`$
$`{\displaystyle \frac{1}{12i(kk_0)}}[f_{10}^{(1)}(\overline{k})+f_{01}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}}{kk_0}}+`$
$`(f_{10}^{(2)}(\overline{k}_0)+f_{01}^{(2)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}}{kk_0}})\text{sgn}[\text{Re}(\overline{k})]]+`$
$`+{\displaystyle \frac{1}{2}}f_{20}^{(1)}(\overline{k})+f_{11}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}}{kk_0}}+{\displaystyle \frac{1}{2}}f_{02}^{(1)}{\displaystyle \frac{\overline{(kk_0)^2}}{(kk_0)^2}}+`$
$`({\displaystyle \frac{1}{2}}f_{20}^{(2)}(\overline{k}_0)+f_{11}^{(2)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}}{kk_0}}+`$
$`{\displaystyle \frac{1}{2}}f_{02}^{(2)}(\overline{k}_0){\displaystyle \frac{\overline{(kk_0)^2}}{(kk_0)^2}})\text{sgn}[\text{Re}(\overline{k})]+o(1).`$
$`\underset{0}{\overset{}{\stackrel{2}{n}}}={\displaystyle \frac{1}{24\pi \overline{(kk_0)^2}}}\times `$
$`[\pi if(\overline{k}_0)\text{sgn}[\text{Re}(\overline{p})]+V.P.{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{d\lambda f(i\lambda )}{\lambda ik_0}}]+`$
$`{\displaystyle \frac{1}{12\pi \overline{(kk_0)}}}[\pi if_{01}^{(1)}(\overline{k}_0)+\pi if_{10}^{(1)}(\overline{k}_0)+(\pi if_{10}(2)(\overline{k}_0)+`$
$`+\pi if_{01}^{(2)}(\overline{k}_0))\text{sgn}[\text{Re}(\overline{k}_0)]+`$
$`+V.P.{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{d\lambda f(i\lambda )}{(\lambda ik_0)^2}}]+{\displaystyle \frac{1}{12i}}[\pi if_{20}^{(1)}(\overline{k}_0)+`$
$`+2\pi if_{11}^{(1)}(\overline{k}_0)+\pi if_{02}^{(1)}(\overline{k}_0)+(\pi if_{20}^{(2)}(\overline{k}_0)+2\pi if_{11}^{(2)}(\overline{k}_0)+`$
$`+\pi if_{02}^{(2)}(\overline{k}_0))\text{sgn}[\text{Re}(\overline{k})]+V.P.{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{d\lambda f(i\lambda )}{(\lambda ik_0)^3}}]+o(1).`$
$`\stackrel{s}{\stackrel{1}{m}}={\displaystyle \frac{\overline{(kk_0)}}{12i(kk_0)^2}}f(\overline{k}_0)f(k_0)`$
$`{\displaystyle \frac{\overline{kk_0}}{12i(kk_0)}}[(f_{10}^{(1)}(\overline{k}_0)\text{sgn}[\text{Re}(\overline{k})]+`$
$`+f_{10}^{(2)}(\overline{k}_0))f(k_0)+f(\overline{k}_0)(f_{10}^{(1)}(k_0)\text{sgn}[\text{Re}(k)]+f_{10}^{(2)}(k_0))]`$
$`{\displaystyle \frac{\overline{(kk_0)^2}}{12i(kk_0)^2}}[f(k_0)(f_{01}^{(1)}(\overline{k}_0)\text{sgn}[\text{Re}(\overline{k})]+f_{01}^{(2)}(\overline{k}_0))+`$
$`+f(\overline{k}_0)(f_{01}^{(1)}(k_0)\text{sgn}[\text{Re}(k)]+f_{01}^{(2)}(k_0))]+o(1).`$
Let us do matching of the external and internal asymptotic expansions of the function $`\stackrel{~}{\nu }`$. The matching condition for $`\stackrel{~}{\nu }`$ in domain $`t^{1/3}|kk_0|1`$ when $`|\theta |1`$ has the form:

$`\left(t^{2/3}\underset{0}{\overset{1}{n}}+t^1\underset{0}{\overset{2}{n}}+t^1\stackrel{1}{n_1}\mathrm{exp}(itS)\right)`$
$`(t^{1/3}\stackrel{1}{𝒩}+t^{2/3}\stackrel{2}{𝒩}+t^1\stackrel{3}{𝒩})=o(t^1).`$
Let us equate the coefficients with equal powers of the large parameter $`t`$. As a result we obtain:

$$\underset{0}{\overset{1}{n}}=\frac{1}{\overline{kk_0}}f(\overline{k}_0)[\varphi _{00}^+(v^2)\varphi _{00}^{}(v^2)],$$
(86)
$`\underset{0}{\overset{2}{n}}=`$ $`{\displaystyle \frac{1}{\overline{kk_0}^2}}f(\overline{k})[\varphi _{01}^+(v^2)\varphi _{01}^{}(v^2)]+`$ (87)
$`{\displaystyle \frac{1}{\overline{kk_0}}}f_{10}^{(1)}(\overline{k}_0)[\varphi _{10}^+(v^2)\varphi _{10}^{}(v^2)]+{\displaystyle \frac{1}{\overline{kk_0}}}f_{10}^{(2)}(\overline{k}_0)\varphi _{10}(v^2)+`$
$`{\displaystyle \frac{1}{\overline{kk_0}}}f_{01}^{(1)}(\overline{k}_0)[\varphi _{01}^+(v^2)\varphi _{01}^{}(v^2)]+`$
$`+{\displaystyle \frac{1}{\overline{kk_0}}}f_{01}^{(2)}(\overline{k}_0)\varphi _{01}(v^2),`$
where $`\varphi `$, $`\psi `$, $`\varphi ^\pm `$ $`\psi ^\pm `$ are defined by (85);

$`\stackrel{1}{𝒩}=𝒫[g_1],\stackrel{2}{𝒩}=𝒫[g_2],`$
$`\stackrel{3}{𝒩}=𝒫[g_3]f_0^2(\overline{k}_0)f(k_0)\mathrm{exp}(i\omega (p))\mathrm{\Psi }(v^2),`$ (88)
where $`g_j`$ are the right-hand sides of the equations for $`\stackrel{j}{𝒩}`$, the operator $`𝒫[g]`$ is defined by (72), the function $`\mathrm{\Psi }(v^2)`$ is smooth and uniformly bounded with respect to $`v^2`$ when $`v^2\text{}`$.

Matching condition for the function $`\stackrel{~}{\mu }`$ in the domain $`t^{1/3}|kk_0|1`$ when $`|\theta |1`$ has the form:

$$(1+t^1\stackrel{1}{m})(1+t^{2/3}\stackrel{1}{}+t^1\stackrel{2}{})=o(t^1).$$

As a result we obtain:

$$\stackrel{1}{}=\stackrel{s}{\stackrel{1}{}};\stackrel{2}{}=\stackrel{s}{\stackrel{2}{}}\mathrm{\Phi }(v^2),$$
(89)
where $`\stackrel{s}{\stackrel{1}{}}`$ is defined by (73), the function $`\stackrel{s}{\stackrel{2}{}}`$ is defined by (74) and the function $`\mathrm{\Phi }(v^2)`$ is smooth and uniformly bounded when $`v^2\text{}`$.

$`\stackrel{1}{m}=\stackrel{s}{\stackrel{1}{m}}{\displaystyle \frac{1}{kk_0}}(12i\varphi _{01}(v^2)\overline{\varphi _{10}(v^2)}+iv^2\varphi _{00}(v^2)\overline{\varphi _{00}(v^2)})`$
$`{\displaystyle \frac{2}{kk_0}}[iv^2(\varphi _{00}^+(v^2)\psi _{00}^{}(v^2)+\varphi _{00}^{}(v^2)\psi _{00}^+(v^2))+`$
$`+12i(\varphi _{01}^+(v^2)\psi _{01}^{}(v^2)+\varphi _{01}^{}(v^2)\psi _{01}^+(v^2)){\displaystyle \frac{1}{2}}\psi _{01}].`$ (90)
Here $`\stackrel{s}{\stackrel{1}{m}}`$ is defined by (84) and the functions $`\psi ^\pm `$, $`\varphi ^\pm `$, $`\varphi `$ and $`\psi `$ are defined by formulas (85).

Thus we have matched internal and external expansions of $`\stackrel{~}{\mu }`$ and $`\stackrel{~}{\nu }`$. The lemma is proved.

These expansions are ununiform with respect to $`k`$. Now we can construct uniform asymptotic expansion with respect to $`k\text{}`$. Following the matching method the uniform expansion is:

$$\left(\begin{array}{c}\widehat{\mu }\\ \widehat{\nu }\end{array}\right)=\left(\begin{array}{c}\stackrel{~}{m}\\ \stackrel{~}{n}\end{array}\right)+\left(\begin{array}{c}\stackrel{~}{}_1\\ \stackrel{~}{𝒩}_1\end{array}\right)A_{1,k}\left(\begin{array}{c}\stackrel{~}{}_1\\ \stackrel{~}{𝒩}_1\end{array}\right).$$
(91)
Here the operator $`A_{n,k}`$ processes on the function $`\stackrel{~}{}_1`$ in the formula for $`_1`$ by followed manner. One must change the variable $`p`$ into the variable $`k`$ using the formula (64) and write all terms of the asymptotic expansion with respect to $`t`$ with the powers are equal to $`m`$, where $`0mn`$. For example, one can obtain for the function $`\stackrel{~}{}(p,\xi ,\eta ,t)`$:

$`A_{1,k}[\stackrel{~}{}(p,\xi ,\eta ,t)]=1+t^1(f(k_0)f(\overline{k}_0)[{\displaystyle \frac{\overline{kk_0}}{12i(kk_0)^2}}+`$
$`+{\displaystyle \frac{1}{kk_0}}\left(12i\varphi _{01}(v^2)\overline{\varphi _{10}(v^2)}+iv^2\varphi _{00}(v^2)\overline{\varphi _{00}(v^2)}\right)`$
$`{\displaystyle \frac{2}{kk_0}}[iv^2(\varphi _{00}^+(v^2)\psi _{00}^{}(v^2)+\varphi _{00}^{}(v^2)\psi _{00}^+(v^2))+`$
$`+12i(\varphi _{01}^+(v^2)\psi _{01}^{}(v^2)+\varphi _{01}^{}(v^2)\psi _{01}^+(v^2)){\displaystyle \frac{1}{2}}\psi _{01}]]`$
$`f(k_0)f_{10}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}}{12i(kk_0)}}+`$
$`f(k_0)f_{10}^{(2)}(\overline{k}_0)\left({\displaystyle \frac{\overline{kk_0}}{12i(kk_0)}}+{\displaystyle \frac{1}{12i}}\right)\text{sgn}[\text{Re}(k)]+`$
$`f(k_0)f_{01}^{(1)}(\overline{k}_0){\displaystyle \frac{\overline{kk_0}^2}{24i(kk_0)^2}}+`$
$`+f(k_0)f_{01}^{(2)}(\overline{k}_0)({\displaystyle \frac{\overline{kk_0}^2}{24i(kk_0)^2}}{\displaystyle \frac{1}{24i}})\text{sgn}[\text{Re}(k)]),`$
The obviously formula for $`A_{1,k}[\stackrel{~}{𝒩}]`$ is more large and doesn’t shown here.

###### Theorem 5

The formula (91) gives the asymptotics solution of the problem (9) with respect to $`mod(O(t^1))`$ as $`t\mathrm{}`$. This asymptotic solution is uniform with respect to $`k\text{}`$ and $`|\theta |1`$.

## 5 Justification of the asymptotics of the solution of $`\overline{D}`$-problem

In this section we prove that the remainder of the asymptotics has order by $`t^{4/3}`$ uniformly with respect to $`k\text{}`$ and this remainder has to be differentiable with respect to $`x`$. We call by the remainder of the asymptotics the difference between solution of the problem (9) and constructed asymptotic solutions (47) when $`\theta ^2t^{2/3}1`$, (63) when $`\theta ^2t^{2/3}1`$ and (91) when $`|\theta |1`$. The differentiability of the remainder will be important when we will construct an asymptotic behaviour of solution of the equation KP-2.

###### Theorem 6

Let $`^\alpha f(k,\overline{k})L_1C^1`$ when $`|\alpha |2`$ when $`k`$ is out of the imaginary axis and

$$\underset{z\text{}}{sup}\left|_{\text{}}\frac{dkd\overline{k}}{|kz|}|F(k)|\right|<\mathrm{\hspace{0.17em}2}\pi ,$$

then the solution of the problem (9) is:

$$\left(\begin{array}{c}\mu \\ \nu \end{array}\right)=\left(\begin{array}{c}\stackrel{~}{\mu }\\ \stackrel{~}{\nu }\end{array}\right)+O(t^{{}_{4}{}^{}/3}),$$
(92)
when $`k\text{},\xi ,\eta \text{}`$. The remainder of the asymptotics has to be differentiable with respect to $`x`$.

The proof. Let us write the system of differential equation for the remainder. Denote the remainder in (97) by $`V`$. Substitute (97) into (9). As a result we obtain:

$$\left(\begin{array}{cc}_{\overline{k}}& 0\\ 0& _k\end{array}\right)V=\left(\begin{array}{cc}0& F(\overline{k})\mathrm{exp}(itS)\\ F(k)\mathrm{exp}(itS)& 0\end{array}\right)V+f,$$
(93)
$$V|_{|k|\mathrm{}}=0.$$
(94)
We denote by the vector $`f`$ the residual which originates in (9) when we substitute the column into this equation $`(\stackrel{~}{\mu },\stackrel{~}{\nu })^T`$:

$$f_1=_{\overline{k}}\stackrel{~}{\mu }+F(\overline{k})\mathrm{exp}(itS)\stackrel{~}{\nu },$$

$$f_2=_k\stackrel{~}{\nu }+F(k)\mathrm{exp}(itS)\stackrel{~}{\mu }.$$

Let us denote by $`X`$ the space of bounded and continuous with respect to $`k`$ vector-functions with the norm:

$$W=\underset{k\text{},(\xi ,\eta )\text{}^2}{sup}|W_1|+\underset{k\text{},(\xi ,\eta )\text{}^2}{sup}|W_2|.$$

Consider a system of integral equations instead of the problem (93), (94):

$$V=G[F]V+H,$$
(95)
where $`G[F]`$ is the integral operator:

$`G[F]V={\displaystyle }{\displaystyle _m\text{}}dmd\overline{m}\times `$
$`\left(\begin{array}{cc}0& \frac{F(\overline{m})}{km}\mathrm{exp}(itS)\\ \frac{F(m)}{\overline{km}}\mathrm{exp}(itS)& 0\end{array}\right)V(m,\xi ,\eta ,t);`$
$$H=_m\text{}𝑑md\overline{m}\left(\begin{array}{c}\frac{f_1(m,\xi ,\eta ,t)}{km}\\ \frac{f_2(m,\xi ,\eta ,t)}{\overline{km}}\end{array}\right).$$

Using obvious forms of the functions $`f_{1,2}`$ one can proof, that

$$H=O(t^{4/3})$$

uniformly with respect to $`\xi ,\eta \text{}`$ when $`k\text{}`$.

The operator $`G[F]`$ is contracting in the space $`X`$, therefore the solution of the integral equation (95) exists in $`X`$ and is evaluated by $`O(t^{4/3})`$ uniformly with respect to $`\xi ,\eta \text{}`$.

Show that the remainder of the asymptotics is differentiable on $`x=t\xi `$. Differentiate with respect to $`x`$ the system of the equation for the remainder. Denote the derivative of the vector $`V`$ by $`\chi `$. Then we obtain:

$$\chi =G[F]\chi +_xG[F]V+_xH.$$

The terms $`_xG[F]V+_xH`$ may be evaluated by order $`O(t^{4/3})`$. The operator $`G[F]`$ is contracting, therefore one can obtain: $`\chi =O(t^{4/3})`$. The theorem is proved.

## 6 Solution of the equation KP-2

An asymptotics of the solution of the problem (5) as $`t\mathrm{}`$ may be written as:

$$\left(\begin{array}{c}\varphi \\ \psi \end{array}\right)=\left(\begin{array}{c}\stackrel{~}{\mu }\\ \stackrel{~}{\nu }\end{array}\right)(k,\overline{k},\xi ,\eta ,t)+\left(\begin{array}{c}\stackrel{~}{\nu }\\ \stackrel{~}{\mu }\end{array}\right)(\overline{k},k,\xi ,\eta ,t)+O(t^{4/3}).$$
(97)
The second term in this formula is the solution of the problem (5) with the boundary condition:

$$\left(\begin{array}{c}\stackrel{~}{\mu }\\ \stackrel{~}{\nu }\end{array}\right)|_{|k|\mathrm{}}=\left(\begin{array}{c}0\\ 1\end{array}\right).$$

The proof of the theorem 1. Let us substitute the function $`\psi (k,\overline{k},\xi ,\eta ,t)`$ into the formula for the solution of the equation KP-2(6). Differentiate the integrand with respect to $`x`$. The main terms of the integrand are the terms which appear after differentiating of the exponent. The derivatives with respect to $`x`$ of others factors of integrand are small because they depend on $`x`$ slowly. The integrals of such terms over the plane are evaluated by the order $`O(t^{4/3})`$ uniformly with respect to $`(x,y)\text{}`$. Rewrite the main term of the integral as the integral over real plane: $`(\kappa ,\lambda )\text{}^2,`$ where $`\kappa =\text{Re}(k),\lambda =\text{Im}(k)`$. As a result we obtain:

$`u(x,y,t)=4{\displaystyle }{\displaystyle _\text{}^2}d\kappa d\lambda |\kappa |f(\kappa +i\lambda )\times `$
$`\times \mathrm{exp}(it(8\kappa ^324\kappa \lambda ^2+2\kappa \xi +4\kappa \lambda \eta ))+O(t^{4/3}).`$ (98)
Thus the main term of the asymptotics of the solution of the Cauchy problem for the equation KP-2 is given by integral with fast oscillating exponent. Let us evaluate the asymptotic behaviour of this integral.

Let $`\eta ^2+12\xi >0`$ and $`t^{1/3}|\eta ^2+12\xi |1`$, then the stationary points of the exponent are: $`\lambda _{1,2}=\frac{\eta }{12}\pm \sqrt{\eta ^2+12\xi },\kappa _{1,2}=0`$. The stationary phase method (see for example ) gives:

$$u(x,y,t)=o(t^1).$$

Let $`\eta ^2+12\xi <0`$ and $`t^{1/3}|\eta ^2+12\xi |1`$, then $`\lambda _{1,2}=\frac{\eta }{12},\kappa _{1,2}=\pm \frac{1}{2}\sqrt{\eta ^212\xi }.`$ The asymptotic behaviour of the integral is:

$`u(x,y,t)=4t^1{\displaystyle \frac{\pi }{12i\sqrt{\eta ^212\xi }}}f_0({\displaystyle \frac{1}{2}}\sqrt{\eta ^212\xi }+{\displaystyle \frac{i\eta }{12}})\times `$
$`\times \mathrm{exp}(11it\sqrt{{\displaystyle \frac{y^2}{t^2}}12{\displaystyle \frac{x}{t}}})+c.c.+o(t^1).`$
To evaluate the main term of the asymptotics of the integral when $`|12\xi +\eta ^2|=o(1)`$, we substitute the scaled variables $`p_1=t^{1/3}\text{Re}(kk_0)`$, $`p_2=t^{1/3}\text{Im}(kk_0)`$ and the parameter $`v^2=t^{2/3}(\eta ^2+12\xi )/\sqrt{12}`$ into integral (98). As a result we obtain:

$`u(\xi ,\eta ,t)=4it^1f(k_0){\displaystyle }{\displaystyle _\text{}^2}dp_1dp_2\times `$
$`p_1\mathrm{exp}(i(8p_1^32v^2p_124p_1p_2^2))+o(t^1).`$
Let us integrate the internal integral with respect to the parameter $`p_2`$, use the even property of the integrand with respect to $`p_1`$. As a result we obtain:

$`u(x,y,t)=8it^1\sqrt{\pi }f(i\eta /12)({\displaystyle _0^{\mathrm{}}}dp_1\sqrt{p_1}\mathrm{cos}\left(8p_1(p_1^28v^2)\right)+`$
$`+{\displaystyle _0^{\mathrm{}}}dp_1\sqrt{p_1}\mathrm{sin}\left(8p_1(p_1^28v^2)\right))+o(t^1).`$
The theorem 1 is proved.

## Appendix A Asymptotic behaviour of double integral with weak singularity of the integrand

Here we obtain the asymptotic behaviour of integrals which are appeared when the asymptotic solution of (9) was studied.

Evaluating of the asymptotic behaviour of one-dimensional integrals with weak singular integrand and fast oscillated exponent was done in  p.26 and  p.332. The asymptotic behaviour of many-dimensional integrals with fast oscillating exponent was studied in , . The asymptotic behaviour of the Cauchy integrals with fast oscillating exponent in one-dimensional case was studied in  and in many-dimensional case was studied in . The asymptotic behaviour of the two-dimensional for some integrals over all complex plane with weak singularity was studied in .

### A.1 Integrals over half-plane with general stationary point of fast oscillated exponent

Here we study an asymptotic behaviour of an integral:

$$I=_{\mathrm{\Omega }^+}\frac{dnd\overline{n}}{ln}\mathrm{exp}(i(n^2+\overline{n}^2)),$$

where $`|l|\mathrm{}`$ and the domain $`\mathrm{\Omega }^+=\{\text{Re}(l+\overline{l}i(l\overline{l}))>0\}`$.

###### Theorem 7

The asymptotic behaviour of the integral $`I`$ as $`|p|\mathrm{}`$ has the form:

$$I=2i\pi \frac{\mathrm{exp}(i(l^2+\overline{l}^2))}{2il}\frac{3i\pi }{2\overline{l}}+O(|l|^2)\text{ѳ }l\mathrm{\Omega }^+;$$

$$I=\frac{i\pi }{2\overline{l}}+O(|l|^2)\text{where}l\mathrm{\Omega }^+.$$

The proof. Let us suppose that $`l\mathrm{\Omega }^+`$. Divide the domain $`\mathrm{\Omega }^+`$ into three domains: first one is $`\mathrm{\Omega }_1^+=\{\mathrm{\Omega }^+\backslash \{[|n||l|/2][|nl|<\epsilon ]\}\}`$. This domain hasn’t stationary point of phase function of the exponent and the singularity of the integrand. Second one is $`\mathrm{\Omega }_2^+=\{\mathrm{\Omega }^+\backslash [|n||l|/2]\}`$. This domain contains the stationary point of the phase function. At last, third domain is $`\mathrm{\Omega }_3^+=\{|nl|<\epsilon \}`$. This domain contains the singularity of the integrand.

Let us integrate by parts over $`\mathrm{\Omega }_1^+`$. As a result we obtain:

$$_{\mathrm{\Omega }_1^+}\frac{dnd\overline{n}}{\overline{ln}}\mathrm{exp}(i(n^2+\overline{n}^2))=_{\mathrm{\Omega }_1^+}\frac{d\overline{n}}{\overline{ln}}\frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2in}$$

$$_{\mathrm{\Omega }_1^+}\frac{dnd\overline{n}}{\overline{ln}}\frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2in^2}.$$
(99)
The boundary of the domain $`\mathrm{\Omega }_1^+`$ includes a large half-circle of radius $`R`$ as $`R\mathrm{}`$, a circle of radius $`\epsilon `$ over the point $`l=n`$, a half-circle of radius $`|l|/2`$ and two segments: $`[R\mathrm{exp}(3i\pi /4),|l|\mathrm{exp}(3i\pi )/2]`$ and $`[|l|\mathrm{exp}(i\pi /4)/2,R\mathrm{exp}(3i\pi )]`$.

Consider the integrals over boundary of the domain $`\mathrm{\Omega }_1^+`$. The integral over large half-circle is equals to zero as $`R\mathrm{}`$. The integral over half-circle of radius $`|l|/2`$ has order $`|l|^3`$ (because of oscillations and small value of the integrand). The integral over the circle at $`n=l`$ equals to a residue of the integrand multiplied by $`2i\pi `$ as $`\epsilon 0`$. The integral over the segments $`[R\mathrm{exp}(3i\pi /4),|l|\mathrm{exp}(3i\pi )/2]`$ and $`[|l|\mathrm{exp}(i\pi /4)/2,R\mathrm{exp}(3i\pi )]`$ (let us denote their union by $`L`$) is reduced to the form:

$$_{nL}\frac{d\overline{n}}{2in\overline{(ln)}}=\frac{\mathrm{exp}(3i\pi /4)}{2i\overline{l}}_{|\lambda ||l|/2}\frac{d\lambda }{\overline{l}\lambda \mathrm{exp}(i\pi /4)}.$$

The second term of (99) has to be evaluated as by $`O(|l|^3)`$.

Let us consider the integral over $`\mathrm{\Omega }_2^+`$. Represent:

$$\frac{1}{\overline{ln}}=\frac{1}{\overline{l}}\left(1+\frac{\overline{n}}{\overline{ln}}\right).$$

Then the integral over $`\mathrm{\Omega }_2^+`$ has to be written as:

$`I_2={\displaystyle _{\mathrm{\Omega }_2^+}\frac{dnd\overline{n}}{\overline{ln}}\mathrm{exp}(i(n^2+\overline{n}^2))}=`$
$`{\displaystyle \frac{1}{\overline{l}}}{\displaystyle _{\mathrm{\Omega }_2^+}\frac{dnd\overline{n}}{\overline{ln}}\mathrm{exp}(i(n^2+\overline{n}^2))}+`$
$`{\displaystyle \frac{1}{\overline{l}}}{\displaystyle _{\mathrm{\Omega }_2^+}\frac{dnd\overline{n}}{\overline{ln}}\overline{n}\mathrm{exp}(i(n^2+\overline{n}^2))}.`$
Here we integrate by parts the second term. After evaluations we obtain:

$`I_2={\displaystyle \frac{1}{\overline{l}}}{\displaystyle _{\mathrm{\Omega }^+}𝑑n}d\overline{n}\mathrm{exp}(i(n^2+\overline{n}^2))+`$
$`{\displaystyle \frac{1}{2i\overline{l}}}{\displaystyle _{|l|\mathrm{exp}(3i\pi /4)/2}^{|l|\mathrm{exp}(i\pi /4)/2}}{\displaystyle \frac{dn}{\overline{ln}}}+O(|l|^2).`$
The integral over $`\mathrm{\Omega }_3^+`$ as $`\epsilon 0`$ equals to zero.

Let us sum the obtained asymptotics:

$$I=2i\pi \frac{\mathrm{exp}(i(l^2+\overline{l}^2))}{2il}+\frac{\mathrm{exp}(i\pi /2)}{2i\overline{l}}_{\mathrm{}}^{\mathrm{}}\frac{d\lambda }{\overline{l}\mathrm{exp}(i\pi /4)\lambda }+$$

$$\frac{1}{2i\overline{l}}_{\mathrm{\Omega }^+}𝑑nd\overline{n}\mathrm{exp}(i(n^2+\overline{n}^2))+O(|l|^2).$$

Thus the first statement of the theorem is proved. The second statement has to be proved by the same way.

The theorem is proved.

### A.2 Asymptotic behaviour of the integral <br>with confluent phase function

In this section we obtain an asymptotic behaviour of an integral as $`|p|\mathrm{}`$

$$W^+=_{\mathrm{\Omega }^+}\frac{drd\overline{r}}{\overline{pr}}\mathrm{exp}(i\omega (r)).$$
(100)
###### Theorem 8

The asymptotic behaviour of the integral (100) where $`\omega (p)`$ $`=4(p^3+\overline{p}^3)v^2(p+\overline{p})`$ as $`|p|\mathrm{}`$ and $`|p\pm \frac{v}{\sqrt{12}}|0`$ has the form:

$$W^+=[\begin{array}{c}\frac{1}{\overline{p}}\varphi _{00}^+(v^2)+\frac{1}{\overline{p}^2}\varphi _{01}^+(v^2)+\frac{\pi i}{12\overline{p}^2}+2\pi i\frac{\mathrm{exp}(i\omega (p))}{12p^2}+\\ +O(|p|^3+|v|^2|p|^2),\text{when}p\mathrm{\Omega }^+;\\ \\ \frac{1}{\overline{p}}\varphi _{00}^+(v^2)+\frac{1}{\overline{p}^2}\varphi _{01}^+(v^2)\frac{\pi i}{12\overline{p}^2}+\\ +O(|p|^3+|v|^2|p|^2),\text{when}p\mathrm{\Omega }^+.\end{array}$$
(101)
Let us prove the theorem. Represent the integral in the form:

$`W^+={\displaystyle \frac{1}{\overline{p}}}{\displaystyle _{\mathrm{\Omega }^+}𝑑n}d\overline{n}\mathrm{exp}(i\omega (n))+`$
$`+{\displaystyle \frac{1}{\overline{p}^2}}{\displaystyle _{\mathrm{\Omega }^+}𝑑n}d\overline{n}\mathrm{exp}(i\omega (n))\overline{n}+`$
$`+{\displaystyle \frac{v^2}{12\overline{p}^2}}{\displaystyle _{\mathrm{\Omega }^+}\frac{dnd\overline{n}}{\overline{pn}}\overline{n}^2\mathrm{exp}(i\omega (n))}+`$
$`+{\displaystyle \frac{1}{12\overline{p}^2}}{\displaystyle _{\mathrm{\Omega }^+}\frac{dnd\overline{n}}{\overline{pn}}(12\overline{n}^2v^2)\mathrm{exp}(i\omega (n))}.`$
Integrate by parts the last term. As a result we obtain:

$`W^+={\displaystyle \frac{1}{\overline{p}}}\varphi _{00}^+(v^2)+{\displaystyle \frac{1}{\overline{p}^2}}\varphi _{01}^+(v^2)+`$
$`+{\displaystyle \frac{1}{12\overline{p}^2}}{\displaystyle _{\mathrm{\Omega }^+}}{\displaystyle \frac{d\overline{n}}{\overline{pn}}}{\displaystyle \frac{(12\overline{n}^2v^2)\mathrm{exp}(i\omega (n)))}{(i)(12n^2v^2)}}+`$
$`O(|p|^3+|v|^2|p|^2).`$
Consider the integral over the boundary of the domain $`\mathrm{\Omega }^+`$. This integral may be considered as a sum of integrals over imaginary axis $`\text{Re}(r)=0`$, over half-circle of radius $`R`$ as $`R\mathrm{},\text{Re}(r)>0`$ and over the circle $`|pr|=\epsilon `$ as $`\epsilon 0`$ (if $`p\mathrm{\Omega }^+`$). We can see that the integral over the large half-circle tends to zero. As a result we obtain the statement of the theorem.

One more double integral which we need is:

$$U^+=_{\mathrm{\Omega }^+}\frac{drd\overline{r}}{pr}\overline{r}\mathrm{exp}(i\omega (r)).$$

Its asymptotics can be evaluated by the same way as the asymptotics of the integral (100) as $`|p|\mathrm{}`$ $`|p\pm v|0`$:

$`U=[\begin{array}{c}\frac{1}{\overline{p}}\varphi _{10}^+(v^2)+2i\pi \frac{1}{12ip}\mathrm{exp}(i\omega (p))+\frac{\pi i}{12\overline{p}}+2\pi i\frac{\mathrm{exp}(i\omega (p))}{12p^2}+\\ +O(|p|^2+|v|^2|p|^2),\text{when}p\mathrm{\Omega }^+;\\ \\ \frac{1}{\overline{p}}\varphi _{10}^+(v^2)\frac{\pi i}{12\overline{p}}+O(|p|^2+|v|^2|p|^2),\text{when}p\mathrm{\Omega }^+.\end{array}`$
## Appendix B Reducing of four-multiply integral into double integral

This section is pure technical. Here we show as the four-multiply integrals have to be reduced into double integrals over half-plane and all complex plane.

### B.1 Four-multiply integral with nondegenerate phase of the exponent.

Let us show that the four-multiply integral may be written as a sum of double integrals. Change the variable: $`nm=r`$, then the integral $`J`$ has the form:

$`J={\displaystyle \frac{1}{2i\pi }}{\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dnd\overline{n}}{l_jn}}\times `$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle _{\text{}}\frac{drd\overline{r}}{\overline{r}}\mathrm{exp}(i(r^2+\overline{r}^2))\mathrm{exp}(2i(rn+\overline{rn}))}.`$
Integrate by part over $`\overline{n}`$. As a result we obtain:

$`J={\displaystyle \frac{1}{4\pi ^2}}\underset{R\mathrm{}}{lim}{\displaystyle _{|n|=R}}{\displaystyle \frac{dn\overline{n}}{l_jn}}\times `$
$`{\displaystyle _{\text{}}\frac{drd\overline{r}}{\overline{r}}\mathrm{exp}(i(r^2+\overline{r}^2))\mathrm{exp}(2i(rn+\overline{rn}))}+`$
$`{\displaystyle \frac{1}{2i\pi }}\underset{\epsilon 0}{lim}{\displaystyle _{|ln|=\epsilon }}{\displaystyle \frac{dn\overline{n}}{l_jn}}\times `$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle _{\text{}}\frac{drd\overline{r}}{\overline{r}}\mathrm{exp}(i(r^2+\overline{r}^2))\mathrm{exp}(2i(rn+\overline{rn}))}`$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle }{\displaystyle _{\text{}}}{\displaystyle \frac{dnd\overline{n}}{l_jn}}\overline{n}\times `$
$`{\displaystyle \frac{1}{2i\pi }}{\displaystyle _{\text{}}𝑑r}d\overline{r}(2i)\mathrm{exp}(i(r^2+\overline{r}^2))\mathrm{exp}(2i(rn+\overline{rn})).`$
First term is equal to zero because of fast oscillating of the integrand, second term is equal to residue with sign minus of the integrand at $`n=l`$, the internal integral in the third term has to be evaluated. Finally we obtain:

$$J=\overline{l}_j\frac{\mathrm{exp}(i(l_j^2+\overline{l}_j^2))}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{\overline{l_jm}}\mathrm{exp}(i(n^2+\overline{n}^2))\mathrm{exp}(i(l_j^2+\overline{l}_j^2)).$$

### B.2 Four-multiply integral over different half-planes with nondegenerate phase function

Let us consider an integral:

$`J_+={\displaystyle }{\displaystyle _\mathrm{\Omega }^{}}{\displaystyle \frac{dnd\overline{n}}{ln}}{\displaystyle \frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}.`$
A partially integrating over $`\overline{n}`$ gives:

$`J_+={\displaystyle _\mathrm{\Omega }_{}}{\displaystyle \frac{dn}{ln}}\overline{n}{\displaystyle \frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}`$
$`{\displaystyle }{\displaystyle _\mathrm{\Omega }^{}}{\displaystyle \frac{dnd\overline{n}}{ln}}\overline{n}{\displaystyle \frac{}{\overline{n}}}{\displaystyle \frac{\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}.`$
The derivative is:

$`J_+={\displaystyle _\mathrm{\Omega }^{}}{\displaystyle \frac{dn}{ln}}{\displaystyle \frac{\overline{n}\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}`$
$`{\displaystyle \frac{3}{4}}{\displaystyle _\mathrm{\Omega }^{}}{\displaystyle \frac{dn}{ln}}\mathrm{exp}(i(n^2+\overline{n}^2)).`$
Let us consider integral over part of the boundary in first term :

$`=\underset{R\mathrm{}}{lim}{\displaystyle _{R\mathrm{exp}(i\pi /4)}^{R\mathrm{exp}(3i\pi /4)}}{\displaystyle \frac{dn}{ln}}{\displaystyle \frac{\overline{n}\mathrm{exp}(i(n^2+\overline{n}^2))}{2i\pi }}\times `$
$`{\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i(m^2+\overline{m}^2))}.`$
Change the order of integration over $`m,\overline{m}`$ and $`n`$. The internal integral over $`n`$ has to be evaluate:

$$\underset{R\mathrm{}}{lim}_{R\mathrm{exp}(i\pi /4)}^{R\mathrm{exp}(3i\pi /4)}dn\frac{\overline{n}}{(ln)\overline{(nm)}}=[\begin{array}{c}\frac{1}{2}+\frac{\overline{m}}{il\overline{m}},\text{when}l\mathrm{\Omega }^+;\\ 1/2\text{when}l\mathrm{\Omega }^+.\end{array}$$

In last expression of $`J_+`$ the integrals over the large half of circle $`|n|=R`$ as $`R\mathrm{}`$ are equal to zero.

Final formulas for $`J_+`$ have to be written as:

$$J_+=[\begin{array}{c}\frac{5}{4}i\pi +il_{\mathrm{\Omega }^+}𝑑md\overline{m}\frac{\mathrm{exp}(i(m^2+\overline{m}^2))}{il\overline{m}},\text{when}l\mathrm{\Omega }^+;\\ \\ \frac{5}{4}i\pi +\frac{3}{2}i\pi \mathrm{exp}(i)l^2+\overline{l}^2))\overline{l}\mathrm{exp}(i(l^2+\overline{l}^2))\times \\ _{\mathrm{\Omega }^+}𝑑md\overline{m}\frac{\mathrm{exp}(i(m^2+\overline{m}^2))}{\overline{lm}},\text{when}l\mathrm{\Omega }^{}.\end{array}$$

### B.3 Four-multiply integral with confluent <br>phase function

Let us reduce an four-multiply integral

$$J_1=_{\text{}}\frac{dnd\overline{n}}{pn}\mathrm{exp}(i\omega (n))_{\text{}}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i\omega (m))$$

into sum of double integrals. Denote:

$$V(n,v^2)=\frac{\mathrm{exp}(i\omega (n))}{2i\pi }_{\text{}}\frac{drd\overline{r}}{\overline{nr}}\mathrm{exp}(i\omega (r)).$$

The function $`V(p,v^2)`$ has to be written as:

$$\frac{V}{\overline{n}}=12i\overline{n}\mathrm{exp}(i\omega (n))\varphi _{00}(v^2)+\mathrm{\hspace{0.17em}12}i\mathrm{exp}(i\omega (n))\varphi _{01}(v^2).$$

We can obtain this formula by changing variable in the integrand $`\rho =nr`$ and differentiating the obtained expression of $`V(n,v^2)`$ with respect to $`\overline{n}`$ and, finally changing the variable back: $`r=n\rho `$.

Partial differentiating of the double integral over $`\overline{n}`$ and $`n`$ gives:

$$_{\text{}}\frac{dnd\overline{n}}{pn}V(n,v^2)=\underset{R\mathrm{}}{lim}_{|n|=R}\frac{dn\overline{n}V(n,v^2)}{pn}+$$

$$+\underset{\epsilon 0}{lim}_{|pn|=\epsilon }\frac{dn\overline{n}}{pn}V(n,v^2)_{\text{}}\frac{dnd\overline{n}}{pn}\overline{n}\frac{V(n,v^2)}{\overline{n}}.$$
(103)
Asymptotic behaviour of $`V(n,v^2)`$ as $`|n|\mathrm{}`$ was obtained above. Using this asymptotics for the first part of the formula (103) gives zero. Second term of the right-hand side of (103) gives the residue of the integrand at $`p`$ multiplied by $`2i\pi `$. Then we obtain:

$$J_1(p,v^2)=\overline{p}V(p,v^2)\varphi _{00}(v^2)\mathrm{exp}(i\omega (p))$$

$$\mathrm{\hspace{0.17em}12}i\varphi _{01}(v^2)\frac{1}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{pr}\overline{n}\mathrm{exp}(i\omega (n))+$$

$$+iv^2\varphi _{00}(v^2)\frac{1}{2i\pi }_{\text{}}\frac{dnd\overline{n}}{pr}\mathrm{exp}(i\omega (n)).$$

On the same way one can evaluate an integral

$$J_1^+=_\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\mathrm{exp}(i\omega (n))_{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{nm}}\mathrm{exp}(i\omega (m)).$$

The difference between the integral $`J_1`$ and the integral $`J_1^+`$ consists in the addition terms over boundaries of the domains $`\mathrm{\Omega }^+`$ and $`\mathrm{\Omega }^{}`$ in the result. The finally forms are:
when $`p\mathrm{\Omega }^+`$:

$`J_1^+=\varphi _{00}^+(v^2)\mathrm{exp}(i\omega (p))+`$
$`+\overline{p}\mathrm{exp}(i\omega (p)){\displaystyle _{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{\overline{pm}}\mathrm{exp}(i\omega (m))}`$
$`{\displaystyle \frac{1}{2}}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\overline{n}\mathrm{exp}(i\omega (n))}+`$
$`+iv^2\varphi _{00}^+(v^2){\displaystyle \frac{1}{2i\pi }}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\mathrm{exp}(i\omega (n))}`$
$`12i\varphi _{01}^+(v^2){\displaystyle \frac{1}{2i\pi }}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\overline{n}\mathrm{exp}(i\omega (n))};`$
when $`p\mathrm{\Omega }^+`$:

$$J_1^+=\varphi _{00}^+(v^2)+p_{\mathrm{\Omega }^+}\frac{dmd\overline{m}}{p+\overline{m}}\mathrm{exp}(i\omega (m))$$

$`{\displaystyle \frac{1}{2}}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\overline{n}\mathrm{exp}(i\omega (n))}`$
$`iv^2\varphi _{00}^+(v^2){\displaystyle \frac{1}{2i\pi }}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\mathrm{exp}(i\omega (n))}`$
$`12i\varphi _{01}^+(v^2){\displaystyle \frac{1}{2i\pi }}{\displaystyle _\mathrm{\Omega }^{}\frac{dnd\overline{n}}{pn}\overline{n}\mathrm{exp}(i\omega (n))}.`$
Acknowledgements.

I am grateful to S.G. Glebov, A.M. Il’in, L.A. Kalyakin and M.M. Shakir’yanov for stimulated discussions.