Next Article in Journal
Statistical Analysis of Maximally Similar Sets in Ecological Research
Previous Article in Journal
Scaling Laws in the Fine-Scale Structure of Range Margins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

On the Riemann Function

School of Mathematics and Statistics, University of Sydney, Sydney, NSW 2006, Australia
Mathematics 2018, 6(12), 316; https://doi.org/10.3390/math6120316
Submission received: 26 October 2018 / Revised: 6 December 2018 / Accepted: 6 December 2018 / Published: 10 December 2018

Abstract

:
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 S L ( 2 , R ) . Alternatively, the equation admits a Lie-Bäcklund symmetry algebra. Combining the results from several methods, a new class of Riemann functions is then derived which admits no symmetries whatsoever.

1. 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,5,6,7,8,9,10,11,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,15,16,17,18], to the solution of the non-linear Schrödinger equation [19,20], and modelling hyperbolic quasi-linear equations [21,22,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 Riemann 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 S L ( 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
L [ U ] = U r s + a ( r , s ) U r + b ( r , s ) U s + c ( r , s ) U = 0 .
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
L * [ V ] = V r s ( a V ) r ( b V ) s + c V ,
and define V = R ( r , s , r 0 , s 0 ) such that
L * [ R ] = 0 ,
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 .
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
U ( r o , s 0 ) = 1 2 [ U R ] A + [ U R ] B + 1 2 A B R U s U R s d r + R U r U R r d s .
An equivalent formulation to that just given can also be obtained by using the canonical variables
r = y + x , s = y x .
With this change of variables, (1) becomes
U y y U x x + a ( x , y ) U y b ( x , y ) U x + c ( x , y ) U = 0 ,
where a = 2 ( a + b ) , b = 2 ( b a ) , c = 4 c . The Riemann function, R ( x , y , x 0 , y 0 ) , must now satisfy
L * [ R ] = R y y R x x ( a R ) y + ( b R ) x + c R = 0 ,
R x + R y = 1 2 ( a + b ) R on y y 0 = x x 0 ,
R x R y = 1 2 ( a b ) R on y y 0 = ( x x 0 ) ,
R ( x 0 , y 0 , x 0 , y 0 ) = 1 .
In the literature, both the characteristic (1) and the canonical form (9) are used.

2. 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:
  • 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 (e.g., 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)
    U y y U x x m ( 1 m ) x 2 U = 0 ,
    which has Riemann function
    R ( x , y , x 0 , y 0 ) = 2 F 1 ( m , 1 m ; 1 ; z 0 ) ,
    where
    z 0 = ( y y 0 ) 2 ( x x 0 ) 2 4 x x 0
    and 2 F 1 is the hypergeometric function (see the Appendix A for details of all the hypergeometric functions used in this article). Equation (15) was derived via a Fourier cosine transform.
  • Hadamard [26] showed that the coefficient of the logarithmic term in his elementary solution is the Riemann function for the adjoint equation.
  • For separable equations, Copson found that it is straightforward to construct an integral equation whose unique solution is the Riemann function.
  • Chaundy [27,28,29,30,31,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
    U r s + m 1 ( 1 m 1 ) ( r + s ) 2 m 2 ( 1 m 2 ) ( r s ) 2 + m 3 ( 1 m 3 ) ( 1 r s ) 2 m 4 ( 1 m 4 ) ( 1 + r s ) 2 U = 0 ,
    for which the Riemann function is given by
    R ( r , s , r 0 , s 0 ) = F B m 1 , m 2 , m 3 , m 4 , 1 m 1 , 1 m 2 , 1 m 3 , 1 m 4 , 1 , z 1 , z 2 , z 3 , z 4
    where F B is a Lauricella hypergeometric function of four variables [33] and
    z 1 = ( r r 0 ) ( s s 0 ) ( r + s ) ( r 0 + s 0 ) , z 2 = ( r r 0 ) ( s s 0 ) ( r s ) ( r 0 s 0 ) ,
    z 3 = ( r r 0 ) ( s s 0 ) ( 1 r s ) ( 1 r 0 s 0 ) , z 4 = ( r r 0 ) ( s s 0 ) ( 1 + r s ) ( 1 + r 0 s 0 ) .
  • Mackie [24] constructed complex integral solutions of certain equations. An appropriate choice of contour results in the Riemann function.
  • Titchmarsh [34] 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.

2.1. Self Adjoint Riemann Functions

In 1979, a summary of all the known self-adjoint Riemann functions was published by Lanckau [35]. For an arbitrary equation,
U r s + C ( r , s ) U = 0 ,
the cases he listed were:
  • C ( r , s ) = C 0 , which has Riemann function
    R ( r , s , r 0 , s 0 ) = J 0 2 z 5 ,
    where C o is an arbitrary constant, J 0 is the modified Bessel function of zero order and
    z 5 = C 0 ( r r 0 ) ( s s 0 ) .
  • C ( r , s ) = m ( 1 m ) ( r + s ) 2 , which has Riemann function
    R ( r , s , r 0 , s 0 ) = P m 1 2 z 1 ,
    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 ) = C 0 + m ( 1 m ) ( r + s ) 2 , which has Riemann function,
    R ( r , s , r 0 , s 0 ) = Ξ 2 m , 1 m , 1 , z 5 , z 1 ,
    where Ξ 2 is a hypergeometric function of two variables [36] and z 5 and z 1 are given by (23) and (19) respectively.
  • C ( r , s ) = m ( 1 m ) ( r + s ) 2 n ( 1 n ) ( r s ) 2 , which has Riemann function
    R ( r , s , r 0 , s 0 ) = F 3 m , n , 1 m , 1 n , 1 , z 1 , z 2 ,
    and F 3 is a hypergeometric function of two variables [36], with z 1 , z 2 given by (19). The representation of its Riemann function in terms of F 3 was first obtained by Henrici [37].
  • C ( r , s ) = C 0 + m ( 1 m ) ( r + s ) 2 n ( 1 n ) ( r s ) 2 , which has Riemann function
    R ( r , s , r 0 , s 0 ) = F B m , n , 1 m , 1 n , 1 , z 5 , z 1 , z 2 ,
    where F B is a Lauricella hypergeometric function [36] 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 m 3 = m 4 = 0 in (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). This is 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
r s = ϵ 2 c 0 + ( r s ) .
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
r = ϵ r s = ϵ s m 3 ( 1 m 3 ) = ϵ 2 C 0 .
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
C ( r , s ) = k = 1 4 q k ( 1 + q k ) f k 2 ( r , s ) ,
with
f k ( r , s ) = a k r s + b k r + c k s + d k ,
where q k , a k , b k , c k , d k are constants subject to the conditions that
a k d k b k c k = 1 , for all k ,
and
a l d k b l c k + a k d l b k c l = 0 , for all k and l k .
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
r ¯ = α 1 r + β 1 γ 1 r + δ 1 , s ¯ = α 2 s + β 2 γ 2 s + δ 2 .
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.

3. Methods 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 [41] developed an iterative technique to obtain the Riemann function. Olevskiĭ [42] also derived an addition theorem for the Riemann function of an equation based on the Riemann functions of two simpler separable equations.

3.1. 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 [43], was still embryonic. To understand what was intended, the expanded 1962 translation is more instructive. Take the Telegrapher’s equation
L [ U ] = U r s + C 0 U = 0 ,
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
V ( r , s ) = R ( r , s , 0 , 0 )
satisfies conditions (3) to (6) imposed on the Riemann function, then
W ( r , s ) = V ( α r , α 1 s )
also satisfies them. Clearly L * [ V ( r , s ) ] = 0 implies L * [ W ] = 0 so that (3) is satisfied. Condition (4) means (since a = b = 0 ) that V = R remains constant along the coordinate axes. If V has this property, then W also has it; 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 r s . For general ( r 0 , s 0 ) , the Riemann function then has the form
R ( r , s , r 0 , s 0 ) = f ( z ) ,
where
z = ( r r 0 ) ( s s 0 ) .
The equation L * [ R ] = 0 then yields the equation z f + f + C 0 f = 0 for f which, if we set λ = 4 C 0 z , becomes Bessel’s equation
d 2 f d λ 2 + 1 λ d f d λ + f = 0 ,
which has solution
f = J 0 ( λ ) .
Consequently the Riemann function is
V ( r , s , r 0 , s 0 ) = J 0 4 C 0 ( r r 0 ) ( s s 0 ) .
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 [44], Daggit [45], Ibragimov [46], and Iwasaki [47], which will all be explored in Section 4.

3.2. Successive Iterations and the Banach Fixed Point Principle

In 1946, Cohn [41] considered equations of the form
U r s + H ( r + s ) U = 0 ,
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
R ( r , s , r 0 , s 0 ) = 1 B [ r , s ] R ( r 1 , s 1 , r 0 , s 0 ) H ( r 1 + s 1 ) d r 1 d s 1 ,
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
R ( r , s , r 0 , s 0 ) = 1 B [ r , s ] H ( r 1 + s 1 ) d r 1 d s 1 + B [ r , s ] H ( r 1 + s 1 ) d r 1 d s 1 B [ r 1 , s 1 ] H ( r 2 + s 2 ) d r 2 d s 2
The series (33) for R begins with the terms 1 Δ + where
Δ = B [ r , s ] H ( r 1 + s 1 ) d r 1 d s 1 .
This iterative process is also known as the Banach fixed point principle [48]. 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
R ( r , s , r 0 , s 0 ) = R ( Δ , r 0 , s 0 ) .
Substituting (34) into (31) and (6), one finds that
Δ r Δ s H d 2 R d Δ 2 + d R d Δ + R = 0 , R ( 0 , r 0 , s 0 ) = 1 .
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 H
d d ω 1 H d 2 H d ω 2 3 2 H 2 d H d ω 2 = 0 ,
which has solution
H ( ω ) = λ ( λ + 1 ) μ 2 sinh 2 μ ( ω + ν ) ,
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
R ( r , s , r 0 , s 0 ) = 2 F 1 ( λ , 1 + λ ; 1 ; Φ ) ,
and
Φ = sinh μ ( r r 0 ) sinh μ ( s s 0 ) sinh μ ( r + s + ν ) sinh μ ( r 0 + s 0 + ν ) .
The form of H immediately suggests a few simpler limiting cases. For example, take λ λ e μ ν and rename λ λ / μ , μ μ / 2 , then if ν we obtain H ( r + s ) = λ e μ ( 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. [49] proved an existence and uniqueness theorem for the iterative process (33). They considered the more general problem
U r s f ( r , s ) U = 0 ,
where
U ( r , s 0 ) = σ ( r ) , U ( r 0 , s ) = τ ( s ) , U ( r 0 , s 0 ) = σ ( r 0 ) = τ ( s 0 ) .
Reformulate Equation (39) in terms of the linear Volterra integral equation
U ( r , s ) = g ( r , s ) + B [ r , s ] f ( α , β ) U ( α , β ) d α d β
for which g ( r , s ) = σ ( r ) + τ ( s ) σ ( r 0 ) and B [ r , s ] = [ r 0 , r ] × [ s 0 , s ] , as defined by Cohn. Writing the solution of (40) as
U ( r , s ) = lim k g ( r , s ) + B [ r , s ] f ( α 1 , β 1 ) g ( α 1 , β 1 ) + = B [ r , s ] B [ r k 1 , s k 1 ] f ( α 1 , β 1 ) f ( α k , β k ) g ( α k , β k ) d α k d β k d α 1 d β 1
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).

3.3. 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 [42]. This was not noted until 1977, in an erratum [50] to Papadakis and Wood’s paper [51], which rederived Olevskiĭ’s result.
In Olevskiĭ’s notation, the Riemann function, R ρ 1 ρ 2 , of the equation
U y y U x x + ρ 1 ( y ) ρ 2 ( x ) U = 0 ,
is
R ρ 1 ρ 2 ( x , y , x 0 , y 0 ) = R ρ 2 ( x , y , x 0 , y 0 ) + y y 0 x x 0 R ρ 2 ( x , t , x 0 , 0 ) t R ρ 1 ( t , y , 0 , y 0 ) d t .
R ρ 1 and R ρ 2 are the Riemann functions for
U y y U x x + ρ 1 ( y ) U = 0
and
U y y U x x ρ 2 ( x ) U = 0 .
Alternatively, integrating by parts yields the equivalent formula
R ρ 1 ρ 2 ( x , y , x 0 , y 0 ) = R ρ 1 ( x , y , x 0 , y 0 ) + x x 0 y y 0 R ρ 1 ( t , y , 0 , y 0 ) t R ρ 2 ( x , t , x 0 , 0 ) d t .
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) satisfies (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
F B ( a 1 , , a p , a 1 , , a q , b 1 , , b p , b 1 , , b q , 1 , x 1 , , x p , x 1 , , x q ) = F B ( a 1 , , a p ; b 1 , , b p ; 1 ; x 1 , , x p ) + 0 1 F B ( a 1 , , a p ; b 1 , , b p ; 1 ; ( 1 t ) x 1 , , ( 1 t ) x p ) × d d t F B ( a 1 , , a q ; b 1 , , b q ; t x 1 , , t x q ) d t .
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
U y y U x x + m ( 1 m ) sin 2 y U = 0 ,
was
R ( x , y , x 0 , y 0 ) = 2 F 1 ( m , 1 m ; 1 ; z 7 ) ,
where
z 7 = cos ( y y 0 ) cos ( x x 0 ) 2 sin y sin y 0 .
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
U y y U x x n ( 1 n ) sin 2 x U = 0 ,
with the associated Riemann function,
R ( x , y , x 0 , y 0 ) = 2 F 1 ( n , 1 n ; 1 ; z 8 ) ,
and
z 8 = cos ( x x 0 ) cos ( y y 0 ) 2 sin x sin x 0 .
Hence applying (44), the Riemann function of
U y y U x x + m ( 1 m ) sin 2 y n ( 1 n ) sin 2 x U = 0
is
R ( x , y , x 0 , y 0 ) = 2 F 1 ( m , 1 m , 1 , z 7 ) + x x 0 y y 0 2 F 1 ( m , 1 m , 1 , u 1 ( t ) ) t 2 F 1 ( n , 1 n , 1 , v 1 ( t ) ) d t ,
where
u 1 ( t ) = cos ( y y 0 ) cos ( t ) 2 sin y sin y 0 , v 1 ( t ) = cos ( x x 0 ) cos ( t ) 2 sin x sin x 0 .
Now make the change of variables
t = cos ( x x 0 ) cos ( t ) cos ( x x 0 ) cos ( y y 0 )
so that (53) becomes
R ( x , y , x 0 , y 0 ) = 2 F 1 ( m , 1 m , 1 , z 7 ) + 0 1 2 F 1 ( m , 1 m , 1 , ( 1 t ) z 7 ) t 2 F 1 ( n , 1 n , 1 , z 8 t ) d t .
Using (45), this is recognised to be
R ( x , y , x 0 , y 0 ) = F 3 ( m , 1 m , n , 1 n ; 1 ; z 7 , z 8 ) ,
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 Section 4.1 and Section 4.5. Again, this shows the historic difficulty in identifying equivalences in the literature.

4. 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 S L ( 2 , R ) or to a Lie-Bäcklund symmetry algebra.

4.1. 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 [44,52]. Subsequently, the approach has become an active area of investigation [45,46,53,54,55]. 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 [56], let
v = i = 1 p ξ i ( x , u ) x i + α = 1 q ϕ α ( x , u ) u α ,
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,
P r ( n ) v = v + α = 1 q J ϕ α J ( x , u ( n ) ) u J α ,
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 j k p , 1 k n . The coefficient functions ϕ α J of P r ( n ) v are given by the following formula:
ϕ α J ( x , u ( n ) ) = D j ϕ α i = 1 p ξ i u i α + i = 1 p ξ i u J , i α .
The total differential, D j , is given by
D j = x i + α = 1 q J u J , i α u J α .
We also have that
u i α = u α x i , u J , i α = u J α x i = k + 1 u α x i x j 1 x j k .
Now suppose that
Δ ( x , u ( n ) ) = 0
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
P r ( n ) v [ Δ ( x , u ( n ) ) ] = 0 , whenever Δ ( x , u ( n ) ) = 0 ,
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
U y y U x x 2 m x U x = 0 ,
which is slightly different to (14) (Any equation U y y U x x f ( x ) U x = 0 , where f ( x ) is an arbitrary function, may be put into self-adjoint form by the transformation U = e 1 2 x f ( t ) d t V ( x , y ) to obtain V y y V x x + 1 2 f + 1 2 f 2 V = 0 . 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 [57] 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 ) = 2 G ( x 0 , y 0 , x , y ) if ( x , y ) lies inside P 0 A B 0 if ( x , y ) lies outside P 0 A B
where P 0 A B was defined in Figure 1.
So given the vector field,
v = ξ ( x , y ) x + η ( x , y ) y + α ( x , y ) U U ,
the symmetry group of (59) is found by taking
P r ( 2 ) v [ U y y U x x 2 m x U x δ ( x x 0 ) δ ( y y 0 ) ] = 0
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
ξ = x m ( y y 0 ) , η = ( x 0 2 x 2 ) ( y y 0 ) 2 2 m , α = ( y y 0 ) .
Hence
d U ( y y 0 ) U = m d x x ( y y 0 ) = 2 m d y ( x 0 2 x 2 ) ( y y 0 ) 2 .
The similarity variables are defined by the integrals
d x d y = 2 x ( y y 0 ) ( y y 0 ) 2 + ( x 2 x 0 2 ) , d U U = m d x x ,
which give
U ( x , y ) = x m V ( z ) , z = x ( y y 0 ) 2 x + x 0 2 x ,
with V ( z ) an arbitrary function of z. Substituting (61) into (59), Bluman found that
( z 2 4 x 0 2 ) V + 2 z V + m ( 1 m ) V = 0 .
Letting ψ = ( 2 x 0 z ) / ( 4 x 0 ) , (62) becomes the hypergeometric equation
ψ ( 1 ψ ) d 2 V d ψ 2 + ( 1 2 ψ ) d V d ψ m ( 1 m ) V = 0 .
The solution displaying the desired properties, (3)–(6), is V ( ψ ) = A 2 F 1 ( m , 1 m ; 1 ; ψ ) , where A is a normalization constant to be determined. Putting this all together and using (60), the Riemann function for (59) is then
R ( x , y , x 0 , y 0 ) = x x 0 m 2 F 1 m , 1 m ; 1 ; ( y y 0 ) 2 ( x x 0 ) 2 4 x x 0 .
In 1970, Daggit [45] 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 [53] also presented similar results to those given by Daggit. Extending Daggit’s results, Wood [58] defined any Riemann function
R ( r , s , r 0 , s 0 ) = M ( r , s , r 0 , s 0 ) f ( G ( r , s , r 0 , s 0 ) , r 0 , s 0 ) ,
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 (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 [59] in 1989). All the results of Daggit [45], Riemann [25] in Section 2 method 1 and Cohn [41] 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 [60] arrived at essentially the same conclusion when they found Cohn’s method [41] 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
U r s + U = 0 ,
U r s m ( 1 m ) ( r s ) 2 U = 0 ,
U r s m ( 1 m ) sinh 2 ( r s ) U = 0 ,
U r s m ( 1 m ) sin 2 ( r s ) U = 0 ,
U r s + m ( 1 m ) cosh 2 ( r s ) U = 0 ,
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 [61], 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 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 , , Z n } if
[ Z α , Z β ] = C α β γ Z γ
and
[ X α , X β ] = C α β γ X γ ,
where
[ X α , X β ] = X α X β X β X α
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 [46], 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 [62] showed that (1) admits a four-dimensional Lie algebra, if and only if, the invariants (Ovsiannikov invariants)
p = k h , q = ln | h | r s h ,
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)
h = a r + a b c , k = b s + a b c .
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 S L ( 2 , R ) [46]. Using standard results from Lie theory [63], the symmetries exponentiate to a local Lie representation of the group S L ( 2 , R ) by operators T ( G ) , where
T ( G ) U ( t , r ) = ( α + γ t ) 2 γ 2 r 2 1 / 2 U ( δ t + β ) ( α + γ t ) γ δ r 2 ( α + γ t ) 2 γ 2 r 2 , r ( α + γ t ) 2 γ 2 r 2
and
G = α β γ δ S L ( 2 , R ) .
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. [55] studied the equation
U r s 1 2 r U s + s U = 0 ,
for which the Riemann function was obtained via a symmetry reduction as
R ( r , s , r 0 , s 0 ) = r 0 r J 0 2 ( r r 0 ) ( s 2 s 0 2 ) .
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 [54]. There, it was shown that if k = 0 and q is a non-zero constant in (72), then (1) can be written as
U r s + 1 2 U r 2 q 1 r + s U t 1 q 1 r + s U = 0 ,
which can be factorized to give
r 2 q 1 r + s s + 1 2 U = 0 .
The Riemann function is then easily obtained from the factorization.

