Discrete harmonic functions in the three-quarter plane

In this article we are interested in finding positive discrete harmonic functions with Dirichlet conditions in three quadrants. Whereas planar lattice (random) walks in the quadrant have been well studied, the case of walks avoiding a quadrant has been developed lately. We extend the method in the quarter plane -- resolution of a functional equation via boundary value problem using a conformal mapping -- to the three-quarter plane applying the strategy of splitting the domain into two symmetric convex cones. We obtain a simple explicit expression for the algebraic generating function of harmonic functions associated to random walks avoiding a quadrant.


Introduction
In the two-dimensional continuous case, a harmonic function is a function for which the standard Laplacian ∆ = ∂ 2 ∂x 2 + ∂ 2 ∂y 2 is zero. For example, in the cone {r exp(it) : 0 r < ∞, 0 t η} of opening angle η, then u(r, t) = r π/η sin π η t is the unique (up to multiplicative constants) positive harmonic function equal to zero on the boundary.
In the two-dimensional discrete case, consider the simplest Laplacian operator The function f (i, j) = i is positive harmonic in the right half plane and equal to zero on y-axis.
In the positive quadrant, the function f (i, j) = ij is positive harmonic with Dirichlet boundary conditions. What about the three-quarter plane? Surprisingly, harmonic functions in the threequarter plan have more complex expressions, see (14) and (15). Unlike the two last examples, the three quadrants is not a convex cone and finding positive harmonic functions is more complicated. Although the difference between the intersection of two half-planes (the quarter plane) and union of two half-planes (the three-quarter plane) seems geometrically insignificant at first sight, Bousquet-Mélou has already noticed this complexification with a paths enumeration point of view: in [6] Bousquet-Mélou notes that studying the simple walks (with cardinal points jumps) is as difficult as studying the Gessel walks in the quadrant (a famous model with jumps {E, NE, W, SW}, solved after ten years of research). In this article, we develop a systematic approach to find positive discrete harmonic functions in the three-quarter plane with Dirichlet conditions.
Context. In the discrete case, planar lattice random walks in cones occupy a central position in probability and combinatorics. Harmonic functions play an important role in probability theory. Doob h-transform is a way to build conditioned random processes in cones from a random process and a positive harmonic function vanishing on the boundary of the cone. Finding positive harmonic functions for random processes is therefore a natural objective in the study of confined random walks. There are very few ways to compute discrete harmonic functions (see [17] and references therein).
Most of walks study have been done on the quadrant, or more generally in convex cones. A natural generalization is to consider other domains of restriction and determine how the framework of (random) walks is different from the quarter plane to this region. Recently, non-convex cones, in particular the three-quarter plane have been examined. Unlike the quarter plane, where generating functions involve positive powers of the variables, in the three-quarter plane we face both positive and negative powers of the variables, making convergence issues arise. A natural strategy consists of cutting the three quadrants in some convex cones in which the generating functions are convergent. In [6], Bousquet-Mélou sees the three-quarter plane as the union of three quadrants and obtains some results for the simple and diagonals walks avoiding a quadrant. Integral expressions for the generating function of walks avoiding a quadrant with symmetric step sets for walks are derived in [18], where the three-quarter plane is seen as the union of two symmetric convex cones of opening angle 3π/4. Asymptotics of the number of excursions of walks with small steps in the three-quadrant is computed in [15] by Mustapha. In this article, following [9,4], Mustapha expresses the critical exponent of harmonic functions in three quadrants as a function of the critical exponent of harmonic functions in a quadrant. When this exponent is not rational, then the generating function of walks is not D-finite and [15,Thm 1.3] proves that the generating function of the walks of the 51 non-singular step sets with infinite group are not D-finite in the three-quarter plane (recall that neither are they in the quadrant).
In this article, we find an explicit expression for generating functions of discrete harmonic functions associated to random walks avoiding a quadrant with a mixed approach of [17] and [18]. We focus on the analytic approach developed in [17], which consists of writing a functional equation for the generating function for a fixed harmonic function, transforming this functional equation into a boundary value problem and finally solving this problem, which results in an explicit expression for the generating function. We begin by making some assumptions on both the random walks and discrete harmonic functions for these random walks.
(H1) The walk is homogeneous inside the cone with transition probabilities {p i,j } −1 i,j 1 to the nearest neighbors; (H2) The transition probabilities are symmetric (p i,j = p j,i ) and p 0,0 = p −1,1 = p 1,−1 = 0; (H3) In the list p 1,1 , p 1,0 , p 0,−1 , p −1,−1 , p −1,0 , p 0,1 , there are no two consecutive zeros; (H4) The drifts are zero: We first suppose with (H1) the random walks to be homogeneous with small steps, which therefore provide us techniques and tools developed in [11]. Moreover, with assumption (H2), we suppose the walks to be symmetric with no anti-diagonal jumps. The third hypothesis (H3) is not essential in the study, but automatically excludes degenerate cases which could have been studied with easier methods. Finally, we assume the walks to have zero drifts with (H4). Note that the zero and non-zero drifts are two very different frameworks [10], and most results are given in the non-zero drift case. Furthermore, (H4) makes the random walks hit the negative axes almost surely. Figure 1 illustrates a possible transition probabilities set which satisfies (H1), (H2) and (H3). We will use the combinatorial step set terminology: the simple model is when p 1,0 = p 1,0 = p −1,0 = p 0,−1 = 1/4, Gouyou-Beauchamps model has p 1,0 = p −1,1 = p −1,0 = p 1,−1 = 1/4 and Gessel model is p 1,0 = p 1,1 = p −1,0 = p −1,−1 = 1/4. In addition, we ask the associated discrete harmonic functions f = (f (i, j)) (i,j)∈C to satisfy four properties: We also make the hypothesis (P4) which, as we shall see (Remark 5), will be automatically satisfied for positive harmonic functions associated to symmetric step sets.
The first property (P1) is the harmonicity condition, the second one (P2) is the zero condition on the boundary, the third one (P3) is the positivity condition within the cone and the last one (P4) is a symmetry condition (coming then from (H2)). In other words, we are interested in symmetric positive discrete harmonic functions for symmetric random walks with small steps constrained to the boundary of C. Thereafter, for such a fixed harmonic function f , we consider its generating function We observe here one of the main difficulties of the three quadrants: the series H(x, y) is not convergent. Thereafter, we may see this object as a formal series. The goal is to find a formal explicit expression of H(x, y), and from this expression, it is not difficult to get the f (i, j) for all (i, j) ∈ C, with a formal power series expansion.
Method. Working in non-convex cones, here in particular in the three-quarter plane, arises convergence problem. Indeed, it is not difficult to write a functional equation from the properties of harmonicity (P1) and vanishing conditions on the boundary (P2). However, even if this functional equation (13) seems closed to the quarter plane case (17) solved in [17], they are in fact fundamentally different: the three quadrants case involves negative powers of x and y making the series not convergent anymore. To remedy this difficulty, we follow the same strategy as in [18]: we divide the three quadrants into two symmetric convex cones and the diagonal (see Figure 5) and write a system of two functional equations (one for each cone). At first sight, this cut increases even more the level of difficulty: we have now two functional equations and more unknowns to deal with, but in the particular case of symmetric harmonic functions and symmetric transition probabilities, the system is composed of twice the same equation and the problem can be seen as a slightly different variation of the quadrant case (see Figure 6). In this symmetric case, we are able to use the tools and methods of [17]: transform the functional equation into a boundary value problem, solve it and write an explicit expression for the generating function. Finally, let us point out that this method of splitting the domain in two octants can also be applied in the quadrant in the symmetric case, and allow us to give alternative proofs of [17].
Main results. Our main result is an expression of the diagonal section with P is a polynomial in R[y], w ϕ a conformal mapping (see Section 3.1, w ϕ satisfies an algebraic differential equation and can be simply expressed in terms of sine and arcsine functions) and G ϕ an explicit function defined from w ϕ . In particular, taking P of degree 1, we get the unique positive discrete harmonic function (satisfying (P1), (P2) and (P3)) and the generating function D(x, y) has exactly the same asymptotics that the ones found in the article of Mustapha [15]. After some simplifications, we have (unique positive harmonic function, P is of degree 1), (3b) with C and C ′ positive constants, W ϕ a conformal mapping and δ ϕ a polynomial. This expression suffices to give a complete solution to the problem and Theorem 15 gives a formal explicit expression for H(x, y). This formal series is expressed with the kernel K(x, y), a polynomial in x and y, and two conformal mappings W ϕ (x) and W ϕ (y). Each quantity depends on the transition probabilities of the random walk.  Figure 2: A random walk with associated angle θ is the three quadrants (resp. the quadrant), can be seen as an inhomogeneous random walk in the half-plane with associated angle π − θ 2 (resp. θ 2 ) in each quadrant.
As mentioned earlier, for the simple step set, the positive discrete harmonic function is defined by f (i, j) = i in a half plane and by f (i, j) = ij in the positive quadrant. We have then a product of degree one in the half plane and of degree two in the quadrant. Along our study, we find similar relations on harmonic functions in various quadrants (see Lemma 11 and Equation (66)). The growth of harmonic functions can be expressed from the angle θ defined by [12, Sec. 1.2] In the three-quarter plane, the critical exponent of the harmonic function is π/(2π − θ) = π/(2θ ϕ ).
After the decomposition of the three quadrants and the changes of variables, the random walk in the three quadrants is equivalent to an inhomogeneous random walk in the half-plane (see Figure 10), with associated angle θ ϕ = π − θ/2 in each quadrant (see Figure 2a). We also recover the factor two between harmonic functions in the quadrant and harmonic functions in the half plane: in the quarter plane, the critical exponent is π/θ = π/(2θ ψ ). With the same reasoning of the three quadrants case, a random walk associated with angle θ can be seen as an inhomogeneous random walk in the half-plane, with associated angle θ ψ = θ/2 in each quadrant (see Figure 2b).
Structure of the paper. In Section 2, we start with underlining an important tool used in this analytic study of harmonic functions of random walks with small steps. This object, denominated by K(x, y), is called the kernel of the walks and is a polynomial of degree 2 in x and y which encodes the transition probabilities of the walk. Further, we set up various functional equations for generating series of harmonic functions involving K(x, y).
In Section 3, we express explicitly generating functions of harmonic functions for some models of random walks. To come to these results, we transform a functional equation (built in the previous section) into a boundary value problem in Subsection 3.2, and the solution process of this problem can be found in Subsection 3.3. The resolution of the boundary value problem involves conformal gluing functions, introduced in Subsection 3.1. We end this part with the application of the results to the example of the simple random walks (Subsection 3.4).
In Appendix A, we suppose that we don't have the symmetry of the transition probabilities of the walks. We are still able to build kernel functional equations, but unfortunately there are too many unknowns and we are not able to solve this difficult problem yet.
In Appendix B, we apply this decomposition of the domain into two convex cones to symmetric harmonic functions of symmetric random walks in the quadrant and find partially the same results as [17]. We invite the reader to refer to this appendix along the consultation of the main body on the three quadrants.

Kernel of the random walks
The kernel of the random walks, denoted by K(x, y), is an important object of the study. It appears in the various functional equations of this section, and we use it later to transform a functional equation into a boundary value problem. The kernel is characterized by the transition probabilities of the associated random walk by The kernel is a polynomial of degree two in x and y (consequence of (H1) and (H3)) and can be written as where We also define the discriminants in the x-plane and in the y-plane: The discriminant δ(x) (resp. δ(y)) is a polynomial of degree three or four (this is a consequence of (H2)). Hence there are four branch points x 1 , x 2 , x 3 , x 4 (resp. y 1 , y 2 , y 3 , y 4 ), with x 4 = ∞ (resp. y 4 = ∞) when δ(x) (resp. δ(y)) is of degree 3. Examples of kernel coefficients are given in Section 3.4.
We follow the notations of [17] which are slightly different from the usual ones in [11]. Indeed, [17] points out that the kernel in (5) is the reciprocal one of the kernel Q(x, y) in [11]: K(x, y) = x 2 y 2 Q(1/x, 1/y).
Lemma 1 (Sec. 2.3 in [11] and Sec. 2.5 in [17]). The branch points x i are real and More precisely, we have x 1 ∈ [−1, 1) and Symmetric results hold for the y i , roots of δ(y).
Lemma 2 (Lem. 6.5.1 in [11]). The curve L in (12) is symmetric with respect to the horizontal axis, smooth except at Y (1) = 1 where it may have a corner point. At this point, the angle between the curve and the segment [y 1 , 1] is given by by (4).
is the coefficient of correlation of the increments of the walk.

A first functional equation
From the properties (P1) and (P2), we deduce a functional equation satisfied by the generating function H(x, y) defined in (2) for random walks under the hypothesis (H1) where This equation looks similar to the functional equation for discrete harmonic functions of random walks in a quadrant written in [17,Sec. 2]. In this case, we assume that random walks satisfy the hypotheses (H1), (H2), (H3), (H4). We also ask the associated discrete harmonic functions to have the properties (P1), The generating function of such harmonic functions f satisfies the equation L(x, y) As noted above, the structure of equations (13) and (17) are the same except that in the case of (13) there are infinitely many terms with positive and negative powers of x and y. As noticed in [18], this difference is not anecdotal, the series are not convergent anymore (for instance, it is not possible to evaluate them) and almost all the methodology of [7,16] can no longer be performed.
This difficulty appears as well in the study of counting walks avoiding a quadrant and in order to avoid this situation, in [6] Bousquet-Mélou views the three-quarter plane as the union of three quadrants and with a combinatorial approach, gets results on the simple and the diagonal walks. In [18], the three quadrants are split in two symmetric convex cones of opening angle 3π 4 . In this article, we follow the same strategy as in [18].

Functional equations in the 3π 4 -cones
The cone C is cut into three parts: the lower part {i 1, j i − 1}, the diagonal {i = j} and the upper part {j 1, i j − 1}. Let L(x, y) (resp. D(x, y) and U (x, y)) be the generating function of harmonic functions evaluated in the lower part (resp. diagonal and upper part), see Figure 5. By construction, we have where Lemma 4. For any random walks with properties (H1) and (H2), the generating function L(x, y) satisfies the following functional equation Remark 5. The symmetry of both the random walks (H2) and the harmonic function (P4) are crucial.
In fact, we can directly observe the symmetry of the positive discrete harmonic function in the case of a symmetric random walk. Let S(n) be a symmetric random walk in the three quadrants and τ (i,j) the exit time of the random walk starting at (i, j) from C. Writing (i, j) = (r cos(t), r sin(t)), the harmonic function for the Brownian motion in the cone {r exp(it) : 0 After the change of variable (i, j) → (−j, i) = r cos t + π 2 , r sin t + π 2 , the harmonic function for the Brownian motion avoiding the negative quadrant is u(r, t) = r 2/3 sin( 2 3 t + π 2 ), t ∈ − π 2 , π . This harmonic function is symmetric with respect to the diagonal (u r, π 4 − t = u r, π 4 + t ). There is a direct link between the positive discrete harmonic function V and the harmonic function u for the Brownian motion (see [9,Lem. 12]): The exit time τ (i,j) and the harmonic function u are both symmetric in the symmetric cone C. Therefore the positive discrete harmonic function V is symmetric as well.
Proof. Thanks to (P1) and (P2), we can easily write a functional equation for L(x, y) and D(x, y). Let D u (x, y) (resp. D ℓ (x, y)) be the generating function of a harmonic function evaluated in the upper (resp. lower) diagonal: We have Due to the symmetry of the cut and the random walks (p i,j = p j,i , hypothesis (H2)), we can simplify the last equation and get Plugging (23) into (21) and multiplying by xy, we get (20).
Equation (20) is transformed into and K ϕ (x, y) = α ϕ (x)y 2 + β ϕ (x)y + γ ϕ (x) = α ϕ (y)x 2 + β ϕ (y)x + γ ϕ (y), The functional equation (25) now has a closer structure to the functional equation in the quadrant (17) (see Figure 6), and to solve this problem, we can use some tools stated in [17]. The significant difference is the mixed factor in x and y in front of D ϕ (y), whereas in the quarter plane case, see Equation (17), the generating function H(x, y) is decomposed into a sum of two univariate functions: one section on the x-axis, namely K(x, 0) H(x, 0), and the other on the y-axis, namely K(0, y) H(0, y). This functional equation (25) holds at least for |x| < 1 and |y| < 1. Indeed, in [17, Sec. 2.1], it is proved that for positive discrete harmonic function f in the quarter plane, i 1 f (i, 1)x i−1 , j 1 f (1, j)y j−1 and further K(x, y) (i,j)∈Q f (i, j)x i−1 y j−1 is bounded at least on {(x, y) ∈ C 2 : |x| < 1, |y| < 1}. We can finally deduce that K ϕ (x, y)L ϕ (x, y), D ϕ (y) and L ϕ (x, 0) are bounded at least for |x| < 1 and |y| < 1 as well.
Remark 6. Note that in (H2), we exclude non-zero probabilities p −1,1 and p 1,−1 . Indeed, as noted in [18,Sec. 2], after the change of variable ϕ, we would end up with random walks with non-zero probabilities p −2,1 and p 2,1 respectively for big jumps of vector (−2, −1) and (2, 1). Note that the analytic theory for walks with big steps is not known.

Expression for the generating functions
We start this section with the introduction of an important conformal gluing function W (x) that maps the domain G M to the complex plane cut along a segment. Then we transform the functional equation (25) into a boundary value problem, which is a problem involving both regularity and boundary conditions. Finally we solve this problem and end up with an explicit expression for the generating function H(x, y). The boundary value problem satisfied by the generating function D ϕ (y) and ultimatly the generating function H(x, y) will be both expressed in terms of w. Further, all quantities with a ϕ-tag are objects defined from the changed random walks with kernel K ϕ (x, y), as in Section 2.3.

Conformal gluing function
Definition 7 (Conformal gluing functions). A function W is said to be a conformal gluing function for the set G M if: • W is continuous on G M \ {1} and W (1) = ∞.
Remark 8. In fact, to be a conformal gluing function, W only needs to be injective in G M , analytic on G M and to satisfy W (x) = W (x) for x ∈ M, as soon as W (x) is defined and finite. Moreover, W has an infinite value on M (or W would be constant, see [17,Lem. 4]). In this article, we choose W to be infinite at 1 (if W m is a conformal gluing function taking the value ∞ anywhere else on M, then we can consider W = 1/(W m − W m (1))).

The expression of such a function is computed in [12, Sec. 2.2]. Let
Then with θ defined in (4). Moreover, there exists c = 0 such that for x in the neighborhood of 1,  (29) is valid only for |T (x)| 1. When |T (x)| 1, writing arcsin(T ) = π/2 ± i ln(T + √ T 2 − 1), we can write Lemma 9. The function W defined in (29) has a zero at x 1 and is equal to 1 at X(y 1 ). Moreover, for x ∈ G M , it satisfies the following differential equation Proof. We start from (29). After differentiation and squaring, we can write The strategy now is to write δ(x) under its factorized form (remember that δ has four roots x 1 , 1, 1, and x 4 with possibly 2δ(x) (and δ ′′ (1) < 0). We also note that T (x 1 ) = 1 which implies that W (x 1 ) = 0 (and because W is injective, it is the only zero). The only zero of W ′ is at X(y 1 ) (see [3,Sec. 5.3]), and from (33), we deduce that W (X(y 1 )) = 1.
Similarly, we define W a conformal gluing for G L in the sense of Definition 7. We show that W and W are strongly related as Indeed, see [16, Thm 6 and Rem 6], W (Y (x)) is a conformal gluing for G M but may not be the one defined by (29). However, we know for sure that W (Y (x)) = aW (x) + b, with a, b ∈ C. By plugging in x = x 1 and x = X(y 1 ), we deduce a = −1 and b = 1. Thereafter, we are interested in a conformal gluing function which vanishes at 0. We consider w defined by The function w is a conformal gluing function in the sense of Definition 7, vanishes at 0, and satisfies the differential equation (which can be easily derived from (32)) Remark 10. Let w (resp. w ϕ and w ϕ ) be a conformal gluing function for G L (resp. G Mϕ and G Lϕ ) which vanishes at 0. Similar results and properties hold for those functions.
The conformal gluing functions w ϕ and w ϕ are defined from θ ϕ which can be computed with (4) and transition probabilities ϕ(p i,j ) = p i−j,i . The angles θ ϕ and θ are simply related (see Figure 2a), as stated in the following lemma. Relations between θ and θ ϕ will be ultimately interpreted as relations between the growths of harmonic functions.
In fact, in [16,Rem. 6], it is said that if W is a suitable mapping, then any aW (x)+b cW (x)+d with a, b, c, d ∈ C such that ad − bc = 0 is also a suitable mapping. In our case, we are interested in mapping with a pole at 1. As Y (1) = 1, by letting x goes to 1, we get c = 0. Hence we only need to find a and b such that W (Y (x)) = aW (x) + b. Lemma 11. The angle between the curve L ϕ and the segment [y ϕ,1 , 1] is given by with θ defined in (4).

Boundary value problem
Lemma 12. The generating function D ϕ (y) can be analytically continued from the unit disc D to the domain G Lϕ ∪ D and is continuous on G Lϕ \ {1}. Moreover, for all y ∈ L ϕ \ {1}, D ϕ (y) satisfies the following boundary condition where g ϕ is defined in (40) and w ϕ is a conformal gluing function defined in Section 3.1.
Remark 13. The boundary value problem in Lemma 12 is very close to the boundary value problem stated in [17,Lem. 3].
Proof. We first assume the analyticity and continuity of D ϕ (y) and we begin to prove the boundary condition (37). We start to evaluate the functional equation (25) at Y 0,ϕ (x) for x close to [x ϕ,1 , 1): By letting x go to any point of [x ϕ,1 , 1] with a positive (resp. negative) imaginary part, we obtain two new equations. We subtract these two equations and get: With (10) this last equation can be simplified as This boundary condition is the homogeneous equation of the boundary condition in [18,Lem. 5], with the same difficulty to deal with a non-meromorphic prefactor in G Lϕ , namely d(y) in (39).
Therefore we follow the same strategy as [18,Sec 3.4] and introduce the function g ϕ defined in G Lϕ by with h(y) = −y + y 2 − 1 and w ϕ defined in Section 3.1. The function g ϕ /g ′ ϕ is analytic in G Lϕ , has finite limits on L ϕ and the condition √ δϕ(y) √ δϕ(ȳ) = gϕ(y)/g ′ ϕ (y) gϕ(ȳ)/g ′ ϕ (ȳ) (the argumentation is very closed to [18,Sec 3.4] then we refer the reader to this article for the details). Finally, the boundary condition (39) can be rewritten as (37).
With the same reasoning as [18,Sec. 3.3], we can prove that D ϕ (y) is analytic in G Lϕ , has finite limits on L ϕ \ {1} and then is continuous on G Lϕ \ {1}. We start with proving the analytic continuation of D ϕ (y) from the unit disc D to G Lϕ . From the functional equation (25), for y ∈ {y ∈ C : |X ϕ,0 (y)| 1} ∩ D}, The latter equation defined on G Lϕ ∩ D can be continued in G Lϕ ∪ D. The generating function D ϕ (y) is analytic on D, and on G Lϕ ∪ D \ D , D ϕ (y) has the same singularities as δ ϕ (y) and X ϕ,0 (y), namely the branch cut [y 1 , y 4 ]. However, this segment does not belong to G Lϕ ∪ D\D , then D ϕ (y) can be analytically continued to G Lϕ . Then we prove that D ϕ (y) has finite limits on L ϕ \ {1}. From equation (38), we only need to study the zeros of . Thanks to the relation X ϕ,0 (Y ϕ,0 (x)) = x valid in G Mϕ (see [11,Cor. 5.3.5]), it recurs to study the zeros of δ ϕ for y in G Lϕ ∪ D \ D. The discriminant δ ϕ vanishes at y ϕ,1 , 1, y ϕ,4 and only 1 belongs to the last set. Then D ϕ (y) has finite limits on L ϕ \ {1}.

Solution of the boundary value problem
In this section we solve the boundary value problem stated in Lemma 12 and obtain in Theorem 14 an explicit expression for the generating function D ϕ (y) = i 1 f (i, i)y i−1 . Finally, Theorem 15 gives a formal expression of the generating function H(x, y) = (i,j)∈C f (i, j)x i−1 y j−1 . Along the proof Theorem 14, we also state Lemma 17 which gives an expression of D ϕ (y) for a family of harmonic function non necessarily positive.
Theorem 14. The generating function D ϕ (y) can be written as with θ ϕ , w ϕ (y) and W ϕ defined in Subsection 3.1 and δ ϕ in (27).
Theorem 15. Let f be a harmonic function associated to a random walk in the three-quarter plane with hypotheses (H1), (H2), (H3) and (H4). The generating function H(x, y) of f can be written as Remark 16. The expression of H(x, y) in Theorem 15 is formal, as the series in (2) is not convergent. We will rather consider H as the sum of generating functions, as in (19).
There exists a polynomial P ∈ R[y] such that Remark 18. More generally, for any polynomial P of degree n, we get a harmonic functions satisfying (P1) and (P2) (but not necessarily (P3)). Remark 20 illustrates this statement.
Proof of Lemma 17. The conformal gluing function w ϕ (y) defined in Subsection 3.1 transforms the domain G Lϕ into the complex plane cut by the segment I ϕ = [ w ϕ (Y ϕ (x ϕ,1 )), 1). Let w −1 ϕ be the inverse function of w ϕ . We denote by w −1 ϕ (z) + and w −1 ϕ (z) − the left and right limits of w −1 ϕ on I ϕ . The latter are complex conjugate on I. Let and The set Z ϕ is finite, otherwise by the principle of isolated zeros, gϕ(y) g ′ ϕ (y) D ϕ (y) would be equal to zero on the whole of G Lϕ , and let N ϕ be the cardinal of Z ϕ . Together with Lemma 12, we have is bounded in C, then by Liouville's theorem, the latter is constant in the whole complex plane. Then there exists a polynomial P of degree N ϕ such that Proof of Theorem 14. Take P of degree one in (43), the generating function D ϕ (y) can be written as The expression of gϕ(y) g ′ ϕ (y) where the function g is defined in (40) can be simplified into . Furthermore, as the function w is decreasing on (X(y 1 ), 1), taking the square root of (35) gives Then, thanks to equation (30), there exists a c g = 0 such that for y in the neighborhood of 1, and finally, there exists a c ′ D = 0 such that The asymptotics of D ϕ (y) in the last equation (45) is exactly the asymptotics found in [15,Eq. (1.4)], see Remark 19. By unicity of the positive harmonic function, P is a polynomial of degree one. Going back to (44) and noticing that g ϕ (0)/g ′ ϕ (0) = 0, we deduce ν = − w ϕ (0) = 0. Then, computing the first term in the expansion of (44) at y = 0, we find µ = f (1,1) w ′ ϕ (0) . Remark 19. In [15], Mustapha computes the asymptotic of the number of small steps walks excursions in the three-quarter plane. In particular, Equation (1.4) shows that the polynomial exponent α C in the three quadrants can be easily expressed as a function of the critical exponent α Q in the quadrant. From this equation we have for y close to 1, a constant c H = 0 such that H(x, 0) = c H +o (1) (1−x) π/(2π−θ) , and we deduce that there exists c D = 0 such that D ϕ (y) = c D +o (1) (1−y) π/(2π−θ)+1 . Let us point out that the positivity of the harmonic function (property (P3)) is crucial in the proof of asymptotic results in [15].
Proof of Theorem 15. The expression of D ϕ (y) in Theorem 14 suffices to calculate H(x, y). Indeed, to have an expression of H(x, y), we only need to have one of D(x, y) and L(x, y) (thanks to the symmetry of the problem we have U (x, y) = L(y, x)). We begin by plugging in Y ϕ,0 (x) to the functional equation (25) and get Together with (25), we can write an expression of L ϕ (x, y) and thanks to the symmetry, we get for free an expression of U ϕ (x, y) = L ϕ (x −1 , y), and with a series expansion of L ϕ (x, y), U ϕ (x, y) and D ϕ (y), the value of f (i, j) for all (i, j) ∈ C. Nevertheless, if we still want to write a formal expression of the generating function H(x, y) defined in (2), we can use the following relations and from them we can write We have now all the ingredients to write the following theorem which gives an expression of H(x, y).

Example of the simple random walk
In this subsection, we apply the results of Subsection 3.3 to the simple random walk. The change of variable ϕ defined in (24) transforms the simple random walk into the Gessel random walk (see Figure 7). We first need to compute the angle θ ϕ , the conformal gluing functions W ϕ , W ϕ and w ϕ . Then we get an expression of D ϕ (y) and L ϕ (x, 0) and end up with a series expansion of H(x, y). The simple walk is defined by p 1,0 = p −1,0 = p 0,−1 = p 0,1 = 1 4 . We easily have θ = π 2 and together with (36), we deduce θ ϕ = 3π 4 . From Equations (5) and (26), the kernel after the change of variable is defined by Conformal gluing functions. On one side, we have On the other side, Series expansion of D ϕ (y) and L ϕ (x, 0). Noticing that w ϕ (0) = −16/9, W ϕ (0) = 0, thanks to Thanks to (46) with K ϕ (x, 0) = 1/4, we have Series expansion of L ϕ (x, y). To begin with, we write an expression of K ϕ (x, y)L ϕ (x, y). With (25), we get Remark 20. The harmonic function f (i, j) = ij is a positive harmonic function of the simple walk in the quarter plane. In the three-quarter plane, this function is still harmonic, but does not satisfy (P3). However, there exists a polynomial P as in (43) such that D ϕ (y) = i 1 f (i, i)y i−1 = 1+y (1−y) 3 . A quick study of exponent in the same idea of Section 3.3 shows that P should be of degree 2. Then there exists (a, b, c) ∈ R such that We have , and evaluating this last equation with y = 0 gives c = 0. Then, dividing by w ϕ (y) and letting y go to 0 gives b = − 9 16 . Finally, examining the first and second term of the expansion of (56) at y = 0 gives a = 3 4 . Vice versa, letting P (y) = 3 4 y 2 − 9 16 y, when we develop D ϕ (y) in series, we find D ϕ (y) = 1 + 4 y + 9 y 2 + 16 y 3 + 25 y 4 + O y 5 . (57) L(x, y) U (x, y) Mixing equations (21), (59) and (22) and multiplying by xy, we get (58).
In this system of two functional equations, the bivariate generating function U (x, y) (resp. L(x, y)) is related to the bivariate generating functions D(x, y), D ℓ (x, y) (resp. D u (x, y)) and the univariate generating function U −0 (x −1 ) (resp. L 0− (y −1 )). In order to simplify this system of functional equations, we perform two changes of variables ϕ L for the lower part and ϕ U for the upper part with ϕ U (x, y) = (x, x −1 y) and ϕ L (x, y) = (xy, x −1 ).
This change of variables transforms the upper part {j 1, i j − 1} into the left quadrant {i −1, j 0} and the lower part {i 1, j i − 1} into the right quadrant {i 1, j 0} (see Figure 8). Note that the diagonal is changed by both ϕ L and ϕ U into the positive y-axis (see Figure 8). The step set is changed as well, and walks in the three-quadrant can be seen as inhomogeneous walks in the half plane with two different step sets in each quadrant and a mixed step set on the positive y-axis (see Figure 9). Other type of inhomogeneous walks are studied in [5,8].
In the symmetric case (symmetry of the random walks and hence the symmetry of the unique positive harmonic function), the study of harmonic functions of random walks in the three quadrant is then equivalent to the study of harmonic functions of random walks reflected on the y-axis and constrained by the x-axis in the positive quadrant (see Figure 10). Articles [2,19] work on problem in the same vein: there authors study walks in the quadrant with different weights on the boundary.
Unfortunately, due to number of unknown functions, we are not able to solve this system of functional equations yet. Let us end this section by pointing out that the split cone along the diagonal can also be related to the Join-the-Shortest-Queue model (JSQ). Supposed that there are two lines (see Figure 11, left), each of them with a service time exponentially distributed of rate r 1 where L(x, y) = In order to simplify the last functional equation (61), we apply the following change of variables ψ(x, y) = (x, x −1 y).
The equation (61) is changed into