Bessel function type solutions of the ultradiscrete Painlevé III equation with parity variables

The discrete Painlevé III equation ($$\mathrm {dP_{III}}$$dPIII) possesses a class of special solutions with determinantal structure whose entries are written by the q-difference Bessel function. In this paper, an ultradiscrete analogue of the q-Bessel function is constructed by ultradiscretization with parity variables (p-ultradiscretization). Based on this result, special solutions for the p-ultradiscrete Painlevé III equation are derived from those of $$\mathrm {dP_{III}}$$dPIII. The ultradiscrete solutions capture oscillating behaviour of the (differential) Painlevé III equation.


Introduction
Ultradiscretization [23] is a limiting procedure for reducing a given difference equation to a piecewise linear equation that is written by addition, subtraction, and the max operation among dependent variables. In this procedure, we first transform a dependent variable x n in the given difference equation by B Shin Isojima isojima@hosei.ac.jp 1 Department of Industrial and Systems Engineering, Hosei University, 3-7-2 Kajino-cho, Koganei, Tokyo 184-8584, Japan where ε > 0 is a parameter. Then, we apply ε log to both sides of the equation and take the limit ε → +0. By using the exponential laws and the identity lim ε→+0 ε log e X ε + e Y ε = max(X, Y ), (2) multiplication, division, and addition among x n are replaced by addition, subtraction, and the max operation among X n , respectively. The resulting piecewise linear equation, which is called an ultradiscrete system, is often considered to be the time evolution rule of a cellular automaton [24] since its dependent variables as well as its independent variables will take discrete values. We can regard the transformation (1) as an approximation of the solution x n by its leading term e X n /ε . In this context, the ultradiscrete system is the balancing relation among the exponents of leading terms, say X n . In spite of this drastic approximation, an ultradiscrete system may inherit essential properties of the original equation. In particular, if we have an exact solution of the original equation, we can constract a solution of the corresponding ultradiscrete system. The box-and-ball system (BBS) [20], which is an ultradiscrete analogue of the Korteweg-de Vries (KdV) equation, is a good example to observe this feature of ultradiscrete systems. The BBS has an infinite number of preserving quantities and therefore keeps the integrability of the KdV equation. Moreover, the BBS admits the N -soliton solutions, which illustrate soliton interaction in a simple manner. The direct relationship between the soliton solutions of BBS and those of the KdV equation is also clarified through the limiting procedure. Because of this interesting feature, ultradiscrete systems are faithfully studied from viewpoints of mathematical interest and applied possibility (See, for example, [19]).
However, x n must be positive to apply (1). Moreover, the equation that we consider must be subtraction-free since it is not straightforward to find a meaningful limit of ε log e X/ε − e Y /ε . These conditions strongly restrict a class of equations to which we can apply ultradiscretization procedure. For example, it is difficult to ultradiscretize a equation describing oscillation. Some attempts have been made to overcome this difficulty. Instead of (1), the ansatz x n = sinh(X n /ε) has been studied in [7,10]. Kasman and Lafortune investigated the limit of the form ε log γ i e w i (ε)/ε (γ i ∈ C) in detail [6]. An algebraic approach by means of a non-archimedean valuation was presented by Ormerod and it was applied to the discrete Painlevé III equation [16]. The present author proposed a new procedure called ultradiscretization with parity variables (p-ultradiscretization) with coworkers [14]. This procedure keeps track of the sign of the original variables by introducing the parity (sign) variable ξ n = x n /|x n |. We can study ultradiscrete analogues of a wider class of equations by virtue of these attempts.
The Painlevé equations are important nonlinear differential equations that have rich mathematical structures, for example, they admit a class of special solutions written in terms of a determinant whose entries are given by a hypergeometric-type function. Their q-difference analogues [17] are simply called the q-Painlevé equations. They also have rich structures and are actively studied (See, for example, [4,18]). Ultradiscrete analogues of the q-Painlevé equations and their special solutions are also investigated vigourously. Ultradiscrete Painlevé equations by traditional procedure are studied in [13,15,21], and so on. In [16], as mentioned above, hypergeometric solutions of the ultradiscrete Painlevé III equation and an ultradiscrete Bessel function were constructed by an algebraic approach. The procedure of p-ultradiscretization has been applied to the q-Painlevé II equation [3] in [9] and some successive works have been reported in [8,11,12]. Furthermore, the p-ultradiscrete analogue of the Painlevé VI equation and its special solutions have been proposed and the asymptotic behaviour of the solutions was discussed in [22]. It is an interesting problem to construct p-ultradiscrete analogues of other Painlevé equations and study their mathematical structures including special solutions.
In this study, we focus on a discrete analogue of the Painlevé III equation (dP III ), which is written as [5] where α, β, γ , δ, and λ are parameters. For a special set of parameters (3) admits a class of special solutions as follows. A q-difference analogue of the Bessel equation is written as where q is a multiplicative difference interval that satisfies |q| < 1. We write and consider a function with determinantal structure, Then, functions defined by solve (3) with parameters (4). Note that (7) satisfies the bilinear equations Since (3) is subtraction-free, it is ultradiscretizable in the traditional procedure. However, the above special solutions may not be captured directly, because (5) shows oscillating behaviour and the determinant (7) itself includes subtraction. Hence, it is reasonable to introduce p-ultradiscretization. This is another approach to ultradiscretization of the Painlevé III equation by means of a different method from [16]. Our aim is to derive a class of special solutions for p-ultradiscrete analogues of (9)-(12) from (7). Moreover, we construct some solutions of p-ultradiscrete analogues of (3) and study their behaviour. This paper is organized as follows. In Sect. 2, we review the procedure of p-ultradiscretization. In Sect. 3, we evaluate a q-difference analogue of the Bessel function as preparation for the main content. In Sect. 4, we evaluate (7) into which we substitute the q-difference Bessel function. In Sect. 5, we construct an ultradiscrete solution for the p-ultradiscrete analogue of (3) and compare its behaviour with the corresponding special solution of the Painlevé III (differential) equation. Finally, concluding remarks are given in Sect. 6.

Ultradiscretization with parity variables
In this section, we review the procedure of p-ultradiscretization. As an example, we consider a simple difference equation Note that (13) possesses the general solution One can formally obtain an usual ultradiscrete analogue of (13) by assuming x n > 0. If we deform (13) as x n+1 + ax n = b, all terms become positive. Hence, by putting a = exp(A/ε), b = exp(B/ε) and x n = exp(X n /ε), the deformed equation is ultradiscretized as However, we find no solution of (15) for some specific values of A, B and X n . For example, when A = 1, B = 0 and X n = 1, (15) does not satisfied for any values of X n+1 . In this meaning, (15) is not well-defined. This failure of ultradiscretization trivially comes from unreasonable assumption x n > 0. Introduction of p-ultradiscretization improves this situation. We assume x n = 0 and introduce the parity (or sign) variable ξ n ∈ {1, −1} by which gives the sign number of x n . We consider the following four cases for (13): (i) ξ n+1 = 1 and ξ n = 1, (ii) ξ n+1 = 1 and ξ n = −1, (iii) ξ n+1 = −1 and ξ n = 1, (iv) ξ n+1 = −1 and ξ n = −1. For each case, we deform (13) as Since all terms are positive in (17), each equation is ultradiscretized by introducing the amplitude variable X n by |x n | = exp(X n /ε) and the formal replacement 0 = exp(−∞/ε) as respectively. Now, we regard the obtained set of equations (18) as an ultradiscrete analogue of (13). A pair (ξ n , X n ) includes richer information of x n than usual ultradiscrete variable. Let us study (18) with A = 1, B = 0 and (ξ n , X n ) = (1, 1). For these values, the cases (ii) and (iv) do not apply to ξ n = 1 and the case (i) gives no solution. However, the case (iii) admits the unique solution X n+1 = 2 and this means ξ n+1 should be −1. Hence, we have found the solution (ξ n+1 , X n+1 ) = (−1, 2), which we could not for (15). We comment on indeterminacy in solving a p-ultradiscrete equation under specific conditions. Let us study (18) with A = 1, B = 2 and (ξ n , X n ) = (1, 1). Although the cases (ii) and (iv) do not apply again, the other two cases give the same equation max(X n+1 , 2) = 2. From this equation, we must have the result that the parity is arbitrary (ξ n+1 = ±1) and the amplitude is not unique (X n+1 ≤ 2). We proceed to the next step by choosing one of the indeterminate values. Therefore, we have an infinite number of solutions in this case. We note that a guideline to choose a 'nice value', whose meaning depends on property of the solution which one tries to capture, is still under study. Through this example, we find that uniqueness of solution is lost for some specific values of A, B and X n but a solution always exists. Although this indeterminacy of solution looks troublesome, the situation becomes better comparing from (15), in which a solution may not exist. Noticing the existence of solution, we claim that (18) is well-defined. This example also shows that p-ultradiscretization does not make the max operation invertible. Hence, a p-ultradiscrete system is not equivalent to the original equation but its approximation.
In order to give 'simpler' representation of (18), we use the following notation. We define a function Note that the sign ξ is represented as ξ = s(ξ ) − s(−ξ). Then, we rewrite x n as Xn ε (20) and substitute it into (13). In other words, we employ (20) instead of (1). Then, transposing the negative terms to the other side of the equation, we have We apply ε log to both sides of (21) and then take the limit ε → +0. If we define a function S by the following identity holds: Here, the term including −∞ vanishes from max. We may regard S(ξ ) as the formal ultradiscrete analogue of s(ξ ) with s(ξ ) = e S(ξ ) ε . By utilizing this identity, (21) reduces to max(S(ξ n+1 ) + X n+1 , S(−ξ n ) + X n + A) = max(S(−ξ n+1 ) + X n+1 , S(ξ n ) + X n + A, B).
By considering all cases for the parity variables, we recover the explicit equations (18) from the implicit equation (24). That is, we represent the four cases (i)-(iv) in terms of the function S. We call (24) as well as (18) a p-ultradiscrete analogue of (13). We hereafter present p-ultradiscrete equations in the implicit expression by the following reasons. Firstly, the implicit expression has shorter form than its explicit counterpart. Secondly, one can construct both of forward and backward schemes from the implicit expression, if necessary. We comment that it is not always reasonable to treat the implicit expression as an equation on max-plus algebra for characters X and S(ξ ).
is not well-defined for ξ = −1, which directly corresponds to zero-dividing for the original variable. We should treat (24) as equations on max-plus algebra with cases represented by S(ξ ), not on a new algebra.
In closing this section, we review the following lemma [12], which provides a useful sufficient condition for obtaining a solution of a p-ultradiscrete equation.

Lemma 1 If a solution of a given difference equation x n (ε)
, where ε is an arbitrary positive parameter, is evaluated as the pair of the sign (−1)ξ n and the amplitude X n solve the corresponding pultradiscrete equation.
We may simply refer to the pair ((−1)ξ n , X n ) obtained by this lemma as p-ultradiscrete limit of x n . For example, we consider (14) and assume n ≥ 1, A > 0, x 0 = e X 0 /ε > 0 and X 0 > B − A for simplicity. Then, we find that ((−1) n , n A + X 0 ) is the pultradiscrete limit of (14).

q-Bessel function and its ultradiscrete limit
In this section, we study a special solution of (5) with ν ∈ Z ≥0 , which we refer to as the q-Bessel function. It is given by where Jackson's q-Bessel function, which is defined even for ν ∈ C, is well known, where r φ s is the basic hypergeometric series The relationship between (26) and (28) for ν ∈ Z ≥0 is readily found.

Evaluation of the q-Bessel function
Our aim in this subsection is to prove the following proposition.

Proposition 1 (i) The q-Bessel function J ν (x) is deformed as follows:
(ii) The evaluation as q → 0 holds for n ∈ Z.
To prove Proposition 1 (i), we introduce the well known formula [2] 1 and the following lemma.
(ii) The formulaD Proof of Proposition 1 (i) We consider a function Note that If we put we haveĴ from (33). Hence, we obtainB Next, we proveD We introduce an identity which appears in [2], Exercises 3.2 (iii). Note that The left hand side of (49) is deformed as follows: The right hand side of (49) is rewritten as follows: We have (48) by comparing these two series. Now, applying the q-Euler transformation,Ĵ ν is deformed as follows: which gives by (44). Moreover, (54) reduces to (31) using (30).
Proof of Proposition 1 (ii) Substituting x = q n into (26) and (31), we obtain respectively, wherẽ Note that If n ≥ 1, the term with j = 0 in (55) is dominant. Hence, we readily obtain Since it is difficult to evaluate (55) for n ∈ Z ≤0 , we focus on the other expression (56) and study the minimum value off ν,n (k) for k ∈ Z ≥0 . If 0 ≥ n ≥ −ν, we find that minf ν,n (k) =f ν,n (0) = 0 and that the term with k = 0 in the summation in (56) is dominant. Hence, we obtain Next, if n < −ν and ν + n is even, we find that Hence, the term with k = − ν+n 2 is dominant and then holds. Finally, we consider the case where n < −ν and ν + n is odd. We put ν + n = 2μ Then, the leading term in the summation in (56) may be given by the sum of the terms with k = −μ and −μ − 1, We also find from |f ν,n (−μ)−f ν,n (−μ+1)| = 4 that the other terms do not contribute to the leading term. Therefore, we have Now, (32) has been proved.