4.2. 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 [64]. In 1977, Wahlberg [65] also used the technique to find the Riemann function for
U y y U x x + ( 1 + δ x ) U = 0 .
To solve (77), first make the change of variables,
r = 1 2 λ 1 / 3 ( x 0 + t 0 ( x + t ) , s = 1 2 λ 1 / 3 ( x t ) ( x 0 t 0 ) .
Equation (77) then reduces to
2 W r s + σ + s r W = 0 ,
W = 1 when r s = 0 ,
where σ = λ 2 / 3 ( 1 + λ x 0 ) .
Note that from (72), p = 1 but q c o n s t . Hence, (79) is not part of the equivalence class of equations admitting a symmetry group isomorphic to S L ( 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
W ˜ ( r , p ) 0 e p s W ( r , s ) d s .
The transformed version of (79) is
p W ˜ r W ˜ p + σ r W ˜ = 0
after applying the boundary condition (80). The solution to (81) is easily obtained as
W ˜ ( r , p ) = F r + 1 2 p 2 exp ( σ r ) p 1 3 p 3 .
The function F is determined from the condition that W ( 0 , s ) = 1 , or W ˜ ( 0 , p ) = 1 / p , which yields
W ˜ ( r , p ) = 1 2 r + P 2 exp ( σ r ) p 1 3 p 3 σ 2 r + P 2 + 1 3 2 r + p 2 3 / 2 .
Let p = ξ + i η , the solution to (79) and (80) is then formally given by
W ( r , s ) = 1 2 π i ξ 0 i ξ 0 + i W ˜ ( r , p ) d p ,
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
R ( ξ , η , ξ 0 , η 0 ) = 1 2 π i Γ exp 1 + δ 2 ( ξ + ξ 0 ) z + ( ξ ξ 0 ) 2 ( η η 0 ) 2 4 z δ 2 12 z 3 d z z ,
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 [66]
J n ( α ) = 1 2 π i Γ z n 1 exp 1 2 α z 1 z d z
for the Bessel function of order n, yields
R ( x , y , x 0 , y 0 ) = n = 0 Ω n n ! J 3 n ( Δ ) ,
with
Ω = δ 2 96 ( y y 0 ) 2 ( x x 0 ) 2 1 + δ 2 ( x + x 0 ) 3 / 2 ,
Δ = 1 + δ 2 ( x + x 0 ) ( y y 0 ) 2 ( x x 0 ) 2 1 / 2 .

4.3. The Multiplication Formula

In 1981, Xin Hua Du [67] derived a multiplication theorem for Riemann functions. Returning to the method of successive approximations, first used by Cohn [41] and Olevskiĭ [42] as outlined in Section 3, Du wrote the Riemann function for the equation
U y y U x x + A ( x , y ) U = 0
as
R A ( x , y , x 0 , y 0 ) = k = 1 G k ( ξ , η ; ξ 0 , η 0 ) = k = 1 1 4 k 1 Q 0 Q Q 0 Q k 1 Q 0 Q 1 i = 1 k 1 A ( Q i ) d Q 1 d Q 2 d Q k 1 ,
where G 1 = 1 and
ξ 0 ξ η 0 η A p 1 + q 1 2 , p 1 q 1 2 d p 1 d q 1 = Q 0 Q A ( Q 1 ) d Q 1 .
For simplicity write
R A = T K ( A ) .
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
U y y U x x + ln ρ ( x ) U = 0
is
R ln ρ ( x , y , x 0 , y 0 ) = [ R ρ ] ln ρ * = T k ln ( y y x x ) R ρ R ρ .
Conversely, if R ln ρ is given, then similarly we have
R ρ ( x , y , x 0 , y 0 ) = [ R ln ρ ] * ρ = T k e ( y y x x ) R ln ρ R ln ρ .
Again for convenience, write [ R ρ ] ln ρ * and [ R ln ρ ] * ρ as [ R ρ ] * and [ R ln ρ ] * . We are now in a position to state Du’s main result. If ρ 1 = ρ 1 ( x ) and ρ 2 = ρ 2 ( y ) , then the Riemann function, R ρ 1 ρ 2 , for the equation
U y y U x x + [ ρ 1 ( x ) ρ 2 ( y ) ] U = 0
can be given by
R ρ 1 ρ 2 ( x , y , x 0 , y 0 ) = R ρ 1 + ρ 2 R ρ 1 = [ R ρ 1 ] * , R ρ 2 = [ R ρ 2 ] * * ,
where ( , , , ) ( , , , ) represents a convolution. Furthermore R ρ 1 and R ρ 2 must possess the properties
R ρ 1 ( x , y , x 0 , y 0 ) = R ρ 1 ( x , y y 0 , x 0 , 0 ) = R ρ 1 ( x , y 0 y , x 0 , 0 ) ,
R ρ 2 ( x , y , x 0 , y 0 ) = R ρ 2 ( x x 0 , y , 0 , y 0 ) = R ρ 2 ( x 0 x , y , 0 , y 0 ) .
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
U y y U x x + ln ρ 1 ( x ) ρ 2 ( y ) U = 0 .
This is built using Olevskiĭ’s addition Formula (43) and the equations
U y y U x x + ln ρ 1 ( x ) U = 0 , U y y U x x + ln ρ 2 ( y ) U = 0 ,
for which the required Riemann functions can be obtained from (92). Then applying (93) to (98) yields the required Riemann function for (94).

4.4. Finite Groups and a Link to Appell’s F 4

First, return to the equation
L [ U ] = U r s + m ( m 1 ) ( r s ) 2 n ( n 1 ) ( r + s ) 2 U = 0 .
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
R ( r , s , r 0 , s 0 ) = F 3 ( n , 1 n , m , 1 m ; 1 ; z 1 , z 2 ) ,
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 [47] 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
  • ( r , s , r 0 , s 0 ) ( λ r , λ s , λ r 0 , λ s 0 ) , λ constant ,
  • ( r , s , r 0 , s 0 ) ( 1 / r , 1 / s , 1 / r 0 , 1 / s 0 ) ,
  • ( r , s , r 0 , s 0 ) ( s , r , s 0 , r 0 ) ,
  • ( r , s , r 0 , s 0 ) ( r 0 , s 0 , r , s ) .
Taking advantage of transformation 1, define the new variables
X = r s 0 , Y = s s 0 , Z = r 0 s 0 .
The Riemann function must now satisfy
L ( X , Y ) [ R ] = 0 ,
( X Z ) ( Y 1 ) = 0 implies R = 1 .
The group G generated by the transformations of ( X , Y , Z ) have the generators
X * : ( X , Y , Z ) ( X , X / Z , X / Y ) , Y * : ( X , Y , Z ) ( Y / Z , Y , Y / X ) , Z * : ( X , Y , Z ) ( Z / Y , Z / X , Z ) .
As is easily seen, X * , Y * , and Z * are involutions and commutative with each other. Hence the group G = < X * , Y * , Z * > is G Z 2 × Z 2 × Z 2 and | G | = 8 .
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
p i = X i + 1 / X i + X k / X j + X j / X k ,
q i = X i / ( X j X k ) + ( X j X k ) / X i + 2 ,
s i = X i + 1 / X i X k / X j X j / X k ,
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
r = p 1 / q 3 , s = p 2 / q 3 , t = q 3 .
Substitution of (106) into (101) and (102) gives
{ s 3 M 1 + r s t M 1 + 2 M 2 + 2 M 3 } R = 0 ,
where
M 1 = L ( r , s ) = r s + n ( 1 n ) ( r s ) 2 m ( 1 m ) ( r + s ) 2 , M 2 = ( 1 r 2 ) r 2 + ( 1 s 2 ) s 2 2 r s r s 2 2 r r 2 s s 2 n ( 1 n ) ( r s ) 2 2 m ( 1 m ) ( r + s ) 2 , M 3 = t { ( t 4 ) t ( t t r r s s 1 ) + 2 t } ,
and s 3 is defined by (105). Since M 1 R and r s t M 1 R + 2 M 2 R + 2 M 3 R are G-invariant and s 3 is not an element of k , (107) splits into two parts:
M 1 R = 0 , M 2 R + M 3 R = 0 .
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
M 1 U = M 2 U = 0 ,
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
p = r s 2 2 = ( r s ) ( r 0 s 0 ) 2 ( r s + r 0 s 0 ) 2 , q = r + s 2 2 = ( r + s ) ( r 0 + s 0 ) 2 ( r s + r 0 s 0 ) 2 , u = p n / 2 q m / 2 v ,
transform (109) 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
R ( r , s , r 0 , s 0 ) = j = 0 3 C j ( a , b ) u j ( p , q ; n , m ) ,
where
u 0 = u ( p , q ; n , m ) , 1 1 , C 0 = C ( n , m ) , u 1 = u ( p , q ; 1 n , m ) , 1 , C 1 = C ( 1 n , m ) , u 2 = u ( p , q ; n , 1 m ) , 1 , C 2 = C ( n , 1 m ) , u 3 = u ( p , q ; 1 n , 1 m ) , C 3 = C ( 1 n , 1 m ) ,
and
u ( p , q ; n , m ) = p n / 2 q m / 2 F 4 m + n 2 , m + n + 1 2 , n + 1 2 , m + 1 2 ; p , q , C ( n , m ) = 2 m n Γ 1 2 n Γ 1 2 m π Γ 1 m n .

4.5. Lie-Bäcklund Symmetries

In 2004, Zeitsch [68,69] was able to establish a new equivalence class of Riemann functions by extending the point symmetry ideas of Section 4.1 to Lie-Bäcklund symmetries. The key result was given in [68]. There it was shown that Chaundy’s Equation (17) admitted a two-dimensional generalized symmetry algebra. A straightforward calculation shows that (17) admits no non-trivial point symmetries. Following [70], define the second-order operator
S = f 1 r r + f 2 s s + f 3 r + f 4 s + f 5 .
We say that (110) is a Lie-Bäcklund symmetry operator for (17) provided
S , L = Q L ,
where
L = r s + m 1 ( 1 m 1 ) ( r + s ) 2 m 2 ( 1 m 2 ) ( r s ) 2 + m 3 ( 1 m 3 ) ( 1 r s ) 2 m 4 ( 1 m 4 ) ( 1 + r s ) 2
and
Q = h 1 r + h 2 s + h 3
is a first-order differential operator (here Q may vary with S) and f 1 , , h 3 are arbitrary functions of r and s.
Solving (111) and equating the coefficients of the derivatives to zero yields a two-dimensional vector space of operators. The symmetries, S 1 and S 2 , form a basis and have differential part
S 1 = ( r 4 + 1 ) r r + ( s 4 + 1 ) s s + 2 r 3 r + 2 s 3 s ,
S 2 = r 2 r r + s 2 s s + r r + s s .
Following [70], for each inequivalent orbit, a separable coordinate system exists. Taking a linear combination of (114) and (115), there are then four orbits, namely:
  • S 1 + 2 q S 2 , q > 1 .
  • S 1 + 2 S 2 .
  • S 1 2 S 2 .
  • S 1 + 2 q S 2 , q < 1 .
Analysing each case, Zeitsch then found:
System 1: The separable coordinates are
r = b sn a ( ξ + η ) cn a ( ξ + η ) , s = b sn a ( ξ η ) cn a ( ξ η )
where 2 q = ( 1 + b 4 ) / b 2 , k 2 = 1 b 4 , 0 < b < 1 and k is the modulus of the Jacobian elliptic functions. See [71] for more detail on the elliptic functions sn and cn as well as the other elliptic functions such as dn . When (116) is substituted into (17), the following separable equation is obtained
U ξ ξ U η η + { a 2 m 1 ( 1 m 1 ) dn 2 a ξ sn 2 a ξ cn 2 a ξ k 4 sn 2 a η cn 2 a η dn 2 a η m 2 ( 1 m 2 ) dn 2 a η sn 2 a η cn 2 a η k 4 sn 2 a ξ cn 2 a ξ dn 2 a ξ + 4 a 2 b 2 m 3 ( 1 m 3 ) dn 2 a ξ ( cn 2 a ξ b 2 sn 2 a ξ ) 2 + dn 2 ( a η ) ( cn 2 a η + b 2 sn 2 a η ) 2 1 m 4 ( 1 m 4 ) dn 2 a ξ ( cn 2 a ξ + b 2 sn 2 a ξ ) 2 + dn 2 a η ( cn 2 a η b 2 sn 2 a η ) 2 1 } U = 0 .
The Riemann function for (117) is found by substituting (116) into (18)–(20). Hence
R ( ξ , η , ξ 0 , η 0 ) = F B m 1 , m 2 , m 3 , m 4 , 1 m 1 , 1 m 2 , 1 m 3 , 1 m 4 , 1 , z 11 , z 12 , z 13 , z 14 ,
where
z 11 = ( sn a η cn a η 0 dn a ξ 0 sn a η 0 cn a η dn a ξ ) 2 ( sn a ξ cn a ξ 0 dn a η 0 sn a ξ 0 cn a ξ dn a η ) 2 4 sn a ξ sn a ξ 0 cn a ξ cn a ξ 0 dn a η dn a η 0 ,
z 12 = ( sn a ξ cn a ξ 0 dn a η 0 sn a ξ 0 cn a ξ dn a η ) 2 ( sn a η cn a η 0 dn a ξ 0 sn a η 0 cn a η dn a ξ ) 2 4 sn a η sn a η 0 cn a η cn a η 0 dn a ξ dn a ξ 0 ,
z 13 = b 2 ( sn a η cn a η 0 dn a ξ 0 sn a η 0 cn a η dn a ξ ) 2 ( sn a ξ cn a ξ 0 dn a η 0 sn a ξ 0 cn a ξ dn a η ) 2 ( cn 2 a ξ b 2 sn 2 a ξ ) ( cn 2 a η + b 2 sn 2 a η ) × ( cn 2 a ξ 0 b 2 sn 2 a ξ 0 ) ( cn 2 a η 0 + b 2 sn 2 a η 0 ) ,
z 14 = b 2 ( sn a ξ cn a ξ 0 dn a η 0 sn a ξ 0 cn a ξ dn a η ) 2 ( sn a η cn a η 0 dn a ξ 0 sn a η 0 cn a η dn a ξ ) 2 ( cn 2 a ξ + b 2 sn 2 a ξ ) ( cn 2 a η b 2 sn 2 a η ) × ( cn 2 a ξ 0 + b 2 sn 2 a ξ 0 ) ( cn 2 a η 0 b 2 sn 2 a η 0 ) .
System 2: For the symmetry S 1 + 2 S 2 we find the separable coordinate system
r = tan a ( ξ + η ) 2 , s = tan a ( ξ η ) 2 .
Substituting (123) into (17) yields
U ξ ξ U η η + a 2 m 1 ( 1 m 1 ) sin 2 a ξ m 2 ( 1 m 2 ) sin 2 a η + m 3 ( 1 m 3 ) cos 2 a ξ m 4 ( 1 m 4 ) cos 2 a η U = 0 .
The Riemann function for (124) is then
R ( ξ , η , ξ 0 , η 0 ) = F B m 1 , m 2 , m 3 , m 4 , 1 m 1 , 1 m 2 , 1 m 3 , 1 m 4 , 1 , z 15 , z 16 , z 17 , z 18 ,
where
z 15 = cos a ( ξ ξ 0 ) cos a ( η η 0 ) 2 sin a ξ sin a ξ 0 , z 16 = cos a ( η η 0 ) cos a ( ξ ξ 0 ) 2 sin a η sin a η 0 ,
z 17 = cos a ( ξ ξ 0 ) cos a ( η η 0 ) 2 cos a ξ cos a ξ 0 , z 18 = cos a ( η η 0 ) cos a ( ξ ξ 0 ) 2 cos a η cos a η 0 .
System 3: For the symmetry S 1 2 S 2 , we find the separable coordinate system
r = tanh a ( ξ + η ) 2 , s = tanh a ( ξ η ) 2 .
Substituting (128) into (17) yields
U ξ ξ U η η + a 2 m 1 ( 1 m 1 ) sinh 2 a ξ m 2 ( 1 m 2 ) sinh 2 a η + m 3 ( 1 m 3 ) cosh 2 a η m 4 ( 1 m 4 ) cosh 2 a ξ U = 0 .
The Riemann function for (129) is now
R ( ξ , η , ξ 0 , η 0 ) = F B m 1 , m 2 , m 3 , m 4 , 1 m 1 , 1 m 2 , 1 m 3 , 1 m 4 , 1 , z 19 , z 20 , z 21 , z 22
and
z 19 = cosh a ( η η 0 ) cosh a ( ξ ξ 0 ) 2 sinh a ξ sinh a ξ 0 , z 20 = cosh a ( ξ ξ 0 ) cosh a ( η η 0 ) 2 sinh a η sinh a η 0 ,
z 21 = cosh a ( η η 0 ) cosh a ( ξ ξ 0 ) 2 cosh a η cosh a η 0 , z 22 = cosh a ( ξ ξ 0 ) cosh a ( η η 0 ) 2 cosh a ξ cosh a ξ 0 .
System 4: For the symmetry S = S 1 + 2 q S 2 where q < 1 , the separable coordinate system is
r = b sn a ( ξ + η ) , s = b sn a ( ξ η ) ,
where k = b 2 , 2 q = ( 1 + b 4 ) / b 2 and 0 < b < 1 . As in system 1, a is arbitrary and k is the modulus of the elliptic functions. Substituting (133) into (17) yields
U ξ ξ U η η + { a 2 m 1 ( 1 m 1 ) cn 2 a ξ dn 2 a ξ sn 2 a ξ ( 1 k 2 ) 2 sn 2 a η cn 2 a η dn 2 a η m 2 ( 1 m 2 ) cn 2 a η dn 2 a η sn 2 a η ( 1 k 2 ) 2 sn 2 a ξ cn 2 a ξ dn 2 a ξ + 4 a 2 b 2 m 3 ( 1 m 3 ) cn 2 a ξ dn 2 a ξ ( 1 b 2 sn 2 a ξ ) 2 + cn 2 a η dn 2 a η ( 1 + b 2 sn 2 a η ) 2 1 m 4 ( 1 m 4 ) cn 2 a ξ dn 2 a ξ ( 1 + b 2 s n 2 a ξ ) 2 + cn 2 a η dn 2 a η ( 1 b 2 s n 2 a η ) 2 1 } U = 0 .
The Riemann function for (134) is then
R ( ξ , η , ξ 0 , η 0 ) = F B m 1 , m 2 , m 3 , m 4 , 1 m 1 , 1 m 2 , 1 m 3 , 1 m 4 , 1 , z 27 , z 28 , z 29 , z 30 ,
where
z 27 = sn a η cn a ξ dn a ξ 0 sn a η 0 cn a ξ 0 dn a ξ 2 sn a ξ cn a η dn a η 0 sn a ξ 0 cn a η 0 dn a η 2 ( 4 sn a ξ sn a ξ 0 cn a η cn a η 0 dn a η dn a η 0 ) ,
z 28 = sn a ξ cn a η dn a η 0 sn a ξ 0 cn a η 0 dn a η 2 sn a η cn a ξ dn a ξ 0 sn a η 0 cn a ξ 0 dn a ξ 2 ( 4 sn a η sn a η 0 cn a ξ cn a ξ 0 dn a ξ dn a ξ 0 ) ,
z 29 = b 2 sn a η cn a ξ dn a ξ 0 sn a η 0 cn a ξ 0 dn a ξ 2 sn a ξ cn a η dn a η 0 sn a ξ 0 cn a η 0 dn a η 2 ( 1 b 2 sn 2 a ξ ) ( 1 + b 2 sn 2 a η ) × ( 1 b 2 sn 2 a ξ 0 ) ( 1 + b 2 sn 2 a η 0 ) ,
z 30 = b 2 sn a ξ cn a η dn a η 0 sn a ξ 0 cn a η 0 dn a η 2 sn a η cn a ξ dn a ξ 0 sn a η 0 cn a ξ 0 dn a ξ 2 ( 1 + b 2 sn 2 a ξ ) ( 1 b 2 sn 2 a η ) × ( 1 + b 2 sn 2 a ξ 0 ) ( 1 b 2 sn 2 a η 0 ) .
Returning to (124), setting m 3 = m 4 = 0 , we recover (52). Hence (52) is in fact a sub-case of Chaundy’s equation, obtainable via a Lie-Bäcklund symmetry. A quick calculation shows that (124) only admits one first order symmetry. As such, the group is not isomorphic to S L ( 2 , R ) . In this way, a new equivalence class of Riemann functions has been established. Likewise, take (129) and let m 2 = m 3 = 0 . The Riemann function for this case was first given in [51] by using the techniques from Section 3.3. It is therefore also part of the equivalence class admitting the two-dimensional generalized symmetry algebra (114) and (115).

5. 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
U ξ ξ U η η λ 1 2 m 2 ( 1 m 2 ) sin 2 λ 1 η + m 4 ( 1 m 4 ) cos 2 λ 1 η U = 0 ,
U ξ ξ U η η + λ 2 2 m 1 ( 1 m 1 ) sinh 2 λ 2 ξ m 3 ( 1 m 3 ) cosh 2 λ 2 ξ U = 0 ,
U ξ ξ U η η + n ( 1 n ) ξ 2 U = 0 ,
U ξ ξ U η η + ( 1 + δ ξ ) U = 0 .
The Riemann functions for the first three equations are respectively
R 1 ( ξ , η , ξ 0 , η 0 ) = F 3 m 2 , m 4 , 1 m 2 , 1 m 4 , 1 , z 16 , z 18 ,
R 2 ( ξ , η , ξ 0 , η 0 ) = F 3 m 1 , m 3 , 1 m 1 , 1 m 3 , 1 , z 19 , z 22 ,
R 3 ( ξ , η , ξ 0 , η 0 ) = 2 F 1 n , 1 n ; 1 ; z 0 ,
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. Call R 4 ( ξ , η , ξ 0 η 0 ) the Riemann function for (), which was defined by Equations (86)–(88) (again interchanging the roles of x and y with ξ and η as required).
Start by combining (140) and (141), which produces the equation
U ξ ξ U η η + λ 1 2 m 2 ( 1 m 2 ) sin 2 λ 1 η + m 4 ( 1 m 4 ) cos 2 λ 1 η λ 2 2 m 1 ( 1 m 1 ) sinh 2 λ 2 ξ m 3 ( 1 m 3 ) cosh 2 λ 2 ξ U = 0 .
The Riemann function for (147) is then
R ( ξ , η , ξ 0 , η 0 ) = F 3 ( m 2 , m 4 , 1 m 2 , 1 m 4 , 1 , z 33 , z 34 ) + ξ ξ 0 η η 0 F 3 m 2 , m 4 , 1 m 2 , 1 m 4 , 1 , u 2 ( t ) , u 3 ( t ) × t F 3 m 1 , m 3 , 1 m 1 , 1 m 3 , 1 , v 2 ( t ) , v 3 ( t ) d t ,
where
u 2 ( t ) = cos λ 1 ( η η 0 ) cos λ 1 t 2 sin λ 1 η sin λ 1 η 0 , u 3 ( t ) = cos λ 1 ( η η 0 ) cos λ 1 t 2 cos λ 1 η cos λ 1 η 0 , v 2 ( t ) = cosh λ 2 t cosh λ 2 ( ξ ξ 0 ) 2 sinh λ 2 ξ sinh λ 2 ξ 0 , v 3 ( t ) = cosh λ 2 ( ξ ξ 0 ) cosh λ 2 t 2 cosh λ 2 ξ cosh λ 2 ξ 0 .
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].
The next obvious choice is to combine Equations (142) with (140) and (141) to produce
U ξ ξ U η η + m 1 ( 1 m 1 ) ξ 2 m 2 ( 1 m 2 ) sin 2 η m 4 ( 1 m 4 ) cos 2 η U = 0 ,
U ξ ξ U η η + m 1 ( 1 m 1 ) ξ 2 m 2 ( 1 m 2 ) sinh 2 η + m 3 ( 1 m 3 ) cosh 2 η U = 0 .
The Riemann function for (150) was first given in [51] by using the addition formula. However, both (149) and (150) can be obtained as confluent reductions of (147). To see this, firstly let m 3 ( 1 m 3 ) = 0 in (147). Then if we take the limit as λ 2 0 , we obtain (149). In a completely analogous way, to recover (150), start with (147), set m 4 ( 1 m 4 ) = 0 , and again take the limit as λ 1 0 . The result is (150) as desired.
In order to obtain truly new Riemann functions, Equations (140)–(142) need to be combined with (143). Hence, for the equations
U ξ ξ U η η m 2 ( 1 m 2 ) sin 2 η + m 4 ( 1 m 4 ) cos 2 η + ( 1 + δ ξ ) U = 0 ,
U ξ ξ U η η m 2 ( 1 m 2 ) sinh 2 η m 3 ( 1 m 3 ) cosh 2 η + ( 1 + δ ξ ) U = 0 ,
U ξ ξ U η η m 1 ( 1 m 1 ) η 2 + ( 1 + δ ξ ) U = 0 ,
the Riemann functions are respectively
R ( ξ , η , ξ 0 , η 0 ) = R 1 ( ξ , η , ξ 0 , η 0 ) + ξ ξ 0 η η 0 R 1 ( t , η , 0 , η 0 ) t R 4 ( ξ , t , ξ 0 , 0 ) d t , R ( ξ , η , ξ 0 , η 0 ) = R 2 ( ξ , η , ξ 0 , η 0 ) + ξ ξ 0 η η 0 R 2 ( t , η , 0 , η 0 ) t R 4 ( ξ , t , ξ 0 , 0 ) d t , R ( ξ , η , ξ 0 , η 0 ) = R 3 ( ξ , η , ξ 0 , η 0 ) + ξ ξ 0 η η 0 R 3 ( t , η , 0 , η 0 ) t R 4 ( ξ , t , ξ 0 , 0 ) d t ,
where R 1 ( ξ , η , ξ 0 , η 0 ) , R 2 ( ξ , η , ξ 0 , η 0 ) , and R 3 ( ξ , η , ξ 0 , η 0 ) are given by (144)–(146). R 4 ( ξ , η , ξ 0 , η 0 ) is defined by (86)–(88) and the roles of ξ and η are interchanged as required. The final possibility is to use (143) twice. That is, to take the equation
U ξ ξ U η η ( δ ξ ω η ) U = 0 .
The Riemann function for (154) is given by
R ( ξ , η , ξ 0 , η 0 ) = R 4 ( ξ , η , ξ 0 , η 0 ) + ξ ξ 0 η η 0 R 5 ( t , η , 0 , η 0 ) t R 4 ( ξ , t , ξ 0 , 0 ) d t ,
where R 4 ( ξ , η , ξ 0 , η 0 ) remains unchanged. R 5 is obtained by interchanging the roles of ξ and η in (86)–(88). The Riemann function for (154) first appeared in [69] for the case where δ = ω . The full equation has not been published before.
Equations (151) to (154) admit no point symmetries. Likewise, they have no non-trivial Lie-Bäcklund symmetries. Hence, they fall outside the equivalence classes of Riemann functions studied in Section 4.1 and Section 4.5. In this way, (151)–(154) represent a completely new class of Riemann functions.
In principle, all the Equations (140)–(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.

6. Conclusions

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.
  • The addition formula.
  • Laplace transforms.
  • The multiplication formula.
  • Finite groups.
  • Lie-Bäcklund symmetries.
Furthermore, the equivalences generated by point symmetries, where the governing equation admits a symmetry algebra isomorphic to S L ( 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,73,74]. It is largely based on the work of Vekua [75]. Qualitative properties of the Riemann function have also been explored [76,77,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].

Funding

This research was funded by an Australian Postgraduate Award.

Acknowledgments

The referees are thanked for their feedback. The author is also grateful to Edward D. Fackerell and Christopher M. Cosgrove for their guidance on this work.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Hypergeometric Functions

Here, the definitions and key properties of the hypergeometric functions seen in this paper are included for reference. Greater detail can be found in [33,36]. The functions are included in the order in which they appeared.
  • 2 F 1 ( a , b ; c ; z )
    2 F 1 ( a , b ; c ; z ) = n = 0 ( a ) n ( b ) n z n ( c ) n n ! ,
    where
    ( a ) n = Γ ( a + n ) Γ ( a ) , ( a ) 0 = 1 , ( a ) n = a ( a + 1 ) ( a + n 1 ) , n = 1 , 2 , 3 ,
    This function is a particular solution of the equation
    z ( 1 z ) d 2 U d z 2 + c ( a + b + 1 ) z d U d z a b U = 0 .
    Of particular note is the relationship between the Legendre function and the hypergeometric function
    P m ( 1 2 z ) = 2 F 1 ( m , 1 m , 1 , z ) .
  • Ξ 2 ( a , b , c , z 1 , z 2 )
    Ξ 2 ( a , b , c , z 1 , z 2 ) = ( a ) m ( b ) n ( c ) m + n m ! n ! z 1 m z 2 n ,
    where ( a ) m etc. are given by A1). This function is the solution to the system of equations
    z 1 ( 1 z 1 ) 2 U z 1 2 + z 2 2 U z 1 z 2 + c ( a + b + 1 ) z 1 U z 1 a b U = 0 , z 2 2 U z 2 2 + z 1 2 U z 1 z 2 + c U z 2 U = 0 .
  • F 3 ( a , a , b , b , c , z 1 , z 2 )
    F 3 ( a , a , b , b , c , z 1 , z 2 ) = ( a ) m ( a ) n ( b ) m ( b ) n ( c ) m + n m ! n ! z 1 m z 2 n ,
    where ( a ) m etc. are given by (A1). This function is the solution to the system of equations
    z 1 ( 1 z 1 ) 2 U z 1 2 + z 2 2 U z 1 z 2 + c ( a + b + 1 ) z 1 U z 1 a b U = 0 , z 2 ( 1 z 2 ) 2 U z 2 2 + z 1 2 U z 1 z 2 + c ( a + b + 1 ) z 2 U z 2 a b U = 0 .
  • F B ( a 1 , , a n , b 1 , , b n , c , z 1 , , z n )
    F B ( a 1 , , a n , b 1 , , b n , c , z 1 , , z n ) = ( a 1 ) m 1 ( a n ) m n ( b 1 ) m 1 ( b n ) m n ( c ) m 1 + + m n m 1 ! m n ! z 1 m 1 z n m n .
    where ( a 1 ) m 1 etc. are given by (A1). This function satisfies the system of equations
    z j ( 1 z j ) 2 U z j 2 + k j z k 2 U 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 )
    J n ( z ) = r = 0 ( 1 ) n z n + 2 r 2 r ! ( n + r ) ! .
    This function is a particular solution of the equation
    d 2 U d z 2 + 1 z d U d z + 1 n 2 z 2 U = 0 .
  • F 4 ( a , b , c , c , z 1 , z 2 )
    F 4 ( a , b , c , c , z 1 , z 2 ) = ( a ) m + n ( b ) m + n ( c ) m ( c ) n m ! n ! z 1 m z 2 n ,
    where ( a ) m + n etc. are given by (A1). This function is the solution of the system of equations
    z 1 ( 1 z 1 ) 2 U z 1 2 z 2 2 2 U z 2 2 2 z 1 z 2 2 U z 1 z 2 + c ( a + b + a ) z 1 U z 1 ( a + b + 1 ) z 2 U z 2 a b U = 0 , z 2 ( 1 z 2 ) 2 U z 2 2 z 1 2 2 U z 1 2 2 z 1 z 2 2 U z 1 z 2 + c ( a + b + a ) z 2 U z 1 ( a + b + 1 ) z 1 U z 1 a b U = 0 .

