On The Riemann Function

Riemann’s method is one of the definitive ways of solving Cauchy’s problem for a second order linear hyperbolic partial differential equation in two variables. The first review of Riemann’s method was published by E. T. Copson in 1958. This study extends that work. Firstly, three solution methods were overlooked in Copson’s original paper. Secondly, several new approaches for finding Riemann functions have been developed since 1958. Those techniques are included here and placed in the context of Copson’s original study. There are also numerous equivalences between Riemann functions that have not previously been identified in the literature. Those links are clarified here by showing that many known Riemann functions are often equivalent due to the governing equation admitting a symmetry algebra isomorphic to SL(2, R). Alternatively, the equation admits a LieBäcklund symmetry algebra. Combining the results from several methods, a new class of Riemann functions is then derived which admits no symmetries whatsoever.


Introduction
The interest in Riemann's method is long-standing.The reason is that once the Riemann function is determined, the governing equation can be solved for Cauchy data on any other non-characteristic curve.The value of such a property means that Riemann's method continues to draw the attention of investigators today.Some applications include solving electromagnetic problems exhibiting rotational symmetry [1], finding existence criteria for the eigenvalues of the solution of focal point problems [2], solving for the solution of transient plane waves [3] and the inverse problem of scattering theory [4]- [12].More recently, Riemann's method has been applied to the solution of coupled Korteweg-de Vries equations [13], to boundary value problems for the non-homogeneous wave equation [14]- [18], to the solution of the non-linear Schrödinger equation [19]- [20] and modelling hyperbolic quasi-linear equations [21]- [23].
Copson [24] wrote the first review of Riemann's method in 1958.In total, he listed six different techniques to solve for the Riemann function.However, three other approaches were missed by Copson.They are included here for completeness.Since 1958, another four methodologies have emerged for finding Rieman functions.Often the resulting Riemann functions are not new but rather a reduction of some more general Riemann function or obtainable from another Riemann function by a change of variables.This fact is frequently overlooked in the literature.Here the new methods are outlined to extend Copson's review and the equivalences are clarified.One finding is that the Riemann functions are equivalent because the governing equation admits a symmetry algebra isomorphic to SL (2, R).Alternatively, the equation admits a Lie-Bäcklund symmetry algebra.By combining several of the solution techniques, a new class of Riemann functions is then obtained that admits no symmetries whatsoever.Consider the partial differential equation (PDE) in characteristic variables The aim is to represent a solution U at a point P 0 in terms of the initial data, which are the values of U and one outgoing derivative of U on the initial curve C. Thus both U r = c 1 and U s = c 2 , where c 1 and c 2 are constants, are known on C from such data as shown in Figure 1.Take the adjoint and define V = R(r, s, r 0 , s 0 ) such that ∂ ∂r R(r, s, r 0 , s 0 ) = b(r, s 0 )R(r, s, r 0 , s 0 ), on s = s 0 , ∂ ∂s R(r, s, r 0 , s 0 ) = a(r 0 , s)R(r, s, r 0 , s 0 ), on r = r 0 , R(r 0 , s 0 , r 0 , s 0 ) = 1.( 6) R(r, s, r 0 , s 0 ) is called the Riemann function.It is the solution of the characteristic boundary value problem for the adjoint equation denoted by (3)- (6).R(r, s, r 0 , s 0 ) will hold for arbitrary initial values given along an arbitrary noncharacteristic curve C. The Cauchy problem then has a unique solution given by An equivalent formulation to that just given can also be obtained by using the canonical variables With this change of variables, (1) becomes where a ′ = 2(a + b), b ′ = 2(b − a), c ′ = 4c.The Riemann function, R(x, y, x 0 , y 0 ), must now satisfy R(x 0 , y 0 , x 0 , y 0 ) = 1.(13) In the literature, both the characteristic (1) and the canonical form (9) are used.

Copson's Review
Copson [24] published the first review of Riemann's method in 1958.In all, Copson listed six different methods for finding the Riemann function.They were: 1. Riemann's original method [25], which was based on the fact that the Riemann function does not depend in any way on the curve carrying the Cauchy data.Solving the Cauchy problem by some other means for a special curve (eg. a straight line) then yields the Riemann function by a comparison of the two solutions.Riemann only gave explicit formulae for the Riemann function in the two cases which interested him from gas dynamics.The most famous of these is the Euler-Poisson-Darboux equation (EPD) which has Riemann function where z 0 = (y − y 0 ) 2 − (x − x 0 ) 2 4xx 0 (16) and 2 F 1 is the hypergeometric function (see the appendix for details of all the hypergeomtric functions used in this article).Equation ( 15) was derived via a Fourier cosine transform.2. Hadamard [26] showed that the coefficient of the logarithmic term in his elementary solution is the Riemann function for the adjoint equation.3.For separable equations, Copson found that it is straightforward to construct an integral equation whose unique solution is the Riemann function.4. Chaundy [27] - [32] was able to construct the Riemann function for a number of equations by the use of symbolic operators and power series.Of particular note was for which the Riemann function is given by where F B is a Lauricella hypergeometric function of four variables [33] and 5. Mackie [24] constructed complex integral solutions of certain equations.An appropriate choice of contour results in the Riemann function.6. Titchmarsh [36] gave a direct solution for the Riemann function of the equation of damped waves by means of a complex Fourier integral.
Equations ( 14) and ( 17) have been singled out here as they occur repeatedly in the literature.Equation ( 14) is a subcase of (17) in canonical variables, where m 2 = m 3 = m 4 = 0.In fact in 1958, (17) was the most general self-adjoint equation for which the Riemann function was known.

Self Adjoint Riemann Functions
In 1979, a summary of all the known self-adjoint Riemann functions was published by Lanckau [37].For an arbitrary equation, the cases he listed were: 1. C(r, s) = C 0 , which has Riemann function where C o is an arbitrary constant, J 0 is the modified Bessel function of zero order and and P −m is the Legendre function of order −m, with z 1 given by (19).This characteristic form is also equivalent to the canonical representation ( 14) and ( 15) of the EPD equation.

C(r, s)
where Ξ 2 is a hypergeometric function of two variables [34] and z 5 and z 1 are given by ( 23) and ( 19) respectively. 4.
and F 3 is a hypergeometric function of two variables [34], with z 1 , z 2 given by (19).The representation of its Riemann function in terms of F 3 was first obtained by Henrici [35].

C(r, s)
where F B is a Lauricella hypergeometric function [34] of three variables z 5 , z 1 and z 2 , which are given by ( 23) and (19).
Each of the above cases is easily obtained from Chaundy's equation (17).Consider case four.If (17), the corresponding terms in the multiple power series (18) are replaced by unity.Consequently, (17) reduces automatically to the equation given in case four and ( 18) becomes (26).Likewise, for case two.Cases one and three require a confluence.Starting with an equation in the form of case two, make the change of variables Taking the limit as ǫ → ∞ yields the equation from case one, with the corresponding Riemann function (22).A similar argument shows that case three is obtainable from case four.For case five, start with equation ( 17) and set m 4 = 0. Make the change of variables and let ǫ → 0 to obtain the equation from case five.The one equation listed by Lanckau, which has not previously been linked to Chaundy's equation and its contractions, is where q k , a k , b k , c k , d k are constants subject to the conditions that and This equation was originally published by Püngel [38,39].Neither Püngel, nor Lanckau, made any connection from this case to Chaundy's equation (17).However, standard calculations show that Püngel's equation is Möbius equivalent to Chaundy via Equation ( 28) acts as a normalization constraint while (29) restricts the parameters to Chaundy.Püngel's case is indicative of the historic difficulty of demonstrating equivalences between Riemann functions.

Method's Not Included by Copson
The six approaches listed by Copson [24] were all that were known to him.Likewise the equations listed in section (2.1) were the only cases he was aware of, apart from trivial changes of the dependent or independent variables.However, another three methods were not included by Copson.In 1937, Courant and Hilbert [40] derived the Riemann function for the Telegrapher's equation via Lie-Point symmetries; although this was only really made clear in a 1962 translation of their benchmark text "Methoden Der Mathematischen Physik II".Their approach represents the first use of symmetry groups for Riemann's method.Cohn [42] developed an iterative technique to obtain the Riemann function.Also Olevskiȋ [45], derived an addition theorem for the Riemann function of an equation based on the Riemann functions of two simpler separable equations.

The Telegrapher's Equation
In the original German text, Courant and Hilbert referred to their method as "der Symmetrie der Differentialgleichung den ansatz".But what this in fact means is a little obscure.In 1937, the process of similarity reduction, instigated by Lie [41], was still embryonic.To understand what was intended, the expanded 1962 translation is more instructive.Take the Telegrapher's equation with constant C 0 .Courant and Hilbert's argument is then the following.Since the operator L has constant coefficients, R(r, s, r 0 , s 0 ) depends only on the relative position of the points (r, s) and (r 0 , s 0 ).Moreover, letting (r 0 , s 0 ) be the origin, we observe that if satisfies conditions (3) to (6) imposed on the Riemann function, then means (since a = b = 0) that V = R remains constant along the coordinate axes.If V has this property, then W has it also; and finally V(0, 0) = 1 implies W(0, 0) = 1.Since these conditions determine Riemann's function uniquely, it follows that W(r, s) = V(r, s) is a function of rs.For general (r 0 , s 0 ), the Riemann function then has the form The equation L * [R] = 0 then yields the equation z f Consequently the Riemann function is Clearly what Courant and Hilbert had in mind was that the PDE for the Riemann function and the associated boundary conditions are invariant under a three-parameter Lie transformation group.The details of their argument reflect the emerging nature of symmetry groups as a method of solution in this area.More systematic use of symmetry groups to find Riemann functions have subsequently been developed by Bluman [48], Daggit [50], Ibragimov [51] and Iwasaki [67], which will all be explored in section 4.

Successive Iterations and the Banach Fixed Point Principle
In 1946, Cohn [42] considered equations of the form where H is an arbitrary function of its argument to be determined.Cohn sought to prove the existence of R(r, s, r 0 , s 0 ) by iteration.Hence, he rewrote (31) and condition (6) as the single equation where B[r, s] is the rectangle formed by the horizontal and vertical lines through (r, s) and (r 0 , s 0 ).By continual recursive substitution of the right hand side of (32) into itself, he obtained the formal iterated series The series (33) for R begins with the terms 1 − ∆ + . . .where This iterative process is also known as the Banach fixed point principle [43].Now ∆ equals zero, whenever r = r 0 or s = s 0 , thus R(r, s, r 0 , s 0 ) is also a constant for these values of r or s.This led Cohn to try a Riemann function of the form Substituting ( 34) into ( 31) and ( 6), one finds that Expanding H(r + s) = H(ω) as a power series in ω, the coefficients of the series must vanish, yielding the non-linear ordinary differential equation for and λ, µ, ν are real valued constants.After a few further calculations, Cohn found that the Riemann function for (31), where H is given by (35), is and The form of H immediately suggests a few simpler limiting cases.For example, take λ → λe µν and rename λ → r+s) .Applying the limiting process of confluence, the Riemann function is then R(r, s, r 0 , s 0 ) = J 0 2 µ λ(e µr 0 − e µr )(e µs − e µs 0 ) .
Similarly, letting µ → 0, the Telegrapher's equation (30), is recovered.Apparently unaware of the work of Cohn, in 1986, Vaz et al. [44] proved an existence and uniqueness theorem for the iterative process (33).They considered the more general problem where Reformulate equation (39) in terms of the linear Volterra integral equation as defined by Cohn.Writing the solution of (40) as the existence and uniqueness of (41) was then proven.Cohn chose a specific functional form for f (r, s).In principle, the approach can be applied to more general functions -hence the usefulness of the theorem.Vaz et al. subsequently used (41) to derive the Riemann function for the Telegrapher's equation (30).

Olevskiȋ's Addition Formula
Copson's method three from his review paper [24] proposed a way to derive Riemann functions that extends Riemann's own contribution to the field.It relies on the separability of the governing equation.However, it turns out that an addition formula, which encompasses Copson's approach, was already in the Russian literature in 1952 [45].This was not noted until 1977, in an erratum [47] to Papadakis and Wood's paper [46], which rederived Olevskiȋ's result.
In Olevskiȋ's notation, the Riemann function, R ρ 1 −ρ 2 , of the equation is Alternatively, integrating by parts yields the equivalent formula Olevskiȋ obtained this result in part by using the method of successive approximations, which mirrors Cohn's approach in Section (3.2), but there are not many details in the paper.Mostly, Olevskiȋ limits himself to showing that ( 43) satisifes ( 42) and the conditions ( 11)-( 13).
To apply this addition formula efficiently, the following identity, which was first employed by Chaundy [30], can be useful.Within the region of uniform convergence of the infinite series that define the indicated functions, we have Olevskiȋ included three specific examples in his paper, illustrating the use of the addition formula.We shall consider one of those examples here.Firstly, Olevskiȋ stated that the Riemann function for where No derivation for this Riemann function was given.In particular, it is not on the list of self-adjoint equations ( 21) from section (2.1).Nevertheless, interchanging the roles of x and y and letting m → n yields with the associated Riemann function, and Hence applying (44), the Riemann function of where Now make the change of variables Using (45), this is recognised to be where z 7 and z 8 are defined by (48) and (51).Equations ( 46) and ( 52) are in fact subcases of ( 17).This will be explored further in sections (4.1) and (4.5).Again, this shows the historic difficulty in identifying equivalences in the literature.

Developments Since 1958
Several new constructive techniques have been proposed since the publication of Copson's paper.The resulting Riemann functions are often equivalent, under a change of variables, although that fact has often been missed in the literature.The aim here is to outline the solution methods and clarify the equivalences.A key finding is that the equivalence can be defined in terms of the governing equation admitting a symmetry algebra isomorphic to SL(2, R) or to a Lie-Bäcklund symmetry algebra.

Lie Point Symmetries
As stated in section (3.1), the use of symmetry groups to find Riemann functions dates to 1937.However, a modern treatment of the technique did not emerge until 1967, in the Ph.D. thesis of Bluman [48,49].Subsequently, the approach has become an active area of investigation [50]- [54].Bluman's method gave a completely algorithmic way of applying Lie point symmetries to Riemann's method and thus further developed the technique pioneered by Courant and Hilbert.
Bluman's treatment employed the (now) well-understood infinitesimal representation for deriving symmetry reductions.Employing the definitions and terminology of Olver [55], let where x = (x 1 , x 2 , x 3 , . . ., x p ) ≡ (x, y, z, . ..), u = (u 1 , u 2 , . . ., u q ) ≡ (u, u x , u y , u z , . ..), be a vector field defined on an open subset, M ⊂ X × U, of the space of independent and dependent variables.The n-th prolongation of v is the vector field, defined on the corresponding jet space, M (n) ⊂ X × U (n) , with the second summation being over all unordered multi-indices J = (j 1 , . . ., j k ) and 1 The coefficient functions φ J α of Pr (n) v are given by the following formula: The total differential, D j , is given by We also have that is an n-th order system of differential equations of maximal rank defined over M ⊂ X × U.If G is a local group of transformations acting on M, and for every infinitesimal generator v of G, then G is a symmetry group of the system.Taking these definitions, Bluman analysed the EPD equation written as which is slightly different to (14) 1 .
1 Any equation U yy − U xx − f (x)U x = 0, where f (x) is an arbitrary function, may be put into self-adjoint form by the Hence, the change of variable U = x −m V(x, y) will convert (59) to ( 14).
Bluman's idea was to follow the 1962 lead of Mackie [56] and find the Green's function of (59).Mackie showed that the relationship between the Riemann and Green's functions can be written as R(x, y, x 0 , y 0 ) = −2G(x 0 , y 0 , x, y) if (x, y) lies inside P 0 AB 0 if (x, y) lies outside P 0 AB (60) where P 0 AB was defined in Figure 1.So given the vector field, the symmetry group of ( 59) is found by taking to obtain a system of seven determining equations for the infinitesimals ξ, η and α.Solving these equations produces a four-dimensional symmetry group.The presence of the delta functions (and therefore of their derivatives) gives rise to the extra condition ξ(x 0 , y 0 ) = η(x 0 , y 0 ) = 0. Applying this condition, Bluman obtained the specific group The similarity variables are defined by the integrals x , with V(z) an arbitrary function of z.Substituting (61) into (59), Bluman found that Letting ψ = (2x 0 − z)/(4x 0 ), (62) becomes the hypergeometric equation The solution displaying the desired properties, (3)-( 6), is where A is a normalization constant to be determined.Putting this all together and using (60), the Riemann function for ( 59) is then In 1970, Daggit [50] also employed symmetry groups in an effort to find Riemann functions.His idea was to obtain conditions on the coefficients a, b and c of (1) which allow a similarity reduction to the Riemann function.Daggit found three distinct cases in which the Riemann function could be found.Kokinasidi [52] also presented similar results to those given by Daggit.Extending Daggit's results, Wood [57] defined any Riemann function where M and G are functions, to be a simple Riemann function when (i) the equation is almost self-adjoint, (ii) the ODE for f has coefficients depending on G, and (iii) the characteristic conditions for the Riemann function become initial conditions for f 2 .All the results of Daggit [50], Riemann [25] in section (1) and Cohn [42] from Section (3.2) are simple Riemann functions.Wood concluded that, for self-adjoint equations, all the simple Riemann functions had been found by these three authors.Geddes and Mackie [59] arrived at essentially the same conclusion when they found Cohn's method [42] does not produce any new results even if one seeks a Riemann function R(z; r 0 , s 0 ) where R satisfies an ODE.
In other words, the equations are essentially the only possibilities.The first three are the Telegrapher's equation, the EPD equation and Cohn's equation.The fourth equation is (46) from section (3.3) that was cited, without derivation, by Olevskiȋ.The last one appears to be new.Wood's statement that these are the only possible cases needs clarifying.In [60], Bluman showed how to construct invertible point transformations, which map a given PDE into another PDE, in the sense that any solution of the given PDE is mapped into a solution of the target PDE.In general such a mapping need not be a group transformation.Although, if the mapping from a given PDE to a target PDE is one-to-one (invertible) then the transformation must establish a one-to-one correspondence between the infinitesimal generators of the given PDE and the target PDE.More precisely, it is necessary that any given Lie algebra of infinitesimal generators of the given PDE be isomorphic to a 2 Wood's assumed form of the solution, when taken in light of the fact that he is essentially discussing similarity solutions, can be thought of as a first step towards the direct method of symmetry calculation derived by Clarkson and Kruskal [58] in 1989.
Lie algebra of infinitesimal generators of the target PDE.Recall that a Lie algebra L x with basis set {X 1 , X 2 , . . ., X n } will be isomorphic to the Lie algebra L z with basis set {Z 1 , Z 2 , . . ., and where is the usual commutator and the structure constants {C γ αβ } are the same for L x and L z .A few calculations show that the symmetry algebras of ( 65) -( 68) satisfy ( 69) -( 71) and ( 64) is obtained via a confluence.In effect, ( 64) -( 68) form an equivalence class.
In 1991, Ibragimov [51], unaware of the research of Daggit et al., published a unifying result connecting all previous symmetry group studies for Riemann's method.Again, Ibragimov looked at predicting the form of the Riemann function based on the coefficients a, b and c of (1).The idea is that any linear hyperbolic PDE can be classified based on the size of the group that it admits.Ovsiannikov [61] showed that (1) admits a four-dimensional Lie algebra, if and only if, the invariants (Ovsiannikov invariants) are constants (if h = 0 then swap k and h).When p or q (or both of them) is not constant then the symmetry algebra is two dimensional.Here h and k are Laplace invariants of ( 1) When p and q are constant and q = 0, one obtains case one as found by Daggit.Whereas Daggit's cases two and three correspond to q = 0.As calculated by Bluman, the symmetry algebra of the EPD equation is four dimensional.Furthermore, the symmetry algebra of the EPD equation is isomorphic to SL(2, R) [51].Using standard results from Lie theory [62], the symmetries exponentiate to a local Lie representation of the group SL(2, R) by operators T(G), where Using this fact, one can show that the fundamental solution can be represented in terms of the hypergeometric function.In effect, all the Riemann functions given in this section are encompassed by this result.
In 2015, Andrey et al. [54] studied the equation for which the Riemann function was obtained via a symmetry reduction as Applying (72) shows that p = 1 and q = 0. Hence, (73) falls into case two of Daggit and it admits a four dimensional symmetry algebra.Equation ( 73) is simply a confluent form of the EPD equation ( 14), under a transformation of variables.
In 2003, an interesting qualitative property linked to (72) was published in [53].There it was shown that if k = 0 and q is a non-zero constant in (72), then (1) can be written as which can be factorized to give The Riemann function is then easily obtained from the factorization.

Laplace Transform for a Klein-Gordon Equation with a Non-Constant Coefficient
The first paper exploring the application of Laplace transforms to Riemann's method was due to Scott [63].In 1977, Wahlberg [64] also used the technique to find the Riemann function for To solve (77), first make the change of variables, Equation ( 77) then reduces to where σ = λ −2/3 (1 + λx 0 ).Note that from (72), p = 1 but q = const.Hence, ( 79) is not part of the equivalence class of equations admitting a symmetry group isomorphic to SL(2, R) that was seen in section (4.1).By extension, there is no invertible transformation from (79) to the EPD equation (14).
Equation ( 79) is then solved by Laplace transform.Write the transform of W(r, s) as The transformed version of ( 79) is after applying the boundary condition (80).The solution to (81) is easily obtained as The function F is determined from the condition that W(0, s) = 1, or W(0, p) = 1/p, which yields Let p = ξ + iη, the solution to ( 79) and ( 80) is then formally given by where W(r, p) was defined in (83).Expressing (84) as a contour integral, Wahlberg showed that the Riemann function for (77) can then be represented by where Γ is defined by |z| = ρ, encircling the origin in the positive direction (since the integrand is analytic everywhere except at z = 0, Γ can be taken to be any closed contour encircling the origin in the positive direction).Alternatively, expanding exp(−δ 2 z 3 /12) of the integrand as a power series then applying the formula [65] for the Bessel function of order n, yields with

The Multiplication Formula
In 1981, Xin Hua Du [66] derived a multiplication theorem for Riemann functions.Returning to the method of successive approximations, first used by Cohn [42] and Olevskiȋ [45]  as where G 1 = 1 and For simplicity write From this starting point, Du supposed that A(x, y) = ρ(x) in (89).In this way the corresponding Riemann function (90) can be denoted by R ρ (x, y, x 0 , y 0 ).Then the Riemann function, R ln ρ , for the equation Conversely, if R ln ρ is given, then similarly we have can be given by where ( ) • ( ) represents a convolution.Furthermore R ρ 1 and R ρ 2 must possess the properties The above notation for (91)-( 95) is exactly that used by Du.However it is somewhat cryptic.So a little explanation may be useful.What Du denotes by [ ] * in (95) is essentially the Riemann function for This is built using Olevskiȋ's addition formula (43) and the equations for which the required Riemann functions can be obtained from (92).Then applying (93) to (98) yields the required Riemann function for (94).

Finite Groups and a Link to Appell's F 4
First, return to the equation Recall that it was seen as case 4 of ( 21), where it was obtained as a reduction of Chaundy's equation (17).The Riemann function in characteristic variables was where z 1 and z 2 are given by (19).Applying (72), it is straight forward to show that the Riemann function for (99) cannot be obtained by means of point symmetries as the group is trivial.However, Iwasaki [67] found in 1988 that, although a symmetry reduction will not yield the Riemann function, a finite group isomorphic to Z 2 × Z 2 × Z 2 acting on its Riemann function permitted him to reduce (99) to Appell's system F 4 .As Iwasaki showed, the Riemann function R(r, s, r 0 , s 0 ) is invariant under the transformations of the independent variables of the form 1. (r, s, r 0 , s 0 ) → (λr, λs, λr 0 , λs 0 ), λ constant, 2. (r, s, r 0 , s 0 ) → (1/r, 1/s, 1/r 0 , 1/s 0 ), 3. (r, s, r 0 , s 0 ) → (s, r, s 0 , r 0 ), 4. (r, s, r 0 , s 0 ) → (r 0 , s 0 , r, s).
Taking advantage of (1), define the new variables The Riemann function must now satisfy The group G generated by the transformations of (X, Y, Z) have the generators As is easily seen, X * , Y * and Z * are involutions and commutative with each other.Hence the group Iwasaki then sought an extension to the field K/k where K = C(X, Y, Z) and k = { f ∈ K; f is G-invariant}.With this in mind, he defined the new variables where (i, j, k) runs over all permutations of (1, 2, 3) and X 1 = X, X 2 = Y and X 3 = Z.Thus p i , q i ∈ k.
After further examination of the field, Iwasaki was led to try Substitution of (106) into ( 101) and (102) gives where and s 3 is not an element of k, (107) splits into two parts: Iwasaki then proposed that if f (r ′ , s ′ , t ′ ) is a solution of (108) (suppose that f makes sense at t ′ = 0), then f (r ′ , s ′ , 0) is also a solution of (108).Hence it is reasonable to expect that R(r, s, r 0 , s 0 ) is a function depending only on (r ′ , s ′ ) and to consider a system of PDEs to which (108) reduces if this expectation is correct.Now from (102), s = 1 implies u = 1.If (109) has a solution satisfying this condition then the assumption is correct and R will be given by such a solution.Unfortunately such an expectation is probably specific to this example and may prevent the possibility of building a general method based on Iwasaki's work.In fact, it is reminiscent of the type of inspired guess for which many classical Riemann function papers have been criticised.In any event, the substitutions ) into a system of partial differential equations associated with Appell's hypergeometric function F 4 (α, β, γ, γ ′ ; p, q), where α etc., depend on m and n.Iwasaki then shows that the general solution can be expressed as a linear combination of four contiguous F 4 functions.After further calculations, he concludes that the Riemann function for (99) is where and

New Riemann Functions
Based on the link between (124) and ( 52), seen at the end of section (4.5), the possibility now exists to construct new Riemann functions, which do not fit into either the point symmetry equivalence class, nor that defined by the generalized symmetries of (17).Taking separable equations and applying them to (43) will result in equations for which there is no invertible transformation to a known Riemann function as the group of the resulting equation will be trivial.Among the separable equations from the previous sections, four possibilities for application to the addition formula arise straight away.If we let a = λ 1 and m 1 = m 3 = 0 in (124), a = λ 2 and m 2 = m 3 = 0 in (129), as well as (x, y) → (η, ξ) in ( 14) and (77), we obtain The Riemann functions for the first three equations are respectively where z 16 and z 22 were given by ( 126) and (127).Similarly, z 0 , z 19 and z 22 were defined by ( 16), ( 131) and ( 132) respectively.Note also that the roles of ξ and η are interchanged, as required.
Start by combining ( 140) and (141), which produces the equation The Riemann function for (147) is then where The constants λ 1 and λ 2 were introduced trivially into (140) and (141).However, after applying these two equations to the addition formula, the ratio λ 1 /λ 2 becomes essential.In other words, it is possible to transform away either λ 1 or λ 2 , but it is no longer possible to eliminate both λ 1 and λ 2 from (147) via a change of variables.Effectively, a five parameter Riemann function has been obtained.It is useful to write the equation as (147) though, which at first glance incorporates six parameters, as the equation is then symmetric.This is now the most general self-adjoint equation for which the Riemann function is known.It incorporates one more essential parameter than Chaundy's equation (17).It was first published in [68].
In principle, all the equations (140) to (143) are applicable to the multiplication formula (94).However, no closed form Riemann function has been found using (94).This is left as an open problem for the interested reader.

Conclusion
Extending the review undertaken by Copson in 1958, seven additional approaches to the original six documented methods, have been detailed.Listing them chronologically, they are: • Lie point symmetries.
• The method of successive approximations.Furthermore, the equivalences generated by point symmetries, where the governing equation admits a symmetry algebra isomorphic to SL(2, R) have been clarified.This has now been complemented by a new equivalence class of Riemann functions which is obtainable only from a generalized symmetry and admits no non-trivial point symmetries.Hence, we have obtained a useful new diagnostic tool in the hunt for equivalences amongst PDEs and their Riemann functions.By combining several of the solution techniques, several new Riemann functions were then derived that admitted no symmetries whatsoever.
To conclude, other papers have also been published since Copson's review.For instance, the idea of the complex Riemann function has been considered [38], [72]- [74].It is largely based on the work of Vekua [75].Qualitative properties of the Riemann function have also been explored [76]- [78].As far as other qualitative properties go, only the "duality" property is known [40].Higher dimensional Riemann functions have also been explored by Zhegalov [79] and Koshcheeva [80].where (a 1 ) m 1 etc. are given by (A1).This function satisfies the system of equations ∂z j ∂z k + c − (a j + b j + 1)z j ∂U ∂z j − a j b j U = 0, for i, j ≤ n. • J n (z) This function is a particular solution of the equation • F 4 (a, b, c, c ′ , z 1 , z 2 ) where (a) m+n etc. are given by (A1).This function is the solution of the system of equations

Figure 1 .
Figure 1.The Initial Curve