arxiv: v2 [math.ca] 19 Oct 2012

Transcription

1 Symmetry, Integrability and Geometry: Methods and Applications Def inite Integrals using Orthogonality and Integral Transforms SIGMA 8 (, 77, pages Howard S. COHL and Hans VOLKMER arxiv:.4v [math.ca] 9 Oct Applied and Computational Mathematics Division, Information Technology Laboratory, National Institute of Standards and Technology, Gaithersburg, MD, , USA howard.cohl@nist.gov URL: Department of Mathematical Sciences, University of Wisconsin-Milwaukee, P.O. Box 43, Milwaukee, WI, 53, USA volkmer@uwm.edu Received July 3,, in final form October 5, ; Published online October 9, Abstract. We obtain definite integrals for products of associated Legendre functions with Bessel functions, associated Legendre functions, and Chebyshev polynomials of the first kind using orthogonality and integral transforms. Key words: definite integrals; associated Legendre functions; Bessel functions; Chebyshev polynomials of the first kind Mathematics Subject Classification: 6A4; 33C5; 33C; 33C45; 35A8 Introduction In [3] and [6] (see also [4], we present some definite integral and infinite series addition theorems which arise from expanding fundamental solutions of elliptic equations on R d in axisymmetric coordinate systems which separate Laplace s equation. We utilize orthogonality and integral transforms to obtain new definite integrals from some of these addition theorems. Def inite integrals from integral transforms. Application of Hankel s transform We use the following result for x (, we define F (r ± := lim x r± F (x; see [5, p. 456]: Theorem. Let F : (, C be such that x F (x dx <, ( and let ν. Then (F (r + + F (r = uj ν (ur xf (xj ν (ux dx du ( provided that the positive number r lies inside an interval in which F (x has finite variation. This paper is a contribution to the Special Issue Superintegrability, Exact Solvability, and Special Functions. The full collection is available at

2 H.S. Cohl and H. Volkmer As an illustration for the method of integral transforms, we give the following example. According to [5, (3..] (see also [7, (6.6.3], we have for Re a >, b, c >, Re ν >, then e ka J ν (kbj ν (kcdk = π bc Q ν / ( a + b + c bc, (3 J ν : C \ (, ] C, for order ν C is the Bessel function of the first kind defined in [, (..] and Q µ ν : C \ (, ] C for ν + µ / N with degree ν and order µ, is the associated Legendre function of the second kind defined in [, (4.3.7, 4.]. The Legendre function of the second kind Q ν : C\(, ] C for ν / N is defined in terms of the zero-order associated Legendre function of the second kind Q ν (z := Q ν(z. If we apply Theorem to the function F : (, C defined by F (k := π c k e ka J ν (kc, (4 then condition ( is satisfied. If we use (4 in ( then we obtain the following result. If Re a >, c >, Re ν >, then J ν (kb Q ν / ( a + b + c bc b π c db = k e ka J ν (kc, which is actually given in [3, (.8.8.]. Hardy [8, (33.6] derives an interesting extension of (3 (see also [5, p. 389] and [, p. 7]. We apply the Whipple formula [, (4.9.7, 4.] to Hardy s extension to obtain ke ka J ν (kbj ν (kcdk ( a a = π + b + c bc (a + (b + c / (a + (b c / Q ν /, (5 bc for Re a >, b, c >, Re ν >. It is mentioned in [5, p. 389] that (5 can be derived from (3 by differentiation with respect to a. Using this integral and Theorem, we prove the following theorem. Theorem. Let Re a >, c >, Re ν >. Then ( J ν (kb a + b + c b π c (a + (b + c / (a + (b c / Q ν / db = bc a e ka J ν (kc. Proof. By applying Theorem to the function F : (, C defined by F (k := π c a e ka J ν (kc, using (5, we obtain the desired result. Now, we give another example of how an integral expansion for a fundamental solution of Laplace s equation on R 3 in parabolic coordinates can be used to prove a new definite integral. Theorem 3. Let m N, λ (,, µ, µ (,, µ µ, k (,. Then Q m / (χj m (kλ λ dλ = π λ µµ J m (kλ I m (kµ < K m (kµ >,

3 Definite Integrals using Orthogonality and Integral Transforms 3 χ = 4λ µ + 4λ µ + (λ λ + µ µ 8λλ µµ > and µ := min max {µ, µ }. Proof. We apply Theorem to the function F : (, C defined by F (k := π λ µµ J m (kλ I m (kµ < K m (kµ >, I ν : C \ (, ] C, for ν C, is the modified Bessel function of the first kind defined in [, (.5.] and K ν : C \ (, ] C, for ν C, is the modified Bessel function of the second kind defined in [, (.7.4]. Again, we see that condition ( is satisfied. We obtain the desired result from [6, (4], namely, for λ (,, J m (kλj m (kλ I m (kµ < K m (kµ > kdk = Q m /(χ π λλ µµ.. Application of Fourier cosine transform Theorem 4. Let a, b (, with b a, k (,, Re ν >. Then Q ν / ( a + b + z ab cos(kz dz = π ab I ν (kbk ν (ka. (6 Proof. According to [7, (6.67.4] we have the integral relation I ν (kbk ν (ka cos(kzdk = ( a ab Q + b + z ν /, ab a, b (, with b < a, z >, Re ν >. We obtain the desired result from Theorem with F : (, C defined such that ab F (k := π k I ν(kbk ν (ka and ν =. Furthermore, if one makes the replacement z z/( a, k ka in [7, (7.6.6], namely Q ν / ( + z cos(kzdz = π ( ( k k I ν K ν, k (, and Re ν >, then we see that (6 holds also for any a, b (, with a = b. 3 Def inite integrals from orthogonality relations 3. Degree orthogonality for associated Legendre functions with integer degree and order We take advantage of the degree orthogonality relation for the Ferrers function of the first kind with integer degree and order, namely (cf. [7, (7..] P m n (cos θp m n (cos θ sin θdθ = n + (n + m! (n m! δ n,n, (7

4 4 H.S. Cohl and H. Volkmer m, n, n N, and m n, m n. We are using the associated Legendre function of the first kind (on-the-cut, P µ ν : (, C, for ν, µ C, the Ferrers function of the first kind, which is defined in [, (4.3.]. The following estimates for the Ferrers function of the first kind with integer degree and order will be useful. If θ [, π] and m, n N then [4, 5.3, (9] P m n (cos θ (m + n! n! and if θ (, π then [, p. 3] P m (n + m! n (cos θ < (πn / (csc θ m+/. (9 n! If µ C, ξ > are fixed and ν +, we also have the following asymptotic formulas for the associated Legendre functions P ν µ / Γ(ν + µ + (cosh ξ = (π sinh ξ Γ(ν + 3 ( π / Q µ Γ(ν + µ + ν (cosh ξ = sinh ξ Γ(ν + 3 P ν µ / Γ(ν + µ + (i sinh ξ = (π cosh ξ Γ(ν + 3 ( π / Q µ Γ(ν + µ + ν (i sinh ξ = cosh ξ Γ(ν + 3 (8 e (ν+ ξ( + O ( ν, ( e (ν+ ξ+iπµ( + O ( ν, ( e (ν+ ξ+iπν/( + O ( ν, ( e (ν+ ξ iπ(ν+/+iπµ( + O ( ν. (3 These asymptotic formulae follow from representations of Legendre functions by Gauss hypergeometric functions; see [, (8.., (8..4 5]. Theorem 5. Let n, m N, with n m, ν C \ {m, m +, m + 4,...}, r, r (,, r r, θ (, π. Then ( χ (ν+/4 Q (ν+/ m / = i π (ν+/ (sin θ ν/ (χp m n (cos θ(sin θ (ν+/ dθ (ν+/ ( Q (ν+/ r + r n rr P m n (cos θ, (4 ( r > r < rr χ = r + r rr cos θ cos θ rr sin θ sin θ (5 and r := min max {r, r }. Proof. We start with the following addition theorem for the associated Legendre function of the second kind (see [3], namely for θ (, π, ( r (ν+3/ (sin θ ν/ > r< (ν+/ rr ( r + r (n + m! Q (ν+/ n rr P m n (cos θp m n (cos θ, (6 ( χ (ν+/4 (sin θ ν/ Q (ν+/ m / (χ = i π n=m (n m! (n + χ > is given by (5. By (8 and ( the infinite series is uniformly convergent for θ [, π]. Therefore, if we multiply both sides of (6 by sin θ P m n (cos θ, n N and integrate over θ (, π we obtain (4.

5 Definite Integrals using Orthogonality and Integral Transforms 5 Corollary. Let n, m N with n m, r, r (,, r r, θ (, π. Then Q m / (χp m n (cos θ sin θdθ = π sin θ n + Pm n (cos θ χ > is given by (5. ( r< r > n+/, Proof. Substitute ν = in (4 and use [, (4.5.7]. Theorem 6. Let m, n N with n m, σ, σ (,, θ (, π. Then Q m / (χp m n (cos θ m (n m! sin θ dθ = π( (n + m! sinh σ sinh σ sin θ P m n (cos θ P m n (cosh σ < Q m n (cosh σ >, (7 χ = cosh σ + cosh σ sin θ sin θ cosh σ cosh σ cos θ cos θ sinh σ sinh σ sin θ sin θ (8 and σ := min max {σ, σ }. Proof. We start with the following addition theorem for the associated Legendre function of the second kind (see [6, (37], namely Q m / (χ = π( m sinh σ sinh σ sin θ sin θ n=m (n + [ ] (n m! (n + m! P m n (cos θp m n (cos θ P m n (cosh σ < Q m n (cosh σ >, (9 χ is defined by (8. Note that χ = only if σ = σ and θ = θ, and in that case Q m / (χ has a logarithmic singularity. If σ σ then (8, (, ( show that the series in (9 is uniformly convergent for θ [, π]. Therefore, if we multiply both sides of (9 by sin θ P m n (cos θ and integrate over θ [, π], then by (7 we have obtained (7. If σ = σ then one may use (9, (, ( and the orthogonality relation (7 to show that the series in (9 (as a series of functions in the variable θ (, π converges in L (, π. That is n=m { } (n m! (n + m! φ n(θ Pn m (cosh σ < Q m n (cosh σ > <, ( φ n L (, π forms an orthonormal basis and is defined as φ n (θ := n + (n m! sin θ (n + m! Pm n (cos θ, for n = m, m +,.... Then by the asymptotics of P m n and Q m n (cf. (, ( P m n (cosh σ < Q m n (cosh σ > = O ( n m as n. Also from the estimate of P m n (9, φ n (θ = O (. Therefore, (n m! (n + m! φ n(θ P m n (cosh σ < Q m n (cosh σ > = O ( n, and this implies ( because n= n <. Therefore, we again obtain (7.

6 6 H.S. Cohl and H. Volkmer Theorem 7. Let m, n N with m n, σ, σ (,, θ [, π]. Then Q m / (χp m n (cos θ m (n m! sin θdθ = πi( (n + m! cosh σ cosh σ sin θ P m n (cos θ P m n (i sinh σ < Q m n (i sinh σ >, ( χ = sinh σ + sinh σ + sin θ + sin θ sinh σ sinh σ cos θ cos θ cosh σ cosh σ sin θ sin θ. ( Proof. We start with oblate spheroidal coordinates on R 3, namely x = a cosh σ sin θ cos φ, y = a cosh σ sin θ sin φ, z = a sinh σ cos θ, a >, σ [,, θ [, π], φ [, π. The reciprocal distance between two points x, x R 3 expanded in terms of the separable harmonics in this coordinate system is given in [9, (4, p. 8], namely x x = i a (n + n= n [ ] (n m! ( m ɛ m cos(m(φ φ (n + m! m= P m n (cos θp m n (cos θ P m n (i sinh σ < Q m n (i sinh σ >, ɛ m := δ m, is the Neumann factor [, p. 744] commonly occurring in Fourier cosine series, with σ [,, θ [, π], φ [, π. Note that the corresponding expression given in [6, 5.] is given incorrectly (see [4]. By reversing the order of summations in the above expression and comparing with the Fourier cosine expansion for the reciprocal distance between two points, namely x x = πa cosh σ cosh σ sin θ sin θ ɛ m cos(m(φ φ Q m / (χ, m= χ > is given by (, we obtain the following addition theorem for the associated Legendre function of the second kind Q m / (χ = iπ( m cosh σ cosh σ sin θ sin θ n=m (n + [ ] (n m! (n + m! P m n (cos θp m n (cos θ P m n (i sinh σ < Q m n (i sinh σ >. (3 If we multiply both sides of (3 by sin θ P m n (cos θ and integrate over θ [, π], then by (7 we have obtained (. We justify the interchange of integral and infinite sum as before by using the asymptotic formulas (, (3. Theorem 8. Let m, n N with m n, σ, σ (,, θ (, π. Then Q m / (χp m n (cos θ sin θdθ = π sin θ n + Pm n (cos θ e (n+/(σ> σ<, (4 if we define s = cosh σ, s = cosh σ, τ = cos θ, τ = cos θ, then χ = sin θ(s τ + sin θ (s τ + [ (s τ sinh σ (s τ sinh σ ] sin θ sin θ (s τ(s τ. (5

7 Definite Integrals using Orthogonality and Integral Transforms 7 Proof. We start with bispherical coordinates on R 3, namely x = a sin θ cos φ a sin θ sin φ, y = cosh σ cos θ cosh σ cos θ, z = a sinh σ cosh σ cos θ, a >, σ [,, θ [, π], φ [, π. The reciprocal distance between two points x, x R 3 expanded in terms of the separable harmonics in this coordinate system is given in [9, (9, p. ], namely x x = (cosh σ cos θ(cosh σ a cos θ n= n (n m! ɛ m (n + m! Pm n (cos θp m n (cos θ cos(m(φ φ, m= e (n+/(σ> σ< σ [,, θ [, π], φ [, π. By reversing the order of summations in the above expression and comparing with the Fourier cosine expansion for the reciprocal distance between two points, namely (cosh σ cos θ(cosh σ x x = cos θ πa ɛ m cos(m(φ φ Q sin θ sin θ m / (χ, χ is given by (5, we obtain the following addition theorem for the associated Legendre function of the second kind Q m / (χ = π sin θ sin θ n=m m= (n m! (n + m! e (n+/(σ> σ< P m n (cos θp m n (cos θ. (6 If σ = σ and θ = θ, then χ =, and Q m / (χ has a logarithmic singularity. Note that the corresponding expression given in [6, 6., (45] is given incorrectly (see [4]. If we multiply both sides of (6 by sin θ P m n (cos θ and integrate over θ [, π], then by (7 we have obtained (4. We justify the interchange of integral and infinite sum in the same way as in the proof of Theorem Order orthogonality for associated Legendre functions with integer degree and order In this subsection we take advantage of the order orthogonality relation for the Ferrers function of the first kind with integer degree and order (cf. [, (4.7.8] with m. P m n (cos θp m n (cos θ sin θ dθ = (n + m! m (n m! δ m,m, (7 Theorem 9. Let m N, n N with m n, θ [, π], φ, φ [, π. Then P n (cos γp m n (cos θ sin θ dθ = m Pm n (cos θ cos(m(φ φ, cos γ = cos θ cos θ + sin θ sin θ cos(φ φ.

8 8 H.S. Cohl and H. Volkmer Proof. We start with the addition theorem for spherical harmonics (cf. [, (4.8.], namely P n (cos γ = n m= n (n m! (n + m! Pm n (cos θp m n (cos θ e im(φ φ, (8 P n : C C, for n N, is the Legendre polynomial which can be defined in terms of the terminating Gauss hypergeometric series (see for instance [, Chapters 5, 8] as follows P n (z := F ( n, n + ; z. We then take advantage of the order orthogonality relation for the Ferrers functions of the first kind with integer degree and order. If we multiply both sides of (8 by (sin θ P m n (cos θ and integrate over θ (, π, by using (7 we obtain the desired result. Theorem 9, originating from (8, is the only example of a definite integral that we could find using the order orthogonality relation for the Ferrers functions of the first kind (7. Therefore we highly suspect that this result is previously known, and include it mainly for completeness sake. It would however be very interesting to find another example using this orthogonality relation. 3.3 Orthogonality for Chebyshev polynomials of the f irst kind Here we take advantage of orthogonality from Chebyshev polynomials of the first kind (cf. [, 8.3] T m (cos θt n (cos θdθ = π ɛ n δ m,n, (9 T n : C C, for n N, is the Chebyshev polynomial of the first kind which can be defined in terms of the terminating Gauss hypergeometric series (see [, Chapter 8] ( n, n T n (z = F ; z. The Chebyshev polynomials of the first kind satisfy the identity [, (8.5.] T n (cos θ = cos(nθ. Theorem. Let m, n Z, σ, σ (,. Then Q m / (χ cos(nψdψ = π( m sinh σ sinh σ Γ ( n m + Γ ( n + m + Pn / m (cosh σ <Q m n / (cosh σ >, (3 χ = coth σ coth σ csch σ csch σ cos ψ. (3 Proof. We start with toroidal coordinates on R 3, namely x = a sinh σ cos φ cosh σ cos ψ, y = a sinh σ sin φ cosh σ cos ψ, z = a sin ψ cosh σ cos ψ,

9 Definite Integrals using Orthogonality and Integral Transforms 9 a >, σ (,, ψ, φ [, π. The reciprocal distance between two points x, x R 3 is given algebraically by x x = (cosh σ cos ψ(cosh σ cos ψ a sinh σ sinh σ [ cosh σ cosh σ cos(ψ ψ ] / sinh σ sinh σ cos(φ φ, (σ, ψ, φ are the toroidal coordinates corresponding to the point x. reciprocal square root identity (see for instance [5, (3.] Using Heine s = z x π ɛ m Q m / (zt m (x, m= z > and x [, ], we can obtain a Fourier cosine series representation for the reciprocal distance between two points in toroidal coordinates on R 3, namely x x = (cosh σ cos ψ(cosh σ cos ψ πa sinh σ sinh σ ɛ m cos(m(φ φ Q m / (χ, χ > is given by (3. We can further expand the associated Legendre function of the second kind using the following addition theorem (cf. [7, (8.795.] Q m / (χ = ( m sinh σ sinh σ n= m= ɛ n cos(n(ψ ψ Γ ( n m + Γ ( n + m + Pn / m (cosh σ <Q m n / (cosh σ >. (3 Note that with the above addition theorem, we have the expansion of the reciprocal distance between two points in terms of the separable harmonics in toroidal coordinates x x = (cosh σ cos ψ(cosh σ πa cos ψ ( m ɛ m cos(m(φ φ m= ɛ n cos(n(ψ ψ Γ ( n m + Γ ( n + m + Pn / m (cosh σ <Q m n / (cosh σ > n= (see also [6, 6.] and [4]. If we relabel ψ ψ ψ and multiply both sides of (3 by cos(nψ and integrate over ψ [, π], then by (9 we have obtained (3. The interchange of infinite sum and integral is justified by (, (. Acknowledgements This work was conducted while H.S. Cohl was a National Research Council Research Postdoctoral Associate in the Information Technology Laboratory at the National Institute of Standards and Technology, Gaithersburg, Maryland, USA. The authors would also like to acknowledge two anonymous referees whose comments helped improve this paper.

10 H.S. Cohl and H. Volkmer References [] Abramowitz M., Stegun I.A., Handbook of mathematical functions with formulas, graphs, and mathematical tables, National Bureau of Standards Applied Mathematics Series, Vol. 55, U.S. Government Printing Office, Washington, D.C., 964. [] Askey R., Orthogonal polynomials and special functions, Society for Industrial and Applied Mathematics, Philadelphia, Pa., 975. [3] Cohl H.S., Fourier, Gegenbauer and Jacobi expansions for a power-law fundamental solution of the polyharmonic equation and polyspherical addition theorems, arxiv: [4] Cohl H.S., Erratum: Developments in determining the gravitational potential using toroidal functions, Astronom. Nachr. 333 (, [5] Cohl H.S., Dominici D.E., Generalized Heine s identity for complex Fourier series of binomials, Proc. R. Soc. Lond. Ser. A 467 (, , arxiv:9.6. [6] Cohl H.S., Tohline J.E., Rau A.R.P., Srivastava H.M., Developments in determining the gravitational potential using toroidal functions, Astronom. Nachr. 3 (, [7] Gradshteyn I.S., Ryzhik I.M., Table of integrals, series, and products, seventh ed., Elsevier/Academic Press, Amsterdam, 7. [8] Hardy G.H., Further researches in the theory of divergent series and integrals, Trans. Cambridge Philos. Soc. (98, 48. [9] MacRobert T.M., Spherical harmonics. An elementary treatise on harmonic functions with applications, nd ed., Methuen & Co. Ltd., London, 947. [] Magnus W., Oberhettinger F., Soni R.P., Formulas and theorems for the special functions of mathematical physics, 3rd ed., Die Grundlehren der mathematischen Wissenschaften, Bd. 5, Springer-Verlag, New York, 966. [] Morse P.M., Feshbach H., Methods of theoretical physics, Vols.,, McGraw-Hill Book Co. Inc., New York, 953. [] Olver F.W.J., Lozier D.W., Boisvert R.F., Clark C.W. (Editors, NIST handbook of mathematical functions, Cambridge University Press, Cambridge,. [3] Prudnikov A.P., Brychkov Y.A., Marichev O.I., Integrals and series. Vol. 3. More special functions, Gordon and Breach Science Publishers, New York, 99. [4] Schäfke F.W., Einführung in die Theorie der speziellen Funktionen der mathematischen Physik, Die Grundlehren der mathematischen Wissenschaften, Bd. 8, Springer-Verlag, Berlin, 963. [5] Watson G.N., A treatise on the theory of Bessel functions, nd ed., Cambridge Mathematical Library, Cambridge University Press, Cambridge, 944.