References

  1. Solimeno, S. Riemann-Green’s Functions for Solving Electromagnetic Problems Exhibiting Rotational Symmetry in Media Moving with Superluminal Velocities. J. Math. Phys. 1975, 16, 218–223. [Google Scholar] [CrossRef]
  2. Iraniparast, N. Green-Riemann Functions for a Class of Hyperbolic Focal Point Problems. SIAM J. Math. Anal. 1989, 20, 408–414. [Google Scholar] [CrossRef]
  3. Asfar, O.R. Riemann-Green Function Solution of Transient Plane Waves in Lossy Media. IEEE Trans. Electromagn. Compat. 1990, 32, 228–231. [Google Scholar] [CrossRef]
  4. Coz, M.; Coudray, C. Existence of Generalized Translation Operators from the Agranovitch-Marchenko Transformation (Jost Solutions). J. Math. Phys. 1973, 14, 1574–1578. [Google Scholar] [CrossRef]
  5. Coz, M.; Coudray, C. The Riemann Solution and the Inverse Quantum Mechanical Problem. J. Math. Phys. 1976, 17, 888–893. [Google Scholar] [CrossRef]
  6. Coz, M. The Riemann Solution in the One-Dimensional Inverse Problem. J. Math. Anal. Appl. 1977, 61, 32–250. [Google Scholar] [CrossRef]
  7. Coz, M.; Rochus, P. Partial Differential Matrix Equations for the Inverse Problem of Scattering Theory. J. Math. Phys. 1976, 17, 894–899. [Google Scholar] [CrossRef]
  8. Coz, M.; Rochus, P. The Translation Kernel in the N-Dimensional Scattering Problem. J. Math. Phys. 1977, 18, 2223–2231. [Google Scholar] [CrossRef]
  9. Coz, M.; Rochus, P. Translation Kernels for Velocity Dependent Interactions. J. Math. Phys. 1977, 18, 2232–2240. [Google Scholar] [CrossRef]
  10. Coz, M.; Rochus, P. Extension of the Marchenko Equation to Non-Hermitian Differential Systems. Anna. Phys. 1980, 126, 460–499. [Google Scholar] [CrossRef]
  11. Gugushvili, H.I.; Mentkovsky, Y.L. Inverse Problem of the Theory of Charged Particle Scattering. Ukr. Phys. J. 1969, 13, 1252–1263. [Google Scholar]
  12. Gugushvili, H.I.; Mentkovsky, Y.L. The Inverse Problem of Scattering Theory and Riemann’s Method. Il Nuovo Cimento 1972, 10, 277–291. [Google Scholar] [CrossRef]
  13. Matsuno, Y. Reduction of Dispersionless Coupled Korteweg-de Vries Equations to the Euler-Darboux Equation. J. Math. Phys. 2001, 42, 1744–1760. [Google Scholar] [CrossRef]
  14. Popivanov, N.; Popov, T. Exact Behaviour of Singularities of Protter’s Problem for the 3-D Wave Equation. In Inclusion Methods for Nonlinear Problems; Herzberger, J., Ed.; Springer: New York, NY, USA, 2003; pp. 213–236. [Google Scholar]
  15. Popivanov, N.; Popov, T.; Scherer, R. Asymptotic Expansions of Singular Solutions for (3 + 1)-D Protter Problems. J. Math. Anal. Appl. 2007, 331, 1093–1112. [Google Scholar] [CrossRef]
  16. Popivanov, N.; Popov, T.; Scherer, R. Protter-Morawetz Multidimensional Problems. Proc. Steklov Inst. Math. 2012, 278, 179–198. [Google Scholar] [CrossRef]
  17. Dechevski, L.; Popivanov, N.; Popov, T. Exact Asymptotic Expansion of Singular Solutions for the (2 + 1)-D Protter Problem. Abstr. Appl. Anal. 2012. [Google Scholar] [CrossRef]
  18. Popivanov, N.; Popov, T.; Tesdall, A. Semi-Fredholm Solvability in the Framework of Singular Solutions for the (3 + 1)-D Protter-Morawetz Problem. Abstr. Appl. Anal. 2014. [Google Scholar] [CrossRef]
  19. Wright, O.C.; Forest, M.G.; McLaughlin, K.T.-R. On the Exact Solution of the Geometric Optics Approximation of the Defocusing Non-Linear Schrödinger Equation. Phys. Lett. A 1999, 257, 170–174. [Google Scholar] [CrossRef]
  20. Forest, M.G.; Rosenberg, C.J.; Wright, O.C. On the Exact Solution for Smooth Pulses of the Defocusing nonlinear Schrödinger Modulation Equations Prior to Breaking. Nonlinearity 2009, 9, 2287–2308. [Google Scholar] [CrossRef]
  21. Shiryaeva, E.V.; Zhukov, M.Y. Hodograph Method and Numerical Integration of Two Hyperbolic Quasilinear Equations. Part I. The Shallow Water Equations; Working Paper; Southern Federal University: Rostov-on-Don, Russia, 2014. [Google Scholar]
  22. Shiryaeva, E.V.; Zhukov, M.Y. Hodograph Method and Numerical Solution of the Two Hyperbolic Quasilinear Equations System. Part II. Zonal Electrophoresis Equations; Working Paper; Southern Federal University: Rostov-on-Don, Russia, 2014. [Google Scholar]
  23. Shiryaeva, E.V.; Zhukov, M.Y. Hodograph Method and Numerical Solution of the Two Hyperbolic Quasilinear Equations System. Part III. Two Beam Reduction of the Dense Soliton Gas Equations; Working Paper; Southern Federal University: Rostov-on-Don, Russia, 2014. [Google Scholar]
  24. Copson, E.T. On the Riemann-Green Function. Arch. Rat. Mech. Anal. 1958, 1, 324–348. [Google Scholar] [CrossRef]
  25. Riemann, B. Collected Works of Bernhard Riemann; Dover: New York, NY, USA, 1953; pp. 156–175. [Google Scholar]
  26. Hadamard, J. Lectures on Cauchy’s Problem in Linear Partial Differential Equations; Dover: New York, NY, USA, 1952. [Google Scholar]
  27. Chaundy, T.W. Hypergeometric Partial Differential Equations (I). Q. J. Math. Oxf. Ser. 1935, 6, 288–303. [Google Scholar] [CrossRef]
  28. Chaundy, T.W. Hypergeometric Partial Differential Equations (II). Q. J. Math. Oxf. Ser. 1936, 7, 306–315. [Google Scholar] [CrossRef]
  29. Chaundy, T.W. Partial Differential Equations With Constant Coefficients (II). Q. J. Math. Oxf. Ser. 1937, 8, 280–302. [Google Scholar] [CrossRef]
  30. Chaundy, T.W. Linear Partial Differential Equations (I). Q. J. Math. Oxf. Ser. 1938, 9, 234–240. [Google Scholar] [CrossRef]
  31. Chaundy, T.W. Hypergeometric Partial Differential Equations (III). Q. J. Math. Oxf. Ser. 1939, 10, 219–240. [Google Scholar] [CrossRef]
  32. Chaundy, T.W. Linear Partial Differential Equations (II). Q. J. Math. Oxf. Ser. 1940, 11, 101–110. [Google Scholar] [CrossRef]
  33. Appell, P.; Kampé de Fériet, J. Fonctions Hypergéométriques et Hypersphériques, Polynomes D’Hermite; Gauthier-Villars: Paris, France, 1926. [Google Scholar]
  34. Titchmarsh, E.C. Theory of Fourier Integrals; Oxford University Press: London, UK, 1937; pp. 297–298. [Google Scholar]
  35. Lanckau, E. Die Riemannfunktion Selbstadjungierter Gleichungen. Wiss. Z. Tech. Hochsch. Karl-Marx-Stadt 1979, 21, 535–540. [Google Scholar]
  36. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions I; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  37. Henrici, P.A. Survey of I. N. Vekua’s Theory of Elliptic Partial Differential Equations with Analytic Coefficients. Z. Angew. Math. Phys. 1957, 8, 169–202. [Google Scholar] [CrossRef]
  38. Püngel, J. Zur Darstellung Von Riemannfunktionen Durch Integro Differentialoperatoren. Sitzungsberichte. Abt. 2. Methematik, Astronomie, Physik, Meteorologie Und Technik 1980, 189, 55–66. [Google Scholar]
  39. Püngel, J. Lineare Abbidungen Zwischen Lösungsmengen Partieller Differentialgleichungen Im Komplexen. Graz Mathematische-Statistiche Sektion; Forschungsgesellschaft Joanneum: Graz, Austria, 1978. [Google Scholar]
  40. Courant, R.; Hilbert, D. Methods of Mathematical Physics II; Interscience Publishers: New York, NY, USA, 1962. [Google Scholar]
  41. Cohn, H. The Riemann function for Uxy + H(x + y)U = 0. Duke Math. J. 1947, 14, 297–304. [Google Scholar] [CrossRef]
  42. Olevskiĭ, M.N. On Riemann’s Function for the Differential Equation (2u/x2) − (2u/t2) + [ρ1(x) + ρ2(t)]u = 0. Dokl. Akad. Nauk SSSR 1952, 87, 337–340. [Google Scholar]
  43. Lie, S. On Integration of a Class of Linear Partial Differential Equations by Means of Definite Integrals. Arch. Math. Bd VI Heft 1881, 3, 328–368. [Google Scholar]
  44. Bluman, G.W. Construction of Solutions to Partial Differential Equations by the Use of Transformation Groups. Ph.D. Thesis, California Institute Of Technology, Pasadena, CA, USA, 1967. [Google Scholar]
  45. Daggit, E.A. The Use of Infinitesimal Transformations in Predicting the Form of the Riemann (-Green) Function. J. Math. Anal. Appl. 1970, 29, 91–108. [Google Scholar] [CrossRef]
  46. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations, Volume 2, Applications in Engineering and Physical Sciences; CRC Press: Boca Raton, FL, USA, 1995; pp. 41–51. [Google Scholar]
  47. Iwasaki, K. Riemann Function of Harmonic Equation and Appell’s F4. SIAM J. Math. Anal. 1988, 19, 903–917. [Google Scholar] [CrossRef]
  48. Hille, E. Lectures on Ordinary Differential Equations; Addison Wesley: Reading, MA, USA, 1969. [Google Scholar]
  49. Vaz, P.T.; Deo, S.G. On a Class of Linear Hyperbolic Differential Equations. Indian J. Pure Appl. Math. 1987, 9, 801–809. [Google Scholar]
  50. Papadakis, J.S.; Wood, D.H. Addendum: “An Addition Formula for Riemann Functions”. J. Differ. Equ. 1978, 29, 474. [Google Scholar]
  51. Papadakis, J.S.; Wood, D.H. An Addition Formula for Riemann Functions. J. Differ. Equ. 1977, 24, 397–411. [Google Scholar] [CrossRef]
  52. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: New York, NY, USA, 1996. [Google Scholar]
  53. Kokinasidi, P.D. Construction of the Riemann Function for Some Hyperbolic Equations. Differ. Uravn. 1983, 19, 703–706. [Google Scholar]
  54. Johnpillai, I.K.; Mahomed, F.M. Goursat Problem for the Factorizable Hyperbolic Equation in two Independent Variables. Math. Comp. Appl. 2003, 8, 55–62. [Google Scholar] [CrossRef]
  55. Andrey, A.; Yevgeniya, K.; Anzhelika, V. Construction of the Solution of the Caushy’s Problem by the Riemann’s Method for a Hyperbolic Equation. Am. Res. J. Math. 2015, 1, 44–48. [Google Scholar]
  56. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  57. Mackie, A.G. Green’s Functions and Riemann’s Method. Proc. Edinb. Math. Soc. 1964, 14, 293–302. [Google Scholar] [CrossRef]
  58. Wood, D.H. Simple Riemann Functions. Bull. Am. Math. Soc. 1976, 82, 737–739. [Google Scholar] [CrossRef]
  59. Clarkson, P.A.; Kruskal, M.D. New Similarity Reductions of the Boussinesq Equation. J. Math. Phys. 1989, 30, 2201–2213. [Google Scholar] [CrossRef]
  60. Geddes, R.L.; Mackie, A.G. Riemann Functions for Self-adjoint Equations. Appl. Anal. 1977, 7, 43–47. [Google Scholar] [CrossRef]
  61. Bluman, G.W. On Mapping Linear Partial Differential Equations to Constant Coefficient Equations. SIAM J. Appl. Math. 1983, 43, 1259–1273. [Google Scholar] [CrossRef]
  62. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982; pp. 105–117. [Google Scholar]
  63. Miller, W., Jr. Lie Theory and Special Functions; Academic: New York, NY, USA, 1968. [Google Scholar]
  64. Scott, E.J. Determination of the Riemann Function. Am. Math. Mon. 1973, 80, 907–909. [Google Scholar] [CrossRef]
  65. Wahlberg, C. Riemann’s Function for a Klein-Gordon Equation with a Non-Constant Coefficient. J. Phys. A Math. Gen. 1977, 10, 867–878. [Google Scholar] [CrossRef]
  66. Ince, E.L. Ordinary Differential Equations; Dover: New York, NY, USA, 1956; p. 466. [Google Scholar]
  67. Du, X.H. A Note on Olevskii’s Formula for Self-Adjoint Hyperbolic Equations of the Second Order. Kexue Tongbao 1981, 26, 869–873. [Google Scholar]
  68. Zeitsch, P.J. Symmetry Operators for Riemann’s Method. J. Math. Phys. 2004, 45, 2993–3000. [Google Scholar] [CrossRef]
  69. Zeitsch, P.J. Riemann Functions and the Group E(1, 1). J. Math. Phys. 2004, 45, 3001–3018. [Google Scholar] [CrossRef]
  70. Miller, W., Jr. Symmetry and Separation of Variables; Addison-Wesley: London, UK, 1977. [Google Scholar]
  71. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1965. [Google Scholar]
  72. Bauer, K.W. On a Class of Riemann Functions. Appl. Anal. 1982, 13, 109–126. [Google Scholar] [CrossRef]
  73. Bauer, K.W. On the Determination of Riemann Functions. Complex Var. Theory Appl. 1987, 8, 195–203. [Google Scholar] [CrossRef]
  74. Markett, C.; Püngell, J.; Wallner, H. Multiple Power Series Representations for Riemann Functions of Self-Adjoint Equations. Appl. Anal. 1990, 38, 179–199. [Google Scholar] [CrossRef]
  75. Vekua, I.N. New Methods for Solving Elliptic Equations; North Holland Publish. Co.: Amsterdam, The Netherlands, 1968. [Google Scholar]
  76. Lerner, M.E. Qualitative Properties of the Riemann Function. Differ. Equ. 1992, 27, 1495–1508. [Google Scholar]
  77. Wallner, H. Die Riemannfunktion der Differentialgleichung W+[z+φ(ζ)]W = 0. Arch. Math. 1981, 37, 435–442. [Google Scholar] [CrossRef]
  78. Zhegalov, V.I.; Kotukhov, M.P. On Integral Equations for Riemann Function. Russ. Math. (Iz. VUZ) 1998, 42, 24–28. [Google Scholar]
  79. Zhegalov, V.I. On the three-dimensional Riemann Function. Sib. Math. J. 1998, 38, 929–934. [Google Scholar] [CrossRef]
  80. Koshcheeva, O.A. Construction of the Riemann Function for the Bianchi Equation in an n-Dimensional Space. Russ. Math. (Iz. VUZ) 2008, 52, 35–40. [Google Scholar] [CrossRef]
Figure 1. The Initial Curve.
Figure 1. The Initial Curve.
Mathematics 06 00316 g001

Share and Cite

MDPI and ACS Style

Zeitsch, P.J. On the Riemann Function. Mathematics 2018, 6, 316. https://doi.org/10.3390/math6120316

AMA Style

Zeitsch PJ. On the Riemann Function. Mathematics. 2018; 6(12):316. https://doi.org/10.3390/math6120316

Chicago/Turabian Style

Zeitsch, Peter J. 2018. "On the Riemann Function" Mathematics 6, no. 12: 316. https://doi.org/10.3390/math6120316

APA Style

Zeitsch, P. J. (2018). On the Riemann Function. Mathematics, 6(12), 316. https://doi.org/10.3390/math6120316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop