Abstract

A method is developed for solutions of two sets of triple integral equations involving associated Legendre functions of imaginary arguments. The solution of each set of triple integral equations involving associated Legendre functions is reduced to a Fredholm integral equation of the second kind which can be solved numerically.

1. Introduction

Dual integral equations involving Legendre functions have been solved by Babloian [1]. He applied these equations to problems of potential theory and to a torsion problem. Later on Pathak [2] and Mandal [3] who considered dual integral equations involving generalized Legendre functions which have more general solution than the ones considered by Babloian [1]. Recently, Singh et al. [4] considered dual integral equations involving generalized Legendre functions, and their results are more general than those in [13].

In the analysis of mixed boundary value problems, we often encounter triple integral equations. Triple integral equations involving Legendre functions have been studied by Srivastava [5]. Triple integral equations involving Bessel functions have also been considered by Cooke [69], Tranter [10], Love and Clements [11], Srivastava [12], and most of these authors reduced the solution into a solution of Fredholm integral equation of the second kind. The relevant references for dual and triple integral equations are given in the book of Sneddon [13].

In this paper, a method is developed for solutions of two sets of triple integral equations involving generalized Legendre functions in Sections 3 and 4. Each set of triple integral equations is reduced to a Fredholm integral equation of the second kind which may be solved numerically. The aim of this paper is to find a more general solution for the type of integral equations given in [15] and to develop an easier method for solving triple integral equations in general.

2. Integral Involved Generalized Legendre Functions and Some Useful Results

We first summarize some known results needed in the paper.

We find from [14, equation (21), page 330] that2𝜋Γ12𝜇sinh(𝛼𝑐)𝜇0𝑃𝜇1/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)cos(𝜏𝑥)𝑑𝜏=𝑐cosh(𝛼𝑐)cosh(𝑥𝑐)𝜇1/2𝐻(𝛼𝑥),(2.1)where 𝜇<1/2 and from [4], we obtain2𝜋3/2Γ12+𝜇sinh(𝛼𝑐)𝜇×0Γ12𝜏𝜇+𝑖𝑐Γ12𝜏𝜇𝑖𝑐sinh(𝑓𝜏)sin(𝑥𝜏)𝑃𝜇1/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=𝑐cosh(𝑥𝑐)cosh(𝛼𝑐)𝜇1/2𝐻(𝑥𝛼),(2.2)where 𝜇>1/2 and 𝐻() denotes the Heaviside unit function. Furthermore, 𝑐=𝜋/𝑓, 𝑓>0 and 𝑃𝜇1/2+𝑖(𝜏/𝑐)(cosh𝛼𝑐) is the generalized Legendre function defined in [15, page 370]. From [4, 16], the generalized Mehler-Fock transform is defined by𝜓=cosh(𝛼𝑐)0𝑃𝜇1/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝐹(𝜏)𝑑𝜏,(2.3)and its inversion formula is𝐹(𝜏)=𝑓𝜏𝜋21sinh(𝑓𝜏)Γ2𝜏𝜇+𝑖𝑐Γ12𝜏𝜇𝑖𝑐×0𝑃𝜇1/2+𝑖(𝜏/𝑐)𝜓cosh(𝛼𝑐)cosh(𝛼𝑐)sinh(𝛼𝑐)𝑑𝛼.(2.4)Equations (2.1) and (2.2) are of form (2.3). From the inversion formula given by (2.4), (2.1), and (2.2), it follows thatcos(𝑥𝜏)𝜏=Γsinh(𝑓𝜏)Γ1/2𝜇+𝑖(𝜏/𝑐)1/2𝜇𝑖(𝜏/𝑐)×2𝜋Γ1/2𝜇𝑥sinh(𝛼𝑐)1+𝜇𝑃𝜇1/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝛼cosh(𝛼𝑐)cosh(𝑥𝑐)𝜇+1/21,𝜇<2,(2.5)sin(𝑥𝜏)𝜏=𝜋21Γ1/2+𝜇𝑥0sinh(𝛼𝑐)1𝜇𝑃𝜇1/2+𝑖(𝜏/𝑐)(cosh𝛼𝑐)𝑑𝛼cosh(𝑥𝑐)cosh(𝛼𝑐)1/2𝜇1,𝜇>2.(2.6)

The inversion theorem for Fourier cosine transforms and the results (2.1) and (2.2) lead to𝑃𝜇1/2+𝑖(𝜏/𝑐)=cosh(𝛼𝑐)2𝜋𝑐sinh𝜇(𝛼𝑐)Γ1/2𝜇𝛼0cos(𝜏𝑠)𝑑𝑠cosh(𝛼𝑐)cosh(𝑠𝑐)𝜇+1/21,𝜇<2,𝑃(2.7)𝜇1/2+𝑖(𝜏/𝑐)=cosh(𝛼𝑐)2𝜋𝑐sinh𝜇Γ×(𝛼𝑐)]Γ1/2+𝜇sin(𝑓𝜏)Γ1/2𝜇+𝑖(𝜏/𝑐)1/2𝜇𝑖(𝜏/𝑐)𝛼sin(𝜏𝑠)𝑑𝑠cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇1,𝜇>2.(2.8)

If (𝑡) is monotonically increasing and differentiable for 𝑎<𝑡<𝑏 and (𝑡)0 in this interval, then the solutions of the equations𝑡𝑎𝑓(𝑥)𝑑𝑥(𝑡)(𝑥)𝛼=𝑔(𝑡),𝑎<𝑡<𝑏,0<𝛼<1,(2.9)𝑏𝑡𝑓(𝑥)𝑑𝑥(𝑥)(𝑡)𝛼=𝑔(𝑡),𝑎<𝑡<𝑏,0<𝛼<1,(2.10)are given by Sneddon [13] as𝑓(𝑥)=sin(𝜋𝛼)𝜋𝑑𝑑𝑥𝑥𝑎(𝑡)𝑔(𝑡)𝑑𝑡(𝑥)(𝑡)1𝛼,𝑎<𝑥<𝑏,(2.11)𝑓(𝑥)=sin(𝜋𝛼)𝜋𝑑𝑑𝑥𝑏𝑥(𝑡)𝑔(𝑡)𝑑𝑡(𝑡)(𝑥)1𝛼,𝑎<𝑥<𝑏,(2.12)respectively, where the prime denotes the derivative with respect to 𝑡.

3. Triple Integral Equations with Generalized Legendre Functions: Set I

In this section, we will find solution of the following triple integral equations:01𝜏𝐴(𝜏)sinh(𝜏𝑓)Γ2𝜇1𝜏+𝑖𝑐Γ12𝜇1𝜏𝑖𝑐𝑃𝜇11/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=0,0<𝛼<𝑎,(3.1)0𝐴(𝜏)𝑃𝜇21/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=𝑓(𝛼),𝑎<𝛼<𝑏,(3.2)01𝜏𝐴(𝜏)sinh(𝜏𝑓)Γ2𝜇3𝜏+𝑖𝑐Γ12𝜇3𝜏𝑖𝑐𝑃𝜇31/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=0,𝑏<𝛼<,(3.3) where 𝐴(𝜏) is an unknown function to be determined, 𝑓(𝛼) is a known function, and 𝑃𝜇1/2+𝑖(𝜏/𝑐)[cosh(𝛼𝑐)] is the generalized Legendre function defined in Section 2 and 1/2<𝜇1<1/2, 1/2<𝜇2<1/2, 𝜇3>1/2.

The trial solution of (3.1), (3.2), and (3.3) can be written as𝐴(𝜏)=𝑏0𝜓(𝑡)cos(𝜏𝑡)𝑑𝑡,(3.4)where 𝜓(𝑡) is an unknown function to be determined. On integrating (3.4) by parts, we get𝐴(𝜏)=𝜓(𝑏)sin(𝜏𝑏)𝜏1𝜏𝑏0𝜓(𝑡)sin(𝜏𝑡)𝑑𝑡,(3.5)where the prime denotes the derivative with respect to 𝑡.

Substituting (3.5) into (3.3), interchanging the order of integrations and using (2.2), we find that (3.3) is satisfied identically. Substituting (3.5) into (3.1) and using the integral defined by (2.2), we obtain𝜓(𝑏)cosh(𝑏𝑐)cosh(𝛼𝑐)1/2𝜇1𝑏𝛼𝜓(𝑡)𝑑𝑡cosh(𝑡𝑐)cosh(𝛼𝑐)1/2𝜇1=0,0<𝛼<𝑎.(3.6)Equation (3.6) is equivalent to the following integral equation:𝑑𝑑𝛼𝑏𝛼𝑐sinh(𝑡𝑐)𝜓(𝑡)𝑑𝑡cosh(𝑡𝑐)cosh(𝛼𝑐)1/2𝜇1=0,0<𝛼<𝑎.(3.7)By substituting (3.4) into (3.2), interchanging the order of integrations and using the integral defined by (2.1) we find that𝑐𝛼0𝜓(𝑡)𝑑𝑡cosh(𝛼𝑐)cosh(𝑡𝑐)1/2+𝜇2=2𝜋Γ12𝜇2sinh(𝛼𝑐)𝜇2𝑓(𝛼),𝑎<𝛼<𝑏,𝜇2<12.(3.8)For obtaining the solution of the problem, we need to solve two Abel's type integral equations (3.7) and (3.8).

We assume that𝑑𝑑𝛼𝑏𝛼𝑐sinh(𝑡𝑐)𝜓(𝑡)𝑑𝑡cosh(𝑡𝑐)cosh(𝛼𝑐)1/2𝜇1=𝜙(𝛼),𝑎<𝛼<𝑏.(3.9)The above equation is of the same form as (3.7) and defined in a different region. Equation (3.9) is of form (2.12). Hence, the solution of the integral equation (3.9) can be written as𝜓(𝑡)=cos𝜋𝜇1𝜋𝑏𝑡𝜙(𝛼)𝑑𝛼cosh(𝑐𝛼)cosh(𝑡𝑐)1/2+𝜇11,2<𝜇1<12,𝑎<𝑡<𝑏.(3.10)

The solution of Abel's type integral equations (2.11) together with (3.7) and (3.9) leads to𝜓(𝑡)=cos𝜋𝜇1𝜋𝑏𝑎𝜙(𝛼)𝑑𝛼cosh(𝑐𝛼)cosh(𝑡𝑐)1/2+𝜇11,2<𝜇1<12,0<𝑡<𝑎.(3.11)

Equations (3.10) and (3.11) mean that (3.7) is satisfied identically. Equation (3.8) can be rewritten in the form𝑎0𝜓(𝑡)𝑑𝑡cosh(𝛼𝑐)cosh(𝑡𝑐)1/2+𝜇2+𝛼𝑎𝜓(𝑡)𝑑𝑡cosh(𝛼𝑐)cosh(𝑡𝑐)1/2+𝜇2=1𝑐2𝜋Γ1𝜇2𝑓(𝛼)sinh(𝛼𝑐)𝜇2,𝑎<𝛼<𝑏.𝑓(3.12)Substituting the expression for 𝜓(𝑡) from (3.11) and (3.10) into the first and second integral of (3.12) we obtain𝛼𝑎𝑆(𝑡)𝑑𝑡cosh(𝛼𝑐)cosh(𝑡𝑐)1/2+𝜇2=𝐹(𝛼)𝑎0𝑑𝑡cosh(𝛼𝑐)cosh(𝑡𝑐)1/2+𝜇2𝑏𝑎𝜙(𝑢)𝑑𝑡cosh(𝑐𝑢)cosh(𝑡𝑐)1/2+𝜇2,𝑎<𝑡<𝑏,(3.13)where𝑆(𝑡)=𝑏𝑡𝜙(𝑢)𝑑𝑡cosh(𝑐𝑢)cosh(𝑡𝑐)1/2+𝜇1,(3.14)𝐹(𝛼)=2𝜋Γ1𝜇2𝑓(𝛼)sinh(𝛼𝑐)𝜇2𝜇𝑐cos1𝜋.(3.15)

Assuming that the right-hand side of (3.13) is a known function of 𝛼 it has the form of (2.9), whose solution is given by𝑆(𝑡)=cos𝜋𝜇2𝜋𝑑𝑑𝑡𝑡𝑎𝑐sinh(𝑐𝛼)𝐹(𝛼)𝑑𝛼cosh(𝑐𝑡)cosh(𝑐𝛼)1/2𝜇21𝐼(𝑡),𝑎<𝑡<𝑏,2<𝜇2<12,(3.16)where𝐼(𝑡)=cos𝜋𝜇2𝜋𝑑𝑑𝑡𝑡𝑎𝑐sinh(𝑐𝛼)𝑑𝛼cosh(𝑐𝑡)cosh(𝑐𝛼)1/2𝜇2𝑎0𝑑𝑝cosh(𝑐𝛼)cosh(𝑐𝑝)1/2+𝜇2×𝑏𝑎𝜙(𝑢)𝑑𝑢cosh(𝑐𝑢)cosh(𝑐𝑝)1/2+𝜇21,𝑎<𝑡<𝑏,2<𝜇2<12.(3.17)From the integral𝑑𝑑𝑡𝑡𝑎𝑐sinh(𝑐𝛼)𝑑𝛼cosh(𝑐𝑡)cosh(𝑐𝛼)1/2𝜇2cosh(𝑐𝛼)cosh(𝑐𝑝)1/2+𝜇2=𝑐sinh(𝑐𝑡)cosh(𝑐𝑡)cosh(𝑐𝑝)cosh(𝑐𝑎)cosh(𝑐𝑝)1/2𝜇2cosh(𝑐𝑡)cosh(𝑐𝑎)1/2𝜇21,𝑝<𝑎<𝑡,2<𝜇2<12,(3.18)we then obtain𝜇𝐼(𝑡)=𝑐cos2𝜋sinh(𝑐𝑡)𝜋cosh(𝑐𝑡)cosh(𝑐𝑎)1/2𝜇2𝑎0cosh(𝑐𝑎)cosh(𝑐𝑝)1/2𝜇2𝑑𝑝×cosh(𝑐𝑡)cosh(𝑐𝑝)𝑏𝑎𝜙(𝑢)𝑑𝑢cosh(𝑐𝑢)cosh(𝑐𝑝)1/2+𝜇1.(3.19)Equation (3.14) is an Abel-type equation. Hence, its solution is𝜇𝜙(𝑢)=cos1𝜋𝜋𝑑𝑑𝑢𝑏𝑢𝑐sinh(𝑐𝑣)𝑆(𝑣)𝑑𝑣cosh(𝑣𝑐)cosh(𝑢𝑐)1/2𝜇11,𝑎<𝑢<𝑏,2<𝜇1<12,(3.20)𝑅(𝑝)=𝑏𝑎𝜙(𝑢)𝑑𝑢cosh(𝑐𝑢)cosh(𝑐𝑝)1/2+𝜇1.(3.21)Substituting the expression for 𝜙(𝑢) from (3.20) into (3.21), integrating by parts, and finally interchanging the order of integrations in second integral, we arrive at𝜇𝑅(𝑝)=𝑐cos1𝜋𝜋1cosh(𝑐𝑎)cosh(𝑐𝑝)1/2+𝜇1𝑏𝑎𝑆(𝑣)sinh(𝑐𝑣)𝑑𝑣cosh(𝑐𝑣)cosh(𝑐𝑎)1/2𝜇112+𝜇1𝑏𝑎×𝑆(𝑣)sinh(𝑐𝑣)𝑑𝑣𝑣𝑎𝑐sinh(𝑐𝑢)𝑑𝑢cosh(𝑐𝑢)cosh(𝑐𝑝)3/2+𝜇1cosh(𝑐𝑣)cosh(𝑐𝑢)1/2𝜇1.(3.22)The integral 𝑣𝑎𝑐sinh(𝑐𝑢)𝑑𝑢cosh(𝑐𝑢)cosh(𝑐𝑝)3/2+𝜇1cosh(𝑐𝑣)cosh(𝑐𝑢)1/2𝜇1=[cosh(𝑐𝑣)cosh(𝑐𝑎)]1/2+𝜇1(𝜇11+1/2)[cosh(𝑐𝑣)cosh(𝑐𝑝)][cosh(𝑐𝑎)cosh(𝑐𝑝)]𝑝<𝑎<𝑣,2<𝜇1<12(3.23)together with (3.22) leads to𝑅(𝑝)=𝑐cos𝜋𝜇1𝜋cosh(𝑎𝑐)cosh(𝑝𝑐)1/2𝜇1×𝑏𝑎𝑆(𝜈)sinh(𝑐𝜈)𝑑𝜈cosh(𝜈𝑐)cosh(𝑝𝑐)cosh(𝜈𝑐)cosh(𝑎𝑐)1/2𝜇1.(3.24)From (3.19), (3.21), and (3.24), we obtain𝐼(𝑡)=𝑏𝑎𝑆(𝜈)𝐾(𝜈,𝑡)𝑑𝜈,(3.25)where𝑐𝐾(𝜈,𝑡)=2cos𝜋𝜇1cos𝜋𝜇2sinh(𝑐𝑡)sinh(𝑐𝜈)𝜋2cosh(𝑐𝑡)cosh(𝑐𝑎)1/2𝜇2cosh(𝑐𝜈)cosh(𝑐𝑎)1/2𝜇1×𝑎0cosh(𝑐𝑎)cosh(𝑐𝑝)1𝜇1𝜇2𝑑𝑝.cosh(𝑐𝑡)cosh(𝑐𝑝)cosh(𝑐𝜈)cosh(𝑐𝑝)(3.26)From (3.25), (3.16) can be written as𝑆(𝑡)+𝑏𝑎𝑆(𝜈)𝐾(𝜈,𝑡)𝑑𝜈=cos𝜋𝜇2𝜋𝑑𝑑𝑡𝑡𝑎𝑐sinh(𝑐𝛼)𝐹(𝛼)𝑑𝛼cosh(𝑐𝑡)cosh(𝑐𝛼)1/2𝜇2,𝑎<𝑡<𝑏.(3.27)Equation (3.27) is a Fredholm integral equation of the second kind with kernel 𝐾(𝜈,𝑡). The kernel is defined by (3.26). The integral in (3.26) cannot be solved analytically, but for particular values of 𝜇1 and 𝜇2 the values of 𝐾(𝜈,𝑡) can be found numerically. Hence, the numerical solution of Fredholm integral equation (3.27) can be obtained for particular value of 𝑓(𝛼), 𝜇1, and 𝜇2 to find numerical values of 𝑆(𝑡). Making use of (3.20), (3.11), and (3.10), the numerical results for 𝜓(𝑡) can be obtained. Finally, making use of (3.4) the numerical results for 𝐴(𝜏) can be obtained.

4. Triple Integral Equations with Generalized Legendre Functions: Set II

In this section, we will find the solution of the following triple integral equations:0𝜏𝐴(𝜏)𝑃𝜇11/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=0,0<𝛼<𝑎,(4.1)01sinh(𝜏𝑓)Γ2𝜇2𝜏+𝑖𝑐Γ12𝜇2𝜏𝑖𝑐𝐴(𝜏)𝑃𝜇21/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=𝑓(𝛼),𝑎<𝛼<𝑏,(4.2)0𝜏𝐴(𝜏)𝑃𝜇31/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=0,𝑏<𝛼,(4.3) where 𝜇1>1/2, 1/2<𝜇2<1/2, 1/2<𝜇3<1/2.

We assume that0𝜏𝐴(𝜏)𝑃𝜇31/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏=𝑀(𝛼),0<𝛼<𝑏.(4.4)The inversion formula for generalized Mehler-Fock transforms (2.4) together with (4.3) and (4.4) implies that𝑓𝐴(𝜏)=𝜋21sinh(𝑓𝜏)Γ2𝜇3𝜏+𝑖𝑐Γ12𝜇3𝜏𝑖𝑐×𝑏0sinh(𝑢𝑐)𝑃𝜇31/2+𝑖(𝜏/𝑐)cosh(𝑢𝑐)𝑀(𝑢)𝑑𝑢.(4.5)