Ultradiscrete solution and supplementary result
If we put x = q n , q = e Q ε (Q < 0), and Applying Lemma 1 to (32), we readily obtain an explicit expression for the pultradiscrete analogue of the q-Bessel function Note that this result is identical to the ultradiscrete Bessel function by Narasaki (see Appendix 1), unless n ≤ −ν − 1 and n + ν is odd. Another supplementary result can be obtained from Proposition 1. It is concerned with the number of restricted partitions p n (k) defined by the generating function We have the series expression of (43) by utilizing (70), From (53), we also obtain Comparing these two expressions, we find that (71) has an infinite number of extra terms for a negative integer n. From this fact, we have the following formula.
Proposition 2 For given ν ∈ Z ≥0 and n ∈ Z <0 , fix m ∈ Z such that holds, where j and k run over all pairs satisfying j + kn = m in the summation.
We do not know a combinatorial meaning of this (formal) formula yet.

Evaluation for special solutions of the discrete Painlevé III equation
In this section, we substitute the q-Bessel function into (7) and study the p-ultradiscrete analogue of the resulting function. If τ ν N (n) is written in the form of (25), its pultradiscrete analogue is readily obtained by Lemma 1. We first explain a useful notation and then summarize the main result.
We rewrite (32) in a simpler form since we often use it. We introduce ψ(n) := n(n + ν − 1) (75) and the binomial coefficient n 2 = (n − 1)n/2. Using (6) and these notations, (32) can be rewritten as follows: For τ ν N (n), we define the sign y ν N ,n and the amplitude Y ν N ,n by respectively. Moreover, we use the formula [22] to obtain simpler expression for a multiplicative term. Then, p-ultradiscrete analogues of (9)-(12) are respectively written as These p-ultradiscrete equations are solved by y and Y , as presented below, which are obtained as the p-ultradiscrete analogue of (7). To give explicit functional forms, we introduce the following six cases: Note that there exist values of n such that all antidiagonal elements in (7) are of the type (80) for ν ≥ N − 1, but such values of n do not exist for ν ≤ N − 2. We also introduce where floor(x) denotes the integer part of x. Now, the p-ultradiscrete analogue of (7) is written as follows. Case (A): Case (B): Case (C-a): Case (C-b): Case (D): Case (E): We study each case in the following subsections.

Case (A)
In this case, all arguments of J ν (n) in (7) are positive. Substituting (55) into (7) and using multi-linearity of the determinant, we obtain where Since the determinant in (103) is the Vandermonde determinant, the terms with j k = j l (k = l) disappear. Moreover, under any permutation of ( j 1 , j 2 , . . . , j N ), the absolute values of all elements but the factor N k=1 q 2(n+2k−2) j k in (103) are invariant. Therefore, noting 0 < q < 1, the largest absolute value of the monomial in (103) is achieved by the j k 's with 0 ≤ j N < · · · < j 2 < j 1 , which maximize N k=1 q 2(n+2k−2) j k . Then, in the Vandermonde determinant, the product of diagonal elements N k=1 q 2(k−1) j k contributes to the evaluation. Hence, we study j k 's that maximize N k=1 q 2(n+2k−2) j k +2(k−1) j k ; i.e., j k 's that minimize under (106). We readily find that such j k 's are given by Moreover, further contributions by the factor N k=1 q (n+3k−3)ν and sign N i=1 (−1) N −i = (−1) ( N 2 ) from P ν N ( j ) must be considered. Accordingly, the leading term of τ ν N (n) is given by

Cases (C-a), (D), and (E)
In these cases, all arguments of J ν (n) in anti-diagonal elements of (7) are negative.
The following proposition plays an important role in the evaluation.

