Hindawi Publishing Corporation International Journal of Analysis Volume 2013, Article ID 652541, 12 pages http://dx.doi.org/10.1155/2013/652541
Research Article Generalized Abel Inversion Using Extended Hat Functions Operational Matrix Manoj P. Tripathi,1,2 Ram K. Pandey,1 Vipul K. Baranwal,1 and Om P. Singh1 1 2
Department of Applied Mathematics, Indian Institute of Technology, Banaras Hindu University, Varanasi 221005, India Department of Mathematics, Udai Pratap Autonomous College, Varanasi 221002, India
Correspondence should be addressed to Om P. Singh;
[email protected] Received 24 September 2012; Accepted 2 April 2013 Academic Editor: FrΒ΄edΒ΄eric Robert Copyright Β© 2013 Manoj P. Tripathi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abel type integral equations play a vital role in the study of compressible flows around axially symmetric bodies. The relationship between emissivity and the measured intensity, as measured from the outside cylindrically symmetric, optically thin extended radiation source, is given by this equation as well. The aim of the present paper is to propose a stable algorithm for the numerical π¦
π½
inversion of the following generalized Abel integral equation: πΌ(π¦) = π(π¦) β«πΌ ((ππβ1 π(π))/(π¦π β ππ )πΎ )ππ + π(π¦) β«π¦ ((ππβ1 π(π))/(ππ β π¦π )πΎ )ππ, πΌ β€ π¦ β€ π½, 0 < πΎ < 1, using our newly constructed extended hat functions operational matrix of integration, and give an error analysis of the algorithm. The earlier numerical inversions available for the above equation assumed either π(π¦) = 0 or π(π¦) = 0.
1. Introduction
The analytical inversion formula for (1) is given as [4]
Abel integral equation [1] occurs in many branches of science and technology, such as plasma diagnostics and flame studies, where the most common problem of deduction of radial distributions of some important physical quantity from measurement of line-of-sight projected values is encountered. For a cylindrically symmetric, optically thin plasma source, the relation between radial distribution of the emission coefficient and the intensity measured from outside of the radial source is described by Abel transform. The challenging task of reconstruction of emission coefficient from its projection is known as Abel inversion. The earliest application, due to Mach [2], arose in the study of compressible flows around axially symmetric bodies. The Abel integral equation is given by 1
π (π) π
π¦
βπ2 β π¦2
πΌ (π¦) = 2 β«
ππ,
0 β€ π¦ β€ 1,
(1)
where π(π) and πΌ(π¦) represent, respectively, the emissivity and measured intensity, as measured from outside the source [3].
ππΌ (π¦) 1 1 1 ππ¦, π (π) = β β« π π βπ¦2 β π2 π (π¦)
0 β€ π β€ 1.
(2)
There are several analytic and numerical inversion formulae available in the literature [1, 5β20]. Singh et al. [19] constructed an operational matrix of integration based on orthonormal Bernstein polynomials and used it to propose a stable algorithm to invert the following form of Abel integral equation: π¦
πΌ (π¦) = 2 β«
0
π (π) π βπ2 β π¦2
ππ,
0 β€ π¦ β€ 1.
(3)
In 2010, Singh et al. [20] constructed yet another operational matrix of integration based on orthonormal Bernstein polynomials and used it to propose an algorithm to invert the Abel integral equation (1).
2
International Journal of Analysis π‘ β (π β 1) β { , { { β { { { { { (π + 1) β β π‘ ππ (π‘) = { , { β { { { { { { {0,
In 2008, Chakrabarti [21] employed a direct function theoretic method to determine the closed form solution of the following generalized Abel integral equation: πΌ (π¦) = π (π¦) β«
π¦
πΌ
π½ ππβ1 π (π) ππβ1 π (π) ππ + π (π¦) β« πΎ πΎ ππ, (π¦π β ππ ) π¦ (ππ β π¦π )
πΌ β€ π¦ β€ π½,
(6)
π = 1, 2, . . . , π β 1, otherwise, (7)
We modify these functions by adding characteristics functions π[βπ,0) and π(1,1+π] to π0 and ππ , respectively, to yield a new class of extended hat functions ππ (π‘) defined over [βπ, 1+ π] for π > 0, and these are given by π0 (π‘) = π[βπ,0) (π‘) + π0 (π‘) , ππ (π‘) = ππ (π‘) ,
(8)
for 1 β€ π β€ π β 1,
(9)
ππ (π‘) = ππ (π‘) + π(1,1+π] (π‘) .
(10)
Thus, the supports of π0 (π‘) and ππ (π‘) are extended to [βπ, β] and [1 β β, 1 + π], respectively. These extended hat functions ππ (π‘) are continuous, linearly independent and are in πΏ2 [βπ, 1 + π]. As π β 0, obviously, the extended hat functions will converge to the traditional hat functions. A function π β πΏ2 [0, 1] may be approximated in vector form as π=π
π
π π (π‘) β βππ ππ (π‘) =πΉπ+1 Ξ¨π+1 (π‘) =Ξ¨π+1 (π‘) πΉπ+1 ,
(11)
π=0
where π
πΉπ+1 β [π0 , π1 , π2 , . . . , ππ ] ,
(i) π(π¦) = 0, π(π¦) = 2, π½ = 1, π = 2; (ii) π(π¦) = 2, π(π¦) = 0, πΌ = 0, π = 1, respectively, in (4). Mostly for π = 1, 2 and πΎ = 1/2 the generalized Abel integral equation models the physical problems but the integral equation for π = 2 can be reduced to the case π = 1, by change of variables. So we restrict ourselves to π = 1 only.
2. Extended Hat Functions and Their Operational Matrices for Abel Inversion
π
The important aspect of using extended hat functions in the approximation of function π(π‘) lies in the fact that the coefficients ππ in (11) are given by ππ = π (πβ) ,
π = 0, 1, 2 . . . , π.
(13)
Taking πΌ = 0, π½ = 1, and π = 2 and by change of variables, the Abel integral equation (4) reduces to π (βπ¦) π¦ π (βπ) π (βπ¦) 1 π (βπ) ππ + β« β« πΎ πΎ ππ, 2 2 0 (π¦ β π) π¦ (π β π¦) 0 β€ π¦ β€ 1, (14) which may be written as πΌ1 (π¦) = π1 (π¦) β«
π¦
0
(5)
(12)
Ξ¨π+1 (π‘) β [π0 (π‘) , π1 (π‘) , π2 (π‘) , . . . , ππ (π‘)] .
πΌ (βπ¦) =
Hat functions are defined on the domain [0, 1]. These are continuous functions with shape of hats, when plotted on two-dimensional planes. The interval [0, 1] is divided into π subintervals [πβ, (π+1)β], π = 0, 1, 2, . . . , πβ1, of equal lengths β where β = 1/π. The hat functionβs family of first (π + 1) hat functions is defined as follows: { β β π‘ , 0 β€ π‘ < β, π0 (π‘) = { β otherwise, {0,
πβ β€ π‘ < (π + 1) β,
{ π‘ β (1 β β) , 1 β β β€ π‘ β€ 1, ππ (π‘) = { β otherwise. {0,
0 < πΎ < 1, (4)
where the coefficients π(π¦) and π(π¦) do not vanish simultaneously. But the numerical inversion is still needed for its application in physical models since the experimental data for the intensity πΌ(π¦) is available only at a discrete set of points, and it may also be distorted by the noise. This motivated us to look for a stable algorithm which can be used for numerical inversion of the Abel integral equation (4) obtained by joining the two integrals (1) and (3). In this paper, we construct extended hat functions operational matrix of integration to invert the generalized Abel integral equation (4). Using hat functions for approximation of emissivity and intensity profiles has an edge over the earlier works of Singh et al. [19, 20], where they have used orthonormal Bernstein polynomials to approximate those physical quantities in the sense that a general formula for π Γ π operational matrix of integration is obtained in the earlier case whereas no such formula is available for the latter case. In Sections 3 and 4, we give the error estimate and the stability analysis followed by numerical examples to illustrate the efficiency and stability of the proposed algorithm. The above two forms (1) and (3) of Abel integral equations are obtained by taking πΎ = 1/2 and
(π β 1) β β€ π‘ < πβ,
1 π (π) π (π) πΎ ππ + π1 (π¦) β« πΎ ππ, (π¦ β π) π¦ (π β π¦)
0 β€ π¦ β€ 1, (15)
International Journal of Analysis
3
where πΌ1 (π¦) = πΌ(βπ¦), π(π) = π(βπ), π1 (π¦) = π(βπ¦)/2, and π1 (π¦) = π(βπ¦)/2. Instead of considering (15), we consider the more general equation of the form: πΌ1 (π¦) = π1 (π¦) β«
π¦
βπ
(i) When βπ β€ π¦ < (π β 1)β, then Ξ¦πΏπ (π¦) = β«
=β«
(π+1)β
πΎ ππ
ππ (π)
+β«
πβ
ππ (π)
πΎ ππ
πΎ ππ
(π β π¦)
πβ
ππ (π)
(πβ1)β
(π β π¦)
(πβ1)β
1+π
ππ (π)
(π+1)β
(π β π¦)
+β«
+β«
πΎ ππ
(π β π¦)
(π+1)β
ππ (π)
πΎ ππ
πΎ ππ
(π β π¦)
πβ
(as the support of ππ (π‘) lies in [(π β 1) β, (π + 1) β])
(17)
=β«
Thus the problem of Abel inversion is reduced to finding the unknown matrix πΆπ+1 . Substituting (17) into (16), we get
πβ
(πβ1)β
(π+1)β π β (π β 1) β (π + 1) β β π πΎ ππ + β« πΎ ππ. β(π β π¦) (π β π¦) πβ (20)
Changing the variable, π β π¦ = π‘, we get
π Ξ¨π+1 (π¦) πΉπ+1
=
ππ (π)
πβ
=β«
πΎ ππ
(π β π¦)
+β«
Using (11), the functions πΌ1 (π¦) and π(π) may be approximated as
π πΆπ+1
(πβ1)β
π¦
0 β€ π¦ β€ 1. (16)
π π (π) β πΆπ+1 Ξ¨π+1 (π) .
ππ (π)
(π β π¦)
π¦
1+π π (π) π (π) πΎ ππ + π1 (π¦) β« πΎ ππ, (π¦ β π) (π β π¦) π¦
π Ξ¨π+1 (π¦) , πΌ1 (π¦) β πΉπ+1
1+π
π¦
1+π Ξ¨π+1 (π) Ξ¨π+1 (π) [π1 (π¦) β« πΎ ππ + π1 (π¦) β« πΎ ππ] , (π β π¦) βπ (π¦ β π) π¦
Ξ¦πΏπ (π¦) =
1 2βπΎ 2βπΎ [((π β 1) β β π¦) β 2(πβ β π¦) β (1 β πΎ) (2 β πΎ) 2βπΎ
+((π + 1) β β π¦)
0 β€ π¦ β€ 1. (18) The integrals in (18) involve, evaluating integrals of the type 1+π π¦ β«π¦ ((Ξ¨π+1 (π))/(π β π¦)πΎ )ππ and β«βπ ((ππ (π))/(π¦ β π)πΎ )ππ. Let Ξ¦πΏπ (π¦) = β«
1+π
π¦
ππ (π) (π β π¦)
πΎ ππ,
(21) (ii) When (π β 1)β β€ π¦ < πβ, then Ξ¦πΏπ (π¦) = β«
πβ
π¦
π¦
ππ (π)
βπ
(π¦ β π)
Ξ¦π π (π¦) = β«
πΎ ππ
for π = 0, 1, 2 . . . , π, (19)
ππ (π)
πΎ ππ
(π β π¦)
+β«
(π+1)β
Theorem 1. The functions Ξ¦πΏπ (π¦) β πΏ2 [βπ, 1 + π] for π = 0, 1, 2, . . . , π. Proof. We prove the theorem for π = 1, 2 . . . , π β 1. The proofs for π = 0 and π = π are skipped as they may be proved on the same pattern. Based on subdivision of interval [βπ, 1 + π], we calculate Ξ¦πΏπ (π¦) by considering the following cases.
2βπΎ
ππ (π)
πΎ ππ.
(π β π¦)
πβ
(22)
Adopting the same procedure as in (i), we get 1 β (1 β πΎ) (2 β πΎ)
Ξ¦πΏπ (π¦) =
(23) 2βπΎ
and compute the two operational matrices of integration to evaluate these integrals. The scheme of derivation of these two operational matrices is based on the following theorems.
].
Γ [((π + 1) β β π¦)
2βπΎ
β 2(πβ β π¦)
].
(iii) When πβ β€ π¦ < (π + 1)β, then Ξ¦πΏπ (π¦) = β«
(π+1)β
π¦
ππ (π)
πΎ ππ
(π β π¦)
(π+1)β
=β«
π¦
(π + 1) β β π πΎ ππ β(π β π¦)
1 2βπΎ = [((π + 1) β β π¦) ] . β (1 β πΎ) (2 β πΎ)
(24)
(iv) When (π + 1)β β€ π¦ β€ 1 + π, then Ξ¦πΏπ (π¦) = 0. Hence
2βπΎ
2βπΎ
((π β 1) β β π¦) β 2(πβ β π¦) + ((π + 1) β β π¦) { { 2βπΎ 2βπΎ { { 1 ((π + 1) β β π¦) β 2(πβ β π¦) , Ξ¦πΏπ (π¦) = 2βπΎ { β (1 β πΎ) (2 β πΎ) { { {((π + 1) β β π¦) , {0,
, βπ β€ π¦ < (π β 1) β, (π β 1) β β€ π¦ < πβ, πβ β€ π¦ < (π + 1) β, otherwise.
(25)
4
International Journal of Analysis
Thus, from (25), we see that βΞ¦πΏπ (π¦)β2 < β and hence Ξ¦πΏπ (π¦) β πΏ2 [βπ, 1 + π] for π > 0 and bounded. This completes the proof.
so πππ = π½ [(π β π β 1)
Therefore, from (11), we get +
π
Ξ¦πΏπ (π¦) β βπππ ππ (π¦) ,
for π = 0, 1, 2 . . . , π.
(26)
π=0
Theorem 2. The coefficients πππ in (26) are given by
β 2(π β π)
2βπΎ
+ (2 β πΎ) (π β π)
2βπΎ
]
1 1βπΎ 1βπΎ [(1 + π β πβ) β (1 β πβ) ] . (1 β πΎ) (33) 1+π
Similarly, πππ = Ξ¦πΏπ (1) = β«1 (ππ (π)/(π β 1)πΎ ) = (1/(1 β πΎ))π1βπΎ , thus, proving the theorem.
(i) for π = 0, 1, 2, . . . , π β 1, π = 0, 1, 2 . . . , π, π½, { { { {π½ [(π β π + 1)2βπΎ β 2(π β π)2βπΎ πππ = { 2βπΎ { { { +(π β π β 1) ] , {0,
2βπΎ
for π = π, Similar arguments prove the following theorem.
(27)
for π < π, for π > π,
2 Theorem 3. The functions Ξ¦π π (π¦) β πΏ [βπ, 1 + π] for π = 0, 1, 2, . . . , π.
(ii) for π = π, π = 0, 1, 2, . . . , π,
From Theorem 3 and (11), we have
πππ
π
2βπΎ
2βπΎ
π½ [(π β π β 1) β (π β π) + (2 β πΎ) (π β π) { { { { + 1 [(1 + π β πβ)1βπΎ β (1 β πβ)1βπΎ ] , ={ (1 β πΎ) { { { 1 π1βπΎ , { (1 β πΎ)
1βπΎ
Ξ¦π π (π¦) β β πππ ππ (π¦) .
] for π < π,
The coefficients πππ βs are given as follows.
for π = π,
(i) For π = 0, 1, 2, . . . , π β 1; π = 0, 1, 2, . . . , π,
(28)
where π½ = β1βπΎ /((1 β πΎ)(2 β πΎ)).
π½, { { { { { 2βπΎ 2βπΎ { { {π½ [(π β π + 1) β 2(π β π) πππ = { 2βπΎ { { +(π β π β 1) ] , { { { { { {0,
Proof. (i) When π = π, πππ = Ξ¦πΏπ (πβ) = π½ follows from (25). When π < π, which is equivalent to π β€ π β 1, we get πππ = Ξ¦πΏπ (πβ) = β«
1+π
πβ
= (β«
πβ
(πβ1)β
+β«
ππ (π) (π β πβ)
(π+1)β
πβ
πΎ ππ
2βπΎ
)
ππ (π) (π β πβ)
πΎ ππ,
2βπΎ
β 2(π β π)
1
1+π
πβ
1
πππ = Ξ¦πΏπ (πβ) = (β« + β«
2βπΎ
+ (π β π β 1)
(πβ1)β
for π < π,
)
ππ (π) (π β πβ)
(35)
(ii) for π = π; π = 0, 1, 2 . . . , π, 1
for π = 0, π1βπΎ , { { { (1 β πΎ) { { { 2βπΎ 2βπΎ 1βπΎ πππ = {π½ [(π β 1) β π + (2 β πΎ) π ] { { { { { + 1 [(π + πβ)1βπΎ β (πβ)1βπΎ ] , for 1 β€ π β€ π. (1 β πΎ) { (36)
] . (30)
πΎ ππ.
(31) Using (12) and (26), the following integral may be written as
From (10), we have πβ
for π > π,
and
Similarly, for π > π, that is, π β₯ π + 1, πππ = 0 (it follows trivially from (19)). (ii) For π < π,
πππ = β«
for π = π,
(29)
since the support of ππ (π‘) lies in [(π β 1)β, (π + 1)β]. Using (6) and (9), we get from change of variable πππ = π½ [(π β π + 1)
(34)
π=0
1+π 1 (π β (π β 1) β) πΎ ππ + β« πΎ ππ, β(π β πβ) (π β πβ) 1
(32)
1+π
β«
π¦
π Ξ¨π+1 (π) πΏ πΏ πΏ πΎ ππ = [Ξ¦0 (π¦) , Ξ¦1 (π¦) , . . . , Ξ¦π (π¦)] . (π β π¦)
(37)
International Journal of Analysis
5
Substituting the approximation of Ξ¦πΏπ (π¦) from (26) in (37), we get 1+π
β«
π¦
Ξ¨π+1 (π) πΏ πΎ ππ = ππ+1 Ξ¨π+1 (π¦) , (π β π¦)
where π = ππ+1
(38)
π0 [0 [ [0 [ [0 Γ[ [0 [β
β
β
[ [0 [0
πΏ where ππ+1 is a (π + 1) Γ (π + 1) matrix whose (π + 1, π + 1)th entry is πππ , given by (27) and (28) for π, π = 0, 1, 2, . . . , π. The πΏ matrix ππ+1 is given as
πΏ ππ+1 =
β1βπΎ (1 β πΎ) (2 β πΎ) 1 0 1 [ π1 [ [ π2 π1 [ π2 [π Γ[ 3 [ π4 π3 [ β
β
β
β
β
β
[ [π πβ1 ππβ2 [ ππ ππβ1
0 0 1 π1 π2 β
β
β
β
β
β
β
β
β
0 0 0 1 π1 β
β
β
π3 π4
0 0 0 0 1 β
β
β
π2 π3
0 0 0 0 0 1 π1 π2
0 0 0 0 0 0 1 π1
0 0] ] 0] ] 0] , 0] ] ] 0] 0] π0 ](π+1)Γ(π+1) (39)
where π0 =
(2 β πΎ) 1βπΎ π , β1βπΎ
(40)
ππ = [(π β 1)2βπΎ β π 2βπΎ + (2 β πΎ) π 1βπΎ ] +
π1 1 0 0 0 β
β
β
β
β
β
β
β
β
π2 π1 1 0 0 β
β
β
β
β
β
β
β
β
π3 π2 π1 1 0 β
β
β
β
β
β
β
β
β
π4 π3 π2 π1 1 β
β
β
β
β
β
β
β
β
β
β
β
β
β
β
π3 π2 π1 β
β
β
0 β
β
β
ππβ1 ππβ2 β
β
β
β
β
β
β
β
β
β
β
β
1 0
ππ ππβ1 ] ] ππβ2 ] ] ππβ3 ] . ππβ4 ] ] ] β
β
β
] π1 ] 1 ](π+1)Γ(π+1) (44)
The various entries π0 , ππ , and ππ are given by (40)β(42), for π = 1, 2, 3, . . . , π β 1 and π = 1, 2, 3, . . . , π. πΏ in four blocks as If we partition the matrix ππ+1 πΏ ππ+1
=
π [ π΄π΅ ππ0 ], then ππ+1
=
σΈ
π΅σΈ ], π π΄π
[ π0
where π΅
=
[ ππ β
β
β
π2 π1 ]1Γπ , π΅ = [ π1 π2 β
β
β
ππ ]1Γπ , π is a π Γ 1 null vector, and 1 0 [ π1 1 [ [ π2 π1 π΄=[ [ π3 π2 [ [ππβ2 β
β
β
[ππβ1 ππβ2
0 0 1 π1 β
β
β
β
β
β
0 0 0 1 π1 β
β
β
0 0 0 0 1 π1
0 0] ] 0] ] . 0] ] 0] 1]πΓπ
(45)
Using (38) and (43), (18) may be written as
(2 β πΎ) [(π + π β)1βπΎ β (π β)1βπΎ ] , β1βπΎ
ππ = (π + 1)2βπΎ β 2π2βπΎ + (π β 1)2βπΎ ,
β1βπΎ (1 β πΎ) (2 β πΎ)
π = 1, 2, 3, . . . , π, (41)
π = 1, 2, 3, . . . , π β 1. (42)
π π π πΏ πΉπ+1 Ξ¨π+1 (π¦) = πΆπ+1 [π (π¦) ππ+1 + π (π¦) ππ+1 ] Ξ¨π+1 (π¦) . (46)
Solving the above system of linear equations, we obtain β1
πΏ The matrix ππ+1 is called extended hat functions lower operational matrix for Abelβs inversion.
π π π πΏ πΆπ+1 = πΉπ+1 + π (π¦) ππ+1 [π (π¦) ππ+1 ] .
Remark 4. It is evident from (40) that when π = 0, then π0 = πΏ becomes a singular 0, and so the lower triangular matrix ππ+1 πΏ makes it matrix. In this case, the singularity of the matrix ππ+1 redundant for numerical computation since the invertibility of the matrix is required to obtain the solution. To make the πΏ invertible, we introduced a positive parameter π matrix ππ+1 and extend the traditional hat function to the interval [βπ, 1+ π]. Similarly, using (19) and (34)β(36), we construct the extended hat functions upper operational matrix for Abelβs π , such that inversion, ππ+1
π Substituting the value of πΆπ+1 from (47) into (17), the approximate emissivity π(π) is given by
β«
π¦
βπ
Ξ¨π+1 (π) π πΎ ππ = ππ+1 Ξ¨π+1 (π¦) , (π¦ β π)
(43)
β1
π π πΏ π (π) = πΉπ+1 + π (π¦) ππ+1 [π (π¦) ππ+1 ] Ξ¨π+1 (π) .
(47)
(48)
3. Error Analysis In this section, an error analysis of our proposed algorithm is given. Let π β πΏ2 [0, 1], and then, using (11), it is approximated as π=π
π=π
π=0
π=0
π (π‘) β βππ ππ (π‘) = βπ (πβ) ππ (π‘) .
(49)
6
International Journal of Analysis
The above approximation gives exact values at nodal points. Μ and then, for We denote the right-hand side of (49) by π(π‘), πβ β€ π‘ < (π + 1)β, π = 0, 1, . . . , π β 1, we have π=π
πΜ (π‘) = βπ (πβ) ππ (π‘) = π (πβ) ππ (π‘) + π ((π + 1) β) ππ+1 (π‘) π=0
= π (πβ) (
(π + 1) β β π‘ π‘ β πβ ) + π (πβ + β) ( ) β β
the exact results are known, and thus making a comparison between inverted results and theoretical data is possible. We have tested our algorithm on several well-known test profiles that are commonly encountered in experimental data and widely used by researchers [7, 15, 20]. The accuracy of the proposed algorithm is demonstrated by calculating the parameters of absolute error Ξπ(ππ ), average deviation π also known as root mean square error (RMS). They are calculated using the following equations: σ΅¨ σ΅¨ Ξπ (ππ ) = σ΅¨σ΅¨σ΅¨ππ (ππ ) β ππ (ππ )σ΅¨σ΅¨σ΅¨ ,
(using (6)) = (π + 1) π (πβ) + π‘
(π (πβ + β) β π (πβ)) β
1/2
ππ+1
β ππ (πβ + β) .
π 1 2 ={ β[ππ (ππ ) β ππ (ππ )] } (π + 1) π=0
(50) ={
Expanding π(π‘) in Taylor series at π‘ = πβ, we obtain π
β
(π‘ β πβ) (π) π (πβ) . π! π=0
π (π‘) = β
(51)
Thus, from (50) and (51), we get σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨π (π‘) β πΜ (π‘)σ΅¨σ΅¨σ΅¨ = σ΅¨ σ΅¨
σ΅¨σ΅¨ (π (πβ + β) β π (πβ)) σ΅¨σ΅¨ σ΅¨σ΅¨πβ σ΅¨σ΅¨ β σ΅¨ (π (πβ + β) β π (πβ)) β σ΅¨σ΅¨ π β σ΅¨σ΅¨ (π‘ β πβ) (π) +β π (πβ)σ΅¨σ΅¨σ΅¨ σ΅¨σ΅¨ π! π=1 σ΅¨
βπ‘
(52)
as β β 0, and we have σ΅¨ σ΅¨σ΅¨ π β σ΅¨σ΅¨ (π‘ β πβ) (π) σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨π (π‘) β πΜ (π‘)σ΅¨σ΅¨σ΅¨ = σ΅¨σ΅¨σ΅¨σ΅¨πβπσΈ (πβ) β π‘πσΈ (πβ) + β π (πβ)σ΅¨σ΅¨σ΅¨ σ΅¨ σ΅¨σ΅¨ σ΅¨ σ΅¨σ΅¨ π! π=1 σ΅¨ σ΅¨ σ΅¨σ΅¨ β σ΅¨ π σ΅¨σ΅¨ σ΅¨σ΅¨ (π‘ β πβ) (π) σ΅¨ = σ΅¨σ΅¨σ΅¨β π (πβ)σ΅¨σ΅¨σ΅¨ σ΅¨σ΅¨π=2 π! σ΅¨σ΅¨ σ΅¨ σ΅¨ 2
β€
(54)
(π‘ β πβ) πβ2 σ΅¨ 3σ΅¨ + π σ΅¨σ΅¨σ΅¨σ΅¨(π‘ β πβ) σ΅¨σ΅¨σ΅¨σ΅¨ β€ + π (πβ3 ) , 2! 2 (53)
as (π‘ β πβ) < πβ1 , thus proving the following theorem. Μ Theorem 5. The absolute error |π(π‘) β π(π‘)| associated with β2 the approximation (49) is of the order π(π ).
4. Numerical Results and Stability Analysis In this section, we discuss the implementation of our proposed algorithm and investigate its accuracy and stability by applying it on test functions with known analytical Abel inverse. For, it is always desirable to test the behaviour of a numerical inversion method using simulated data for which
1/2
π
1 2 β[Ξπ (ππ )] } (π + 1) π=0
(55) σ΅© σ΅© = σ΅©σ΅©σ΅©Ξπσ΅©σ΅©σ΅©2 ,
where ππ (ππ ) is the emission coefficient calculated at point ππ using (48) and ππ (ππ ) is the exact analytical emissivity at the corresponding point. Note that π, henceforth, denoted by ππ+1 (for computational convenience) is the discrete π2 norm of the absolute error Ξπ denoted by βΞπβ2 . Note that the calculation of ππ+1 in (55) is performed by taking different values of π. In all the test profiles, the exact and noisy intensity profiles are denoted by πΌ1 (π¦) and πΌ1πΏ (π¦), respectively, where πΌ1πΏ (π¦) is obtained by adding a noise πΏ to πΌ1 (π¦) such that πΌ1πΏ (π¦π ) = πΌ1 (π¦π ) + πΏππ , where π¦π = πβ, π = 0, 1, 2, . . . , π, πβ = 1, and ππ is the uniform random variable with values in [β1, 1] such that Max0β€πβ€π |πΌ1πΏ (π¦π ) β πΌ1 (π¦π )| β€ πΏ. The following test problems are solved with and without noise to illustrate the efficiency and stability of our method by choosing three different values of the noises πΏπ as πΏ0 = 0, πΏ1 = ππ+1 , and πΏ2 = π% of ππ+1 , where we mean ππ+1 = (1/(π + 1)) βππ=0 πΌ1 (π¦). In each of the test problems given in this section we have taken positive parameter π = 0.0001 and π = 0.1, except for Example 10, where π = 0 has been used. The absolute errors between exact and approximate emissivities, corresponding to different noises πΏπ , π = 0, 1, 2 have been denoted by πΈ0 (π), πΈ1 (π), and πΈ2 (π), respectively. In the text boxes of the figures, the notations πΈ0(π), πΈ1(π), and πΈ2(π) have been used for πΈ0 (π), πΈ1 (π), and πΈ2 (π), respectively. Though the stability of the algorithm is illustrated by various numerical experiments performed in this section, we analyze it also mathematically as follows. The reconstructed emissivities πππΏ (π) (with πΏ noise) and ππ0 (π) (without noise) are obtained with and without noise term in the intensity profile πΌ1 (π¦), and using (48) these are given by β1
πΏπ π πΏ πππΏ (π) = πΉπ+1 [π (π¦) ππ+1 + π (π¦) ππ+1 ] Ξ¨π+1 (π) ,
ππ0 (π) =
π π πΉπ+1 [π (π¦) ππ+1
+
β1 πΏ π (π¦) ππ+1 ] Ξ¨π+1
(56) (π) ,
International Journal of Analysis
7
Table 1: Noise reduction β(π) for π = 100 at different values of πΏ, π = 0, πΎ = 1/2.
Γ10β4 2
π 0
β(π) for πΏ = 0.01 0.019718
β(π) for πΏ = 0.001 0.0019718
β(π) for πΏ = 0.0001 0.00019718
0.1
0.0040696
0.00040696
0.000040696
0.2
0.0035401
0.00035401
0.000035401
0.3
0.0033118
0.00033118
0.000033118
0.4
0.0032047
0.00032047
0.000032047
1
0.5
0.0031726
0.00031726
0.000031726
0.8
0.6
0.0032047
0.00032047
0.000032047
0.6
0.7
0.0033118
0.00033118
0.000033118
0.8
0.0035401
0.00035401
0.000035401
0.9
0.0040696
0.00040696
0.000040696
0.019718
0.0019718
0.00019718
1
(π¦) = πΌ1 (π¦) + πΏππ β
πΏπ πΉπ+1 Ξ¨π+1
1.6
β(π)
1.4 1.2
0.4 0.2
0
0.2
0.4
0.6
0.8
1
π
πΏ where πΉπ+1 and πΉπ+1 are known matrices, and they are obtained from the following equations:
πΌ1πΏ
1.8
(π¦) , (57)
π Ξ¨π+1 (π¦) . πΌ1 (π¦) β πΉπ+1
Figure 1: Noise reduction β(π) for π = 100, πΏ = 0.0001, π = 0, and πΎ = 1/2.
Example 6. Consider the generalized Abel integral equation: 4 3/2 2 π¦ (1 + π¦) + π¦3/2 [2β(1 β π¦) π¦3 + βπ¦ β π¦2 ] 3 3 = π (π¦) β«
Hence
π¦
0
π (π) β(π¦ β π)
1
π (π)
π¦
β(π β π¦)
ππ + π (π¦) β«
πππΏ (π) β ππ0 (π)
ππ,
(60)
0 β€ π¦ β€ 1, β1
πΏπ π π πΏ β πΉπ+1 ) [π (π¦) ππ+1 + π (π¦) ππ+1 ] Ξ¨π+1 (π) . = (πΉπ+1 (58) π πΏπ π = πΉπ+1 β πΉπ+1 and replacing random noise πΏππ Writing π»π+1 by its maximum value πΏ, we get
σ΅© σ΅©σ΅© σ΅©σ΅©πππΏ (π) β ππ0 (π)σ΅©σ΅©σ΅© σ΅© σ΅© β1 σ΅© σ΅© π σ΅©σ΅© σ΅©σ΅©σ΅© σ΅© σ΅©σ΅© π πΏ σ΅©σ΅© σ΅©σ΅©[π (π¦) ππ+1 = σ΅©σ΅©σ΅©σ΅©π»π+1 + π (π¦) ππ+1 ] σ΅©σ΅©σ΅© σ΅©σ΅©σ΅©σ΅©Ξ¨π+1 (π)σ΅©σ΅©σ΅©σ΅© . σ΅©σ΅© σ΅©
(59)
πΏπ π π Let β(π) = πππΏ (π) β ππ0 (π) = (πΉπ+1 β πΉπ+1 )[π(π¦)ππ+1 + πΏ β1 π(π¦)ππ+1 ] Ξ¨π+1 (π), then β(π) reflects the noise reduction capability of the algorithm and its values at various points, and its graph is shown in Table 1 and Figure 1, respectively. Table 1 demonstrates the noise filtering capability of the algorithm for three different noise outputs. From Table 1 and Figure 1 we see that noise reduction is symmetric about the point π = 0.5, and the maximum reduction in noise is achieved at π = 0.5 for all the three levels of noises πΏ = 0.01, 0.001, and 0.0001 introduced in πΌ1 (π¦). The general behaviour of the noise reduction is the same irrespective of the value of πΏ. In the interval [0.02, 0.98] the algorithm is stable, whereas the noise filtering capability decreases continuously and then jumps symmetrically in [0, 0.02) βͺ (0.98, 1].
where π(π¦) = π¦ + 1, π(π¦) = π¦2 , with the exact analytical solution π(π) = π. The absolute errors πΈπ (π) have been calculated for π = 1000 and are given in Table 2. The value of πΏ1 is 2.3829Γ10β15 , for π = 1000. As πΏ2 = 5.3347 Γ 10β4 > 1010 πΏ1 , the absolute error πΈ2 (π) is appreciably higher than πΈ0 (π) and πΈ1 (π). The Figure 2 compares the absolute errors πΈ0 (π) and πΈ1 (π) for noise πΏ1 = 5.9827 Γ 10β16 , π = 100. Example 7. In this example, we consider the following Abelβs integral equation [7, 20]: π¦
β«
0
π (π) β(π¦ β π)
ππ = π¦11/2 ,
with solution π (π) =
0 β€ π¦ β€ 1, (61)
10
2
2 11[Ξ (11/2)] 5 π. 2πΞ (11)
The absolute errors corresponding to different noises are given in Table 3. The values of various parameters are given as: π3001 = 9.4002 Γ 10β8 (= πΏ1 , π = 3000), π2001 = 2.1116 Γ β7 10 (= πΏ1 , π = 2000), π1001 = 8.4162 Γ 10β7 (= πΏ1 , π = 1000), π3001 = 0.1540, π2001 = 0.1540, and π1001 = 0.1542. Taking π = 0.1, the various values of respective πΏ2 are given in order as 1.5396 Γ 10β4 , 1.5402 Γ 10β4 , and 1.5419 Γ 10β4 .
8
International Journal of Analysis Table 2: The absolute errors πΈπ (π), at different nodal points π, for π = 1000, in Example 6.
π πΈ0 (π) πΈ1 (π) πΈ2 (π)
0.0 0 1.132 Γ 10β13 0.011106
0.2 7.494 Γ 10β16 2.4591 Γ 10β14 0.0064263
0.4 1.2212 Γ 10β15 8.3822 Γ 10β15 0.0025148
0.6 1.8874 Γ 10β15 1.7764 Γ 10β14 0.01039
0.8 8.1046 Γ 10β15 2.931 Γ 10β14 0.00029049
1.0 1.5543 Γ 10β15 3.5194 Γ 10β14 0.00014606
Table 3: The absolute errors πΈπ (π), at different nodal points π, for π = 3000, 2000, and 1000. π 3000 2000 1000 3000 2000 1000 3000 2000 1000
πΈ0 (π)
πΈ1 (π)
πΈ2 (π)
π = 0.0 0 0 0 4.6073 Γ 10β6 1.0031 Γ 10β5 1.7061 Γ 10β5 0.0061012 0.0068177 0.0001341
π = 0.2 1.9712 Γ 10β9 4.4178 Γ 10β7 1.7515 Γ 10β8 1.5029 Γ 10β6 1.1159 Γ 10β6 7.7631 Γ 10β6 0.0011742 0.0068177 0.00067286
π = 0.4 1.5849 Γ 10β8 3.5563 Γ 10β8 1.4137 Γ 10β7 2.5454 Γ 10β6 4.5957 Γ 10β6 7.528 Γ 10β6 0.0011595 0.00078732 0.0010734
Γ10β14 Exact and approximate emissivity
1.6 1.4 πΈ0 (π) and πΈ1 (π)
π = 0.8 1.2724 Γ 10β7 2.8575 Γ 10β7 1.138 Γ 10β6 7.3491 Γ 10β7 4.4563 Γ 10β6 1.5261 Γ 10β5 0.00064593 0.0016351 0.00098541
π = 1.0 2.4875 Γ 10β7 5.5872 Γ 10β7 2.2262 Γ 10β6 2.0451 Γ 10β6 9.2247 Γ 10β7 4.2806 Γ 10β6 0.0029774 4.2642 Γ 10β5 0.00074723
1.4
1.8
1.2 1 0.8 0.6 0.4
1.2 1 0.8 0.6 0.4 0.2 0 β0.2
0.2 0
π = 0.6 5.361 Γ 10β8 1.2035 Γ 10β7 4.7899 Γ 10β7 5.3295 Γ 10β7 2.3615 Γ 10β7 4.4903 Γ 10β6 0.0039968 0.00067568 0.00040813
0
0.2
0.4
0.6
0.8
1
π 0
0.2
0.4
0.6
0.8
1
Exact emissivity Approximate emissivity with noise πΏ2
π 10 Γ πΈ0 (π) πΈ1 (π)
Figure 2: Absolute errors πΈ0 (π) and πΈ1 (π) for π = 100, πΏ1 = 5.9827 Γ 10β16 .
In Figure 3, the exact and reconstructed emissivities (with πΏ2 noise) have been shown for π = 50, and the two emissivities match very well even for higher noise πΏ2 introduced in the intensity profile. For π = 100, Figures 4 and 5 show the absolute errors πΈ0 (π), πΈ1 (π) and πΈ0 (π), πΈ2 (π), respectively. Example 8. In this example we consider the following Abelβs integral equation [22]: β«
π¦
0
π (π) β(π¦ β π)
ππ = πΌ1 (π¦) ,
0 β€ π¦ β€ 1,
(62)
Figure 3: Exact and approximate emissivities with noise πΏ2 = 1.6081 Γ 10β4 for π = 50.
where 4 3/2 π¦ , { { {3 πΌ1 (π¦) = { { 4 3/2 8 1 3/2 { π¦ β (π¦ β ) , {3 3 2
1 0β€π¦< , 2 1 β€ π¦ β€ 1. 2
(63)
The exact solution of the integral equation (62) is given by { {π, π (π) = { {1 β π, {
1 0β€π< , 2 1 β€ π β€ 1. 2
(64)
In Table 4, the absolute errors for different noises have been shown. Various parameters used for Table 4 are
International Journal of Analysis
9
Table 4: The absolute errors πΈπ (π), at different nodal points π, for π = 1000, in Example 8. 0.0 0
0.2 4.4409 Γ 10β16
0.4 5.5511 Γ 10β17
0.6 2.7200 Γ 10β15
0.8 2.1649 Γ 10β15
1.0 3.3888 Γ 10β15
πΈ1 (π)
5.5359 Γ 10β14
7.8271 Γ 10β15
2.2093 Γ 10β14
1.1768 Γ 10β14
1.4433 Γ 10β14
1.5344 Γ 10β14
πΈ2 (π)
0.0042764
0.0018814
3.4439 Γ 10β5
0.0014343
0.0024313
0.00018694
Γ10β3 1.2
Γ10β3 3.5
1
3 2.5
0.8
πΈ0 (π) and πΈ2 (π)
πΈ0 (π) and πΈ1 (π)
π πΈ0 (π)
0.6 0.4
1.5 1
0.2 0
2
0.5 0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0
1
0
0.2
0.4
πΈ0 (π) πΈ1 (π)
1
Example 9. Consider the generalized Abel integral equation: πΌ1 (π¦) = ππ¦ sin (π¦) β«
π¦
0
π (π)
ππ 1/3 (π¦ β π)
+ πβπ¦ cos (π¦) β«
1
π¦
π¦
π (π)
ππ, 1/3 (π β π¦)
5/3
(65) 0 < π¦ < 1,
Example 10. For the pair [14, 15, 23]: π (π) = (1 β π) (1 + 12π) ,
0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1
0
where πΌ1 (π¦) = π sin (π¦) π¦ + π cos(π¦)(1 β π¦) (2 + 3π¦). The exact solution of (65) is π(π) = 10π/9. In Figure 9, the comparison between πΈ0 (π) and πΈ1 (π) is shown, for π = 100.
2
0.5
0.05
2/3
βπ¦
Figure 5: Absolute errors πΈ0 (π) and πΈ2 (π) for πΏ2 = 1.5732 Γ 10β4 , π = 100.
Exact and approximate emissivity
π1001 = 1.8576Γ10β15 = πΏ1 , π1001 = 0.3446, and πΏ2 = 3.4462Γ 10β4 . Figure 6 shows the graph of exact and approximate emissivities π(π) (without noise) for π = 50. Absolute errors πΈ0 (π) and πΈ1 (π), for π = 50 and π = 100, are shown in Figures 7 and 8, respectively.
0
0.2
0.4
0.6
0.8
1
π Approximate emissivity Exact emissivity
Figure 6: Exact and approximate solutions for π = 50.
for 0 β€ π β€ 1,
384 7/2 368 5/2 40 3/2 π¦ β π¦ + π¦ + 2βπ¦) 35 15 3 +
0.8
πΈ0 (π) πΈ2 (π)
Figure 4: Absolute errors πΈ0 (π) and πΈ1 (π) for πΏ1 = 8.2435 Γ 10β5 , π = 100.
πΌ1 (π¦) = [(
0.6 π
π
16 5/2 (1 β π¦) (19 + 72π¦)] , 105
(66)
with π(π¦) = π(π¦) = 1, and the various parameters are as follows: πΏ1 = 2.1111 Γ 10β7 (π = 3000), πΏ1 = 4.7384 Γ β7 10 (π = 2000), πΏ1 = 1.8852 Γ 10β6 (π = 1000), and πΏ2 = 0.0036, for all the three chosen values of π.
10
International Journal of Analysis Table 5: The absolute errors πΈπ (π), at different nodal points π, for π = 3000, 2000, and 1000.
πΈ0 (π)
πΈ1 (π)
πΈ2 (π)
π 3000
π = 0.0 5.4281 Γ 10β7
π = 0.2 2.8998 Γ 10β7
π = 0.4 1.5818 Γ 10β7
π = 0.6 2.6121 Γ 10β8
π = 0.8 1.0578 Γ 10β7
π = 1.0 2.9684 Γ 10β7
2000
1.2088 Γ 10β6
6.5113 Γ 10β7
3.5535 Γ 10β7
5.887 Γ 10β8
2.3717 Γ 10β7
6.5839 Γ 10β7
1000 3000
4.7364 Γ 10β6 6.8261 Γ 10β6
2.5925 Γ 10β6 2.1337 Γ 10β6
1.4164π Γ 10β6 2.2043 Γ 10β6
2.3637 Γ 10β7 1.4765 Γ 10β6
9.4123 Γ 10β7 1.4642 Γ 10β6
2.559 Γ 10β6 3.421 Γ 10β6
2000
1.589 Γ 10β6
5.7191 Γ 10β6
2.8577 Γ 10β6
2.5629 Γ 10β6
4.4183 Γ 10β6
3.9074 Γ 10β8
β6
3.6477 Γ 10β6 0.095356
β5
β6
β5
β5
1000 3000
1.9156 Γ 10 0.01495
4.4716 Γ 10 0.020628
1.643 Γ 10 0.021094
1.712 Γ 10 0.065933
3.6022 Γ 10 0.011844
2000
0.022427
0.010769
0.024946
0.014392
0.0068535
0.028831
1000
0.056845
0.033342
0.013319
0.018905
0.028086
0.045314
Γ10β14 1.4
Γ10β14 2.5
1.2
2 πΈ0 (π) and πΈ1 (π)
πΈ0 (π) and πΈ1 (π)
1 0.8 0.6
1.5
1
0.4 0.5 0.2 0
0
0.2
0.4
0.6
0.8
1
0
0
0.2
0.6
0.4
π
0.8
1
π 10 Γ πΈ0 (π) πΈ1 (π)
10 Γ πΈ0 (π) πΈ1 (π)
Figure 7: Absolute errors πΈ0 (π) and πΈ1 (π) for πΏ1 = 3.4218 Γ 10β16 , π = 50.
Figure 8: Absolute errors πΈ0 (π), πΈ1 (π) for πΏ1 = 5.5794 Γ 10β16 , π = 100.
The absolute errors corresponding to different noises are given in Table 5. Figure 10 shows the exact and approximate emissivities (without noise and with noise πΏ2 = 0.0036) whereas, in Figure 11, a comparison between πΈ0 (π) and πΈ1 (π) is shown for πΏ1 = 1.8232 Γ 10β4 , π = 100.
where πΆ(π§) and π(π§) in (67) are called Fresnel integrals. These are defined as
Example 11. Consider the generalized Abel integral equation (15) with πΎ = 1/2, π(π¦) = (3/4) exp(π¦), π(π¦) = exp(2π¦) + (1/β2π), for the pair π(π) = sin π and πΌ1 (π¦) 5 7 π¦2 1 exp (π¦) π¦13/2 πΉ2 (1, , , β ) + β2π (exp (2π¦) + ) [ β2π ] 4 4 4 ], =[ ] [ β2 β 2π¦ β2 β 2π¦ ) sin π¦ + π ( ) cos π¦) (πΆ ( βπ βπ ] [
(67)
π§
πΆ (π§) = β« cos ( 0
ππ‘2 ) ππ‘, 2
π§
π (π§) = β« sin ( 0
ππ‘2 ) ππ‘. 2 (68)
For, π = 3000, different absolute errors are given in Table 6. The various parameters for π = 3000 are πΏ1 = 0.0173 and πΏ2 = 7.2830 Γ 10β4 . In Figure 12, the exact and approximate emissivities (without noise) are shown for π = 100.
5. Conclusions We have constructed operational matrices of integration based on extended hat functions and used them to propose a stable algorithm for the numerical inversion of the generalized Abel integral equation. The earlier numerical
International Journal of Analysis
11 Table 6: The absolute errors πΈπ (π) for π = 3000, for Example 11.
π πΈ0 (π) πΈ1 (π) πΈ2 (π)
0.0 3.389 Γ 10β7 0.38495 0.0066376
0.2 4.7548 Γ 10β7 0.079121 0.0058214
0.4 7.3287 Γ 10β7 0.17527 0.0070449
0.6 1.3451 Γ 10β6 0.096373 0.018136
Γ10β13 7
Γ10β4 8
6
7
πΈ0 (π) and πΈ1 (π)
πΈ0 (π) and πΈ1 (π)
4 3 2
5 4 3 2
1
1 0
0.2
0.4
0.6
0.8
0
1
0
0.2
0.4
π
0.8
1
πΈ0 (π) πΈ1 (π)
Figure 9: Comparison between πΈ0 (π) and πΈ1 (π) for πΏ1 = 1.6911 Γ 10β15 , π = 100.
Figure 11: Comparison between πΈ0 (π) and πΈ1 (π) for πΏ1 = 1.8232 Γ 10β4 , π = 100.
2.5
0.8 Exact and approximate emissivity
2 1.5 1 0.5 0 β0.5
0.6 π
100 Γ πΈ0 (π) πΈ1 (π)
Exact and approximate emissivity
0.9 1.0734 Γ 10β5 0.050512 0.0021626
6
5
0
0.8 3.7972 Γ 10β6 0.40061 0.013987
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
0
0.2
0.4
0.6
0.8
0
0.1
0.2
0.3
0.4
1
π Approximate emissivity without noise Exact emissivity Approximate emissivity with πΏ2 noise
Figure 10: Exact and approximate emissivities (without noise and with noise πΏ2 = 0.0036) for π = 50.
inversions were restricted to a part of the general Abelβs integral equation. We have extended the hat function beyond their domain [0, 1] to avoid the singularity of the matrix at π‘ = 0, 1. These operational matrices are given by the general formulae (39) and (44), thus making the evaluation of
0.5
0.6
0.7
0.8
0.9
π Approximate emissivity Exact emissivity
Figure 12: Exact and approximate emissivities (without noise) for π = 100.
these matrices of any order extremely easy whereas in case of Bernstein operational matrices no such formula was available [19, 20]. The stability with respect to the data is restored and good accuracy is obtained even for high noise levels in the data. An error analysis and stability analysis are also given.
12
Acknowledgments The authors are grateful to the learned reviewer for his valuable suggestions which have led to the improvement of the paper in the present form. Also, the first author acknowledges the financial support from UGC New-Delhi, India, under Faculty Improvement Program (FIP), whereas the second and the third authors acknowledge the financial support from UGC and CSIR New-Delhi, India, respectively, under JRF schemes.
References [1] N. H. Abel, βResolution dβun problem de mechanique,β Journal FΒ¨ur Die Reine Und Angewandte Mathematik, vol. 1, pp. 153β157, 1826. [2] L. Mach, Wien. Akad. Ber. Math. Phys. Klasse, vol. 105, p. 605, 1896. [3] H. R. Griem, Plasma Spectroscopy, McGraw-Hill, New York, NY, USA, 1963. [4] F. G. Tricomi, Integral Equations, Wiley-Interscience, New York, NY, USA, 1975. [5] R. S. Anderssen, βStable procedures for the inversion of Abelβs equation,β Journal of the Institute of Mathematics and its Applications, vol. 17, no. 3, pp. 329β342, 1976. [6] I. Beniaminy and M. Deutsch, βABEL: stable, high accuracy program for the inversion of Abelβs integral equation,β Computer Physics Communications, vol. 27, no. 4, pp. 415β422, 1982. [7] S. Bhattacharya and B. N. Mandal, βUse of Bernstein polynomials in numerical solutions of Volterra integral equations,β Applied Mathematical Sciences, vol. 2, no. 33β36, pp. 1773β1787, 2008. [8] C. J. Cremers and R. C. Birkebak, βApplication of the Abel integral equation to spectroscopic data,β Applied Optics, vol. 5, pp. 1057β1064, 1966. [9] M. Deutsch and I. Beniaminy, βDerivative-free inversion of Abelβs integral equation,β Applied Physics Letters, vol. 41, no. 1, pp. 27β28, 1982. [10] M. Deutsch, βAbel inversion with a simple analytic representation for experimental data,β Applied Physics Letters, vol. 42, no. 3, pp. 237β239, 1983. [11] R. Grenflo, βComputation of rough solutions of Abel integral equation,β in Inverse ILL-Posed Problems, H. W. Engel and C. W. Groetsch, Eds., pp. 195β210, Academic Press, NewYork, NY, USA, 1987. [12] L. M. IgnjatoviΒ΄c and A. A. Mihajlov, βThe realization of Abelβs inversion in the case of discharge with undetermined radius,β Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 72, no. 5, pp. 677β689, 2002. [13] J.-P. Lanquart, βError attenuation in Abel inversion,β Journal of Computational Physics, vol. 47, no. 3, pp. 434β443, 1982. [14] S. Ma, H. Gao, L. Wu, and G. Zhang, βAbel inversion using Legendre polynomials approximations,β Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 109, no. 10, pp. 1745β1757, 2008. [15] S. Ma, H. Gao, G. Zhang, and L. Wu, βAbel inversion using Legendre wavelets expansion,β Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 107, no. 1, pp. 61β71, 2007. [16] G. N. Minerbo and M. E. Levy, βInversion of Abelβs integral equation by means of orthogonal polynomials,β SIAM Journal on Numerical Analysis, vol. 6, pp. 598β616, 1969.
International Journal of Analysis [17] D. A. Murio, D. G. Hinestroza, and C. E. MejΒ΄Δ±a, βNew stable numerical inversion of Abelβs integral equation,β Computers and Mathematics with Applications, vol. 23, no. 11, pp. 3β11, 1992. [18] M. Sato, βInversion of the Abel integral equation by use of simple interpolation formulas,β Contributions to Plasma Physics, vol. 25, pp. 573β577, 1985. [19] V. K. Singh, R. K. Pandey, and O. P. Singh, βNew stable numerical solutions of singular integral equations of Abel type by using normalized Bernstein polynomials,β Applied Mathematical Sciences, vol. 3, no. 5-8, pp. 241β255, 2009. [20] O. P. Singh, V. K. Singh, and R. K. Pandey, βA stable numerical inversion of Abelβs integral equation using almost Bernstein operational matrix,β Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 111, pp. 245β252, 2010. [21] A. Chakrabarti, βSolution of the generalized Abel integral equation,β Journal of Integral Equations and Applications, vol. 20, no. 1, pp. 1β11, 2008. [22] L. Huang, Y. Huang, and X.-F. Li, βApproximate solution of Abel integral equation,β Computers and Mathematics with Applications, vol. 56, no. 7, pp. 1748β1757, 2008. [23] M. J. Buie, J. T. P. Pender, J. P. Holloway, T. Vincent, P. L. G. Ventzek, and M. L. Brake, βAbelβs inversion applied to experimental spectroscopic data with off axis peaks,β Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 55, no. 2, pp. 231β243, 1996.
Advances in
Operations Research Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Advances in
Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Applied Mathematics
Algebra
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Probability and Statistics Volume 2014
The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Differential Equations Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Submit your manuscripts at http://www.hindawi.com International Journal of
Advances in
Combinatorics Hindawi Publishing Corporation http://www.hindawi.com
Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Journal of
Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of Mathematics and Mathematical Sciences
Mathematical Problems in Engineering
Journal of
Mathematics Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Discrete Mathematics
Journal of
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Discrete Dynamics in Nature and Society
Journal of
Function Spaces Hindawi Publishing Corporation http://www.hindawi.com
Abstract and Applied Analysis
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Journal of
Stochastic Analysis
Optimization
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014