Multiplying (4.1) by [sinh(𝛼𝑐)]1𝜇1/[cosh(𝑥𝑐)cosh(𝛼𝑐)]1/2𝜇1, integrating both sides from 0 to 𝑥 and with respect to 𝛼, and then using (2.6) we obtain0𝐴(𝜏)sin(𝑥𝜏)𝑑𝜏=0,0<𝑥<𝑎.(4.6)Substituting the value of 𝐴(𝜏) from (4.5) into (4.6), interchanging the order of integrations, and using the integral (2.2), we get𝑥0sinh(𝑢𝑐)𝑀(𝑢)𝑑𝑢cosh(𝑥𝑐)cosh(𝑢𝑐)1/2𝜇3=0,𝜇31>2,0<𝑥<𝑎.(4.7)Substituting the value of 𝐴(𝜏) from (4.5) into (4.2) and interchanging the order of integrations we arrive at𝑏0sinh(𝑢𝑐)𝑀(𝑢)𝐾2(𝑢,𝛼)𝑑𝑢=𝑓(𝛼),𝑎<𝛼<𝑏,(4.8)where𝐾2(𝑢,𝛼)=0𝑓𝜋2Γ12𝜇2𝜏𝑖𝑐Γ12𝜇2𝜏+𝑖𝑐Γ12𝜇3𝜏+𝑖𝑐Γ12𝜇3𝜏𝑖𝑐×sinh2(𝑓𝜏)𝑃𝜇31/2+𝑖(𝜏/𝑐)𝑃cosh(𝑢𝑐)𝜇21/2+𝑖(𝜏/𝑐)cosh(𝛼𝑐)𝑑𝜏,(4.9)and then (2.8) and (2.2) imply that𝐾2(𝑢,𝛼)=𝑐𝜋Γ1/2+𝜇2Γ1/2+𝜇3sinh(𝛼𝑐)𝜇2sinh(𝑢𝑐)𝜇3×max(𝛼,𝑢)𝑑𝑠cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇2cosh(𝑠𝑐)cosh(𝑢𝑐)1/2𝜇3,𝜇31>2,𝜇21>2.(4.10)

Equation (4.7) is an Abel-type equation and has the form (2.9). Hence, the solution of (4.7) is𝑀(𝑢)=0,0<𝑢<𝑎.(4.11)

Using (4.10) and (2.5), (4.8) can be written in the form𝑏𝑎sinh(𝑢𝑐)1𝜇3𝑀(𝑢)𝑑𝑢max(𝛼,𝑢)𝑑𝑠cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇2cosh(𝑠𝑐)cosh(𝑢𝑐)1/2𝜇3=Γ1/2+𝜇2Γ1/2+𝜇2sinh(𝛼𝑐)𝜇3𝑐𝜋𝑓(𝛼)=𝐹1(𝛼),say,𝑎<𝛼<𝑏.(4.12)Using the formula𝑏𝑎𝑑𝑢max(𝛼,𝑢)𝑑𝑠=𝑏𝛼𝑑𝑠𝑠𝑎𝑑𝑢+𝑏𝑑𝑠𝑏𝑎𝑑𝑢,(4.13)we can write (4.12) in the form𝑏𝛼𝑆1(𝑠)𝑑𝑠cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇2=𝐹1(𝛼)𝑏𝑑𝑠cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇2×𝑏𝑎𝑀(𝑢)sinh(𝑢𝑐)1𝜇3𝑑𝑢cosh(𝑠𝑐)cosh(𝑢𝑐)1/2𝜇2,𝑎<𝛼<𝑏,(4.14)where𝑆1(𝑠)=𝑠𝑎𝑀(𝑢)sinh(𝑢𝑐)1𝜇3𝑑𝑢cosh(𝑠𝑐)cosh(𝛼𝑐)1/2𝜇2,𝑎<𝑠<𝑏.(4.15)

Assuming that the right-hand side of (4.14) is known function equation and (4.14) has the form of (2.10), hence the solution of (4.14) can be written as𝑆1𝑐(𝑠)=𝜋cos𝜋𝜇2𝑑𝑑𝑠𝑏𝑠𝐹1(𝛼)sinh(𝛼𝑐)𝑑𝛼cosh(𝛼𝑐)cosh(𝑠𝑐)1/2+𝜇2+𝐼11(𝑠),𝑎<𝑠<𝑏,2<𝜇2<12,(4.16)where𝐼1𝑐(𝑠)=𝜋cos𝜋𝜇2𝑑𝑑𝑠𝑏𝑠sinh(𝛼𝑐)𝑑𝛼cosh(𝛼𝑐)cosh(𝑠𝑐)1/2+𝜇2×𝑏𝑑𝑝cosh(𝑝𝑐)cosh(𝛼𝑐)1/2𝜇2𝑏𝑎𝑀(𝑢)sinh(𝑐𝑢)1𝜇3𝑑𝛼cosh(𝑝𝑐)cosh(𝑢𝑐)1/2𝜇3,𝑎<𝑠<𝑏.(4.17)