Proposition 3
Write the absolute value of the leading term of J ν (n) as q g ν (n) . For n ∈ Z and k, l ∈ Z >0 , holds.
Moreover, this is true even for the case in which all anti-diagonal elements have nonpositive arguments. We first study the case (E). In this case, all anti-diagonal elements of (7) are of the type (81), and therefore we have where we have used the notation (76)-(78). For the case (C-a), all anti-diagonal elements in (7) are of type (80). The leading term of (7) is given by where ψ is defined by (75). The case (D) gives a 'mixed' situation. The anti-diagonal elements in (7) are of type (80) or (81). Using (75)-(78), the leading term can be represented as follows: (114)

Case (B)
In this case, J ν (n)'s of types (79) and (80) appear as elements in (7). Moreover, there exist some rows (or a row) that comprise only J ν (n)'s with positive arguments. The product of anti-diagonal elements is no longer dominant. We consider M defined by (90). Then, the ith row (M + 1 ≤ i ≤ N ) comprises only J ν (n)'s with positive arguments. When n is even, (7) has the form .
We present the procedure for calculating its leading term. Although we illustrate the procedure only for the case in which some anti-diagonal elements are of type (80), this procedure applies, in a similar manner, for the case in which all anti-diagonal elements are of type (79) or for odd n. We use the notation {n} := (−(1 − q) 2 q 2n ; q 2 ) ∞ for simplicity. Substituting the series expressions of J ν (q n ) and employing multi-linearity of the determinant, we obtain where we have used (104) and (105). Let the (κ, λ)-element for the determinant in (116) be of type (80). It has the form q 2 j κ ( j κ +ν+n+2κ+λ−3) /{n + 2κ + λ − 3} and ν + n + 2κ + λ − 3 ≥ 0 (see Sect. 2). Hence, the minimum of j ( j + ν + n + 2κ + λ − 3) occurs at j = 0. Note that the determinants with j k = j l (1 ≤ k, l ≤ M) do not disappear, unlike those in case (A). Therefore, the largest absolute value of the monomial in this determinant occurs with the j k 's given by Although other sets of j k 's actually produce the largest absolute value, we discuss them later. Noting that (1 − q) k ∼ 1 and P ν N ( j ) ∼ (−1) 2 ) under (117), our aim turns into evaluating . . .
We deformτ ν N (n) according to the Procedure presented in Appendix 1. To represent the resulting expression, we introduce additional notation. For x k = q 2(N −M−k) , k = 1, 2, . . . , N − M − 1, we consider the fundamental symmetric expressionã k among them, i.e., Note thatã as q → 0. Furthermore, we add a new variable x N −M = 1 to x 1 , . . . , x N −M−1 and write the fundamental symmetric expression among them as a k (k = 1, 2, . . . , N − M) and set a 0 = 1. Then, the relations hold. For convenience, we extend the definition of n in (120) to Using this notation, we define a function Then, after applying the Procedure,τ ν N (n) is written as follows: where, for simplicity, we have used the trivial formula (−1) 2(N −M) = 1. Since (q 2 ; q 2 ) k ∼ 1, our aim reduces to calculating the leading term of det(T (n + 2κ + λ − 3)), which, from the following lemma, is given by the product of anti-diagonal elements.
Therefore, the leading term of det (T (n + 2κ + λ − 3)) is written as follows: and then, the leading term ofτ ν n (n) is obtained from (125). However, when n + N + M − 2 ≥ 1, that is, when T (n) with n ≥ 2 appears in the anti-diagonal elements of det (T (n + 2κ + λ − 3)), the leading term of τ ν n (n) is not obtained only fromτ ν n (n). We must consider other determinants that have the same order of leading terms, in other words, other sets of j k 's. We present the following lemma:

Lemma 4
Suppose that an anti-diagonal element in det (T (n + 2κ + λ − 3)) is T (ñ) withñ ≥ 2 and it is in the ith row ( 1 ≤ i ≤ M). In (116), consider the determinants given by the sets of arguments The order of leading terms of these determinants is identical to that obtained from (117).
From m ≤ñ − 1 and k ≤ñ + 1, we readily obtain m − k + 1 < 0 and therefore which is independent of m and moreover of i.
If we change the original set of j k 's (117) into that of (131), an extra sign (−1) j i arises from P ν N ( j ). Hence, forñ ≥ 2, we may replace (129) with We accordingly obtain ) . (136)

Case (C-b)
In this case, all types of J ν (n)'s appear as anti-diagonal elements in (7). Although our discussion is similar to that for case (B), elements of type (81) contribute to the dominant term. The result is represented as follows:

Behaviour of the special solution
The Painlevé III equation (P III ) is written as follows: where a, b, c, d are parameters. Equation (138) can be derived from (3) through the following continuous limit. That is, (3) reduces to (138) if we introduce a parameter < 0 and put u(n) = λ −n/2 w(n), n = z, α = −1/c 2 , β = −d/c, γ = −b/c, δ = a/c, and λ = 1+2 ; then we take the limit → −0 and transform the independent variable by z = e x . For (138) admits special solutions written in terms of a solution of the Bessel equation That is, the function u ν N (x) defined by solves (138) with (139), where v ν (x) is a solution of (140). We consider (3) with (4) and construct its p-ultradiscrete analogue by introducing w ν N (n) = ω n e W n /ε and q = e Q/ε . Applying the procedure of p-ultradiscretization, we obtain (the special case of) the p-ultradiscrete Painlevé III equation (p-uP III ) max S(ω n+1 ω n ω n−1 ) + W n+1 + W n + W n−1 + (ν + N )Q, Moreover, the p-ultradiscrete analogue of (8) is written as follows: We regard this relation as the equation for ω n and W n and solve this equation. Then for almost all cases, we obtain the unique transformation but for the special case in which Y = (ν + N )Q and y = 1, we obtain the indeterminate transformation Employing this transformation, we construct a special solution to (144) from (y ν N ,n , Y ν N ,n ) obtained in Section 4. Here, we study the case in which N = 1 and ν = 1. For n ≤ −3, we uniquely obtain For n ≥ −2, we encounter the indeterminate transformation Therefore, we substitute Though we omit details, we find that ω −2 is arbitrary but (ω −1 , W −1 ) = (1, 4Q) is required to obtain such (ω 0 , W 0 ). Continuing these experiments, we guess a special solution We can check that this function satisfies (144) by direct substitution. Now, we have obtained a special solution given by (151) and (155). From this example, we find that (152) does not give a unique solution to (144), but it does give important information about the solution, which we want to consider. We study the behaviour of this ultradiscrete solution by introducing the variable ω n e W n . This variable shows qualitative behavior similar to that of w(n) since it is obtained from w ν N (n) = ω n e W n /ε with ε = 1. Moreover, it corresponds to the continuous variable e −x u ν N (e −x ), where u ν N (x) is defined by (141). We find this correspondence from variable transformations in the continuous limit. Figure 1 shows a plot of ω n e W n . Its behaviour is similar to that shown in Fig. 2 for e −x u ν N (e −x ). That is, both solutions tend to zero as n or x → ∞, and both show oscillating behaviour for negative values of the independent variables. Since only discrete points are considered in the ultradiscrete equation and scaling by ε is ignored, the divergence in Fig. 2 is unclear in Fig. 1.
the evaluation has become non-trivial because of the special independent variable x = q n . The key to the evaluation for n ≤ 0 is to deform the known series expression for J ν (q n ) by the q-difference Euler transformation.
Based on the evaluation for J ν (q n ), we have derived the explicit functional forms of special solutions for the p-ultradiscrete analogues of the bilinear equations for dP III , for general values of system parameters N ∈ Z ≥1 and ν ∈ Z ≥0 . These special solutions are obtained by evaluating those of dP III represented by the Casorati determinants whose elements are given by the q-Bessel function. The evaluation for determinantal solutions of the q-Painlevé II equation was discussed in [12]. Compared with its calculation, we have needed a new technique to evaluate the determinant for dP III as shown in Sects. 4.3 and 4.4.
We have also constructed a special solution of p-uP III through the variable transformation. We have found that the behaviour of the obtained ultradiscrete solution is qualitatively similar to the corresponding continuous solution. Although the transformation locally takes the indeterminate form (150), it enables us to find a special solution in a fairly specified form by considering its global behaviour. Such application of indeterminate form has newly observed in this research.
It is a future problem to give ultradiscrete special solutions corresponding to the q-Neumann function and to the general solution as studied in [8]. Another future problem is to construct ultradiscrete analogues of other Painlevé equations and to investigate their mathematical structures, such as degeneration structures or Bäcklund transformations. We hope that the results in this paper will contribute to future analysis of discrete Painlevé equations.
This solution was constructed as a solution to the specific initial value problem for (68).

Appendix 2
We present the procedure used to deformτ ν N (n) in Sect. 4.3. Procedure: step 0-1 : Add the N th row × (−1) to the first row. step 0-2 : Add the N th row × (−1) to the second row.