Equation (4.17) is simplified to𝐼1(𝑠)=𝑐cos𝜋𝜇2𝜋sinh(𝑠𝑐)cosh(𝑏𝑐)cosh(𝑠𝑐)1/2+𝜇2𝑏cosh(𝑐𝑝)cosh(𝑏𝑐)1/2+𝜇2𝑑𝑝×cosh(𝑠𝑐)cosh(𝑐𝑝)𝑏𝑎𝑀(𝑢)sinh(𝑐𝑢)1𝜇3𝑑𝑢cosh(𝑐𝑝)cosh(𝑐𝑢)1/2𝜇3,𝑎<𝑠<𝑏.(4.18)Let𝑅1(𝑝)=𝑏𝑎𝑀(𝑢)sinh(𝑐𝑢)1𝜇3𝑑𝑢cosh(𝑝𝑐)cosh(𝑢𝑐)1/2𝜇3.(4.19)Equation (4.15) is of the form of (2.9). Hence, its solution is𝑀(𝑢)sinh(𝑐𝑢)1𝜇3=𝑐cos𝜋𝜇3𝜋𝑑𝑑𝑢𝑢𝑎𝑆1(𝑠)sinh(𝑠𝑐)𝑑𝑠cosh(𝑢𝑐)cosh(𝑠𝑐)1/2+𝜇3,𝑎<𝑢<𝑏.(4.20)

Substituting the expression for 𝑀(𝑢) from (4.20) into (4.19) and integrating by parts and then using the following integral:𝑏𝑠𝑐sinh(𝑢𝑐)𝑑𝑢cosh(𝑝𝑐)cosh(𝑢𝑐)3/2𝜇3cosh(𝑢𝑐)cosh(𝑠𝑐)1/2+𝜇3=cosh(𝑏𝑐)cosh(𝑐𝑠)1/2𝜇31/2𝜇3cosh(𝑐𝑠)cosh(𝑐𝑝)cosh(𝑐𝑝)cosh(𝑐𝑏)1/2𝜇3,1𝑠<𝑏<𝑝,2<𝜇2<12,(4.21)we find that𝑅1𝜇(𝑝)=𝑐cos3𝜋𝜋cosh(𝑐𝑝)cosh(𝑏𝑐)1/2+𝜇3×𝑏𝑎𝑆1(𝑢)sinh(𝑐𝑢)𝑑𝑢cosh(𝑐𝑝)cosh(𝑢𝑐)cosh(𝑏𝑐)cosh(𝑐𝑢)1/2+𝜇3.(4.22)Making use of (4.18), (4.19), and (4.22), we find that𝐼1(𝑠)=𝑏𝑎𝑆1(𝑢)𝐾2(𝑢,𝑠)𝑑𝑢,(4.23)where𝐾2𝑐(𝑢,𝑠)=2cos𝜋𝜇2cos𝜋𝜇3sinh(𝑠𝑐)sinh(𝑢𝑐)𝜋2cosh(𝑏𝑐)cosh(𝑠𝑐)1/2+𝜇2cosh(𝑏𝑐)cosh(𝑐𝑢)1/2+𝜇3×𝑏cosh(𝑐𝑝)cosh(𝑏𝑐)1+𝜇2+𝜇3𝑑𝑝.cosh(𝑠𝑐)cosh(𝑐𝑝)cosh(𝑐𝑝)cosh(𝑐𝑢)(4.24)Using (4.17) and (4.23), (4.16) can be written in the form𝑆1(𝑠)+𝑏𝑎𝑆1(𝑢)𝐾2(𝑢,𝑠)𝑑𝑢=𝑐𝜋cos𝜋𝜇2𝑑𝑑𝑠𝑏𝑠𝐹1(𝛼)sinh(𝛼𝑐)𝑑𝛼cosh(𝛼𝑐)cosh(𝑠𝑐)1/2+𝜇2,𝑎<𝑠<𝑏.(4.25)

Equation (4.25) is a Fredholm integral equation of the second kind with kernel defined by (4.24). The Fredholm integral equation (4.25) may be solved to find numerical values of 𝑆1(𝑠) for particular values of 𝑓(𝛼). And hence from (4.20) and (4.5), the numerical values for 𝐴(𝜏) can be obtained for particular values of 𝑓(𝛼), 𝜇2, and 𝜇3.

5. Conclusions

The solution of the two sets of triple integral equations involving generalized Legendre functions is reduced to the solution of Fredholm integral equations of the second kind which can be solved numerically.