1

Closed-Form Coverage Probability for Downlink
Poisson Network With Double Shadowed Fading

arXiv:1712.09469v1 [cs.IT] 27 Dec 2017

Jingrui Chen and Chaowei Yuan
Abstract—Performances of cellular networks over κ-µ shadowed fading with long-term shadowing has been studied in
the existing literature. However, the impact of κ-µ shadowed
fading with instantaneous shadowing on performances of cellular networks is unknown. Therefore, this letter analyzes the
downlink coverage probability of a Poisson network with double
shadowed fading which is composed of a large-scale fading of
lognormal distribution and κ-µ shadowed fading with integer
fading parameters. The closest base station association rule
without shadowing is considered. For analytical tractability,
the double shadowed fading is approximated as a weighted
sum of κ-µ shadowed distributions based on the GaussianHermit quadrature. As a main theoretical result, a closed-form
expression for the downlink coverage probability of a Poisson
network under double shadowed fading for the desired signal and
arbitrary fading for the interfering signals is successfully derived.
Numerical simulations reveal that the double shadowed fading
provides a pessimistic coverage on a Poisson network compared
with the long-term shadowing which is incorporated into cell
selection.
Index Terms—Stochastic geometry, coverage probability, double shadowed fading, Poisson network.

I. I NTRODUCTION
Stochastic geometry as a useful tool has employed to the
performance analysis of heterogeneous cellular networks (HetNets), where the irregular locations of users and base stations
(BSs) are modeled by Point processes [1], [2]. The singletier cellular network is developed in [3] and then extended to
a multi-tier HetNet [4]. To meet the needs of various fading
scenarios in future communications, a versatile fading model,
i.e. κ-µ shadowed fading is proposed in [5]. However, it is
intractable to the performance analysis of HetNets over this
fading model due to the existence of the confluent hypergeometric function (CHF). To overcome this problem, the moment
matching method with Gamma distribution is utilized in [6]
and the truncated series form of the CHF is used in [7].
In order to accommodate the requirements of various fading
scenarios in fifth-generation communications, a double shadowed fading model is first proposed in [8], where the double
shadowed fading is composed of a large-scale fading and κ-µ
shadowed fading. Similar to [9], [10], the large-scale fading
is incorporated into cell association, which can characterize
the impact of long-term shadowing on the cell selection. As
a result, the effect of κ-µ shadowed fading with long-term
shadowing on performances of cellular networks is studied.
In this paper, we distinguish from the previous related
works [8]-[10] by studying the impact of instantaneous shadowing on the coverage probability of a Poisson network.
Specifically, we consider a double shadowed fading which
Jingrui Chen and Chaowei Yuan are with the School of Information and
Communication Engineering, Beijing University of Posts and Telecommunications, Beijing, 100876, China, (e-mail: xiaojingyun604604@163.com;
yuancw2000@bupt.edu.cn). Last updated: December 29, 2017.

is composed of a lognormal distribution and κ-µ shadowed
fading with integer fading parameters [11]. The nearest BS
selection without shadowing is also considered.
There are three contributions in this paper. Firstly, an
accurate and tractable expression for the double shadowed
fading is obtained according to the Gaussian-Hermit quadrature (GHQ) [12]. Secondly, a closed-form expression for the
downlink coverage probability of a Poisson network is successfully derived, assuming that the desired signal experiences
double shadowed fading and the interfering signals experience
arbitrary fading. The resultant expression avoids complex
integral evaluations or high order differential calculations.
Meanwhile, it is generic due to the application of versatile
model of double shadowed fading in the desired link and
the assumption of an arbitrary fading distribution for all
interfering links. Most importantly, the theoretical analysis of
coverage probability with respect to the interfering fading is
avoided under the independent assumptions of fading models.
Thirdly, compared with the long-term shadowed environment,
numerical simulations show that the double shadowed fading
offers a smaller coverage to a Poisson network.
II. S YSTEM M ODEL
Consider a downlink Poisson network where the locations of
BSs are modeled by a homogeneous Poisson point process ΦB
with density λ. Assuming that users are uniformly distributed
in the network. The received power of a user at origin (termed
as a typical user) from a BS at x ∈ ΦB is modeled as
−α
P (x) = P Hx kxk , where P is the transmit power, α is the
path loss exponent, and Hx is the independent fading gain with
unit power. Consider the nearest BS association policy, thus
−α
the location of serving BS is x∗ = arg maxx∈ΦB P kxk .
According to [3], the probability density function
 (PDF) of
serving distance is fR (r) = 2πλr exp −πλr2 . Assuming
that the channel power of the desired signal Hx∗ follows a double shadowed distribution, i.e. Hx∗ = Hχ ·HS , where the independent random variables Hχ follows a lognormal distribution
with a standard deviation (SD) σ̃S (σ̃S = σS ln (10)/10)
and HS follows a κ-µ shadowed distribution with integer
fading parameters [11]. For notational simplicity, we denote
the distributions Hχ and HS as Hχ ∼ ln N (0, σ̃S ) and
HS ∼ S (κ, µ, m; hS ) respectively. According to [11, Eq. (12),
Eq. (13)], the PDF and complementary cumulative distribution
function (CCDF) of HS are


M
i −1
X
Ci hm
hS
S
fHS (κ, µ, m; hS ) =
exp −
, (1)
Γ (mi ) Ωi mi
Ωi
i=0
F̄HS (κ, µ, m; hS ) =

M m
i −1
X
X
i=0 j=0

j

Ci (hS /Ωi )
, (2)
Γ (j + 1) exp (hS /Ωi )

2

Lemma 1. Given two mutually independent random variables Hχ and HS , where Hχ ∼ ln N (0, σ̃S ) and HS ∼
S (κ, µ, m; hS ), the exact PDF of the double shadowed fading
Hx∗ = Hχ · HS is

fHx∗ (hx∗ ) =

i=0

0

S

Exa.κ=1,µ=3,m=2,σS=0dB
Exa.κ=1,µ=2,m=1,σS=8dB
Exa.κ=1,µ=3,m=2,σS=8dB

1

Approx. Results
Empirical Values

0.8

0.6

0.4

0.2

0
0

0.5

2σ̃

S

i
2πσ̃S Γ (mi ) Ci−1 Ωm
i hx∗

2

2.5

3

Fig. 1. PDF of double shadowed fading for nGHQ = 32.

Lemma 2. The double shadowed distribution described in
Lemma 1 can be expressed as
nGHQ

X

fHx∗ (hx∗ ) =

al fHS (κ, µ, m; hx∗ /bl ),

(5)

l=1
nGHQ
P

wl1 , bl = exp

√


2σ̃S tl1 , wl1 and tl1

are the weights and abscissas of the l1 -th order Hermite
polynomial respectively, and fHS (·) is given by (1).
Proof: The density function of product of two random
variables [13] Hx∗ = Hχ · HS can be rewritten as


Z ∞
1
hx∗
dhχ .
(6)
fHx∗ (hx∗ ) =
fHχ (hχ ) fHS
h
|h
χ
χ|
0
Substituting the distribution of Hχ ∼ ln N (0, σ̃S ) and (1)
into (6), we can obtain

m −1

i

1.5

hx∗

hS i

 dhS
(ln hx∗ −ln hS )2
hS
exp Ω +
2

√

1

l1 =1

In order to study the impact of κ-µ shadowed fading
with instantaneous shadowing on performances of cellular
networks, the exact PDF of the double shadowed distribution
described in section II is given bellow.

M
X

Exa.κ=1,µ=2,m=1,σ =0dB
1.2

where al = wl /

III. D OUBLE SHADOWED FADING

R∞

1.4

fHx∗ (hx∗ )

where the coefficients M , mi , Ωi and Ci are depended
m
on the fading parameters κ, µ and m. Let ω1 = κµ+m
,
κµ
1
ω2 = κµ+m and ω3 = µ(κ+1) . According to [11, Table
I], if µ > m, then M
 = µ, mi = µ − m − i + 1,
m+i−2
m
Ωi = ω3 and Ci =
(−ω1 ) ω21−m−i for
i−1
ω3
i = 1, 2, .. . , µ − m, and m
i = µ − i + 1, Ωi = ω1
i−2
i−µ+m−1 1−i
and Ci =
(−ω1 )
ω2 for i =
i−µ+m−1
µ − m + 1, . . . , µ, otherwise
M

 = m − µ, mi = m − i,
m−µ
ω3
ω1i ω2m−µ−i . Assuming that
Ωi = ω1 and Ci =
i
the channel powers of the interfering signals {Hx } are i.i.d.,
where x ∈ ΦB \x∗. In order to obtain a generic result, we
consider Hx follows an arbitrary distribution. The orthogonal
multiple access technology [3] is employed at each cell to
eliminate intra-cell interference. Thus, by considering the
interference-limited (noise is ignored) network, the signal-tointerference ratio (SIR) of a typical user can be written as
−α
SIR (r) = P Hx∗I r , where the aggregate interference power
P
−α
I = x∈ΦB \x∗ P Hx kxk .

,

(3)

where the coefficients M , mi , Ωi and Ci are given by (1).
Proof: Based on the density function of product of two
random variables [13], the PDF Hx∗ = Hχ · HS is


Z ∞
1
hx∗
fHx∗ (hx∗ ) =
fHS (hS ) fHχ
dhS .
(4)
h
|h
S
S|
0
Substituting the distribution of Hχ ∼ ln N (0, σ̃S ) and (1)
into (4) yields (3). This completes the proof.
It has been known that the κ-µ shadowed fading contains
majority of the fading models proposed in the literature as
special cases [8], including Rayleigh, Rician, Nakagami-m,
Nakagami-q, One-sided Gaussian, κ-µ, η-µ, etc. Therefore,
Lemma 1 encompasses various composite fading models and
shows the instantaneous property of large scale fading of
lognormal distribution. However, it is intractable to the performance analysis of cellular networks due to an integral
existed in (3). To conquer this problem, the GHQ rule [12,
Eq. (25.4.46)] is employed and the derived expression is as
follows.

fHx∗ (hx∗ ) =

M
X
i=0

R∞
0

√



(ln hχ )2
exp − Ωhx∗
2
hχ −
2σ̃
S

i

m +1

hχ i

dhχ

i 1−mi
2πσ̃S Γ (mi ) Ci−1 Ωm
i hx∗

.

(7)

ln h

Using the change of variable t = √2σ̃χ in (7), applying the
S
Pn
√
GHQ rule and combing with the identity π = l1GHQ
=1 wl1
yields (5). This completes the proof.
As illustrated in Fig. 1, the approximated expression provided in (5) exactly matches the exact expression in Lemma 1.
Most importantly, the fHx∗ (hx∗ ) in Lemma 2 is expressed
as a weighted sum of κ-µ shadowed PDFs fHS (·), which
simplifies the analytical expressions significantly in special
cases. When the lognormal distribution is ignored, (5) is
reduced to (1) as σS = 0dB. Another important expression is
the CCDF of Hx∗ , which is directly utilized for the coverage
analysis of a downlink Poisson network. Combing with (2)
and Lemma 2, we can obtain the CCDF of Hx∗ conveniently.
Corollary 1. The CCDF of the double shadowed fading is
nGHQ

F̄Hx∗ (hx∗ ) =

X

al F̄HS (κ, µ, m; hx∗ /bl ),

(8)

l=1

where the coefficients al and bl are given by Lemma 2 and
F̄HS (·) is given by (2).

3

α

IV. C OVERAGE P ROBABILITY
The coverage probability in a downlink Poisson network
is defined as the probability that the received power of a
typical user is bigger than the SIR threshold θ [3]. Consider
the nearest BS association policy, the closed-form expression
for the downlink coverage probability in a Poisson network is
obtained in the next Theorem.
Theorem 1. When the desired signal experiences the double
shadowed fading and the interfering signals experience arbitrary fading, the coverage probability of a typical user is
pc (θ) =

GHQ
M m
i −1 nX
X
X

al Ci Aj Γ (Bj )

X

j

i=0 j=0

Bj

,

(9)

(−1) Gi

l=1 (N1 ,··· ,Nj )∈Tj

where
the
coefficients
M,
mi ,
Ωi ,
Ci ,
al
and
b
are
given
by
Corollary
1,
T
l
j ) =
(
j
j
P
j P
Nq q = j,
Nq = q ,
(N1 , · · · , Nj ) ∈ (N ∪ {0})
q=1

Aj

j
Q

=

q=1

Bj = 1 +

j
P

1
Γ(Nq +1)

q=1



2(−1)q Ehq
αΓ(q+1)

Nq , Gi = 1 + (θ/Ωi bl )

2/α



θ
Ωi bl

2/α Nq

,

Eh0 , the values of the

q=1

expectations withhrespect to the interfering fading distributions
 i
R∞
are Eh0 = EHx (Ωi bl /θ)2/α 1 − exp −hx t−α/2 dt and
i
h

2/α
Ehq = EHx hx γ q − 2/α, Ωθhi bxl , and the functions Γ (·)
and γ (a, x) are the Gamma and lower incomplete Gamma
function [14] respectively.

where L = Ωiθrbl vα . According to the higher order derivatives [14, Eq. (04.30.2)], we obtain

R ∞
∂j
∂z j exp −2πλEHx Rr ∞(1 − exp (−zLhx )) vdv 
= j! exp −2πλEHx r (1 − exp (−zLhx )) vdv

j
P
Q
(14)
1
1
∂q
×
Γ(Nq +1) Γ(q+1) ∂z q (−2
(N1 ,··· ,Nj )∈Tj q=1
R ∞
×πλEHx r (1 −

Since the Hx∗ follows the double shadowed fading, combined
by (8) in Corollary 1, we obtain

j
α
GHQ
M m
i −1 nX
Ci al PθIr
X
X
Ωi bl

.
(11)
pc (θ, r) =
θIr α
i=0 j=0 l=1 Γ (j + 1) exp P Ωi bl
Due to the fact that exp (−α) = (−α)
where g
as

(j)

−j

exp (−αz)

(j)
z=1

,

(z) is the j-th differential of z, (11) can be written

pc (θ, r) =

GHQ
M m
i −1 nX
X
X

∂j
∂z j E

h


i
exp − ΩzsI
i bl

z=1

j

i=0 j=0

l=1

−1
a−1
l Ci (−1) Γ (j + 1)

, (12)

where s = θrα /P . Since Hx has to be identically distributed
as well as for all interfering BSs, substituting I =
P
−α
into (12) and using the result of the
x∈ΦB \{x∗} P Hx kxk
conditional generating function of PPP [3], we get
M mP
i −1 nGHQ
P
P Ci al (−1)j ∂ j
pc (θ, r) =
Γ(j+1) ∂z j exp (−2πλ
R ∞i=0 j=0 l=1

× EHx r (1 − exp (−zLhx )) vdv z=1 ,

(13)

Dq

,

where Tj is given
R ∞ in Theorem 1. Solving the derivative
∂q
−2πλE
(1 − exp (−zLhx )) vdv , plugging in
Hx
∂z q
r
z = 1 and L, and by using the change of variable t =
−1/α −1
(θ/Ωi bl )
r v, we have

R ∞
q
∂
(1 − exp (−zLhx )) vdv
∂z q −2πλEHx 
r
2/α
(15)
2πλr 2 (−1)q Ehx q
θ
=
,
α
Ωi bl
h
i

2/α
where Ehq
=
EHx hx γ q − 2/α, Ωθhi bxl
.
In R
addition, substituting z  =
1 and L into
∞
EHx r (1 − exp (−zLhx )) vdv and by using the change
α
of variable t = Ωθri blhvxα , we obtain
R ∞

exp −2πλE
(1 − exp
(−zLhx )) vdv
Hx
r

2/α

(16)
Eh0 ,
= exp −πλr2 Ωθi bl
hR
 i
∞
where Eh0 = EHx (Ωi bl /θ)2/α 1 − exp −hx t−α/2 dt .
Substituting (14), (15) and (16) into (13), we have
pc (θ, r) =

M mP
i −1 nGHQ
P
P
i=0 j=0

Proof: Substituting the expression of SIR into the definition of the coverage probability, we have


θI1 rα
pc (θ, r) = P hx∗ >
.
(10)
P

exp (−zLhx )) vdv

Wi,j (t) =

P

Ci al (−1)j Wi,j (r 2 )

,

exp(πλr 2 (θ/Ωi bl )2/α Eh0 )
l=1

 2 N q
2πλtEh 
q
α
θ
j
Q
(−1)q αq! Ω b
i l

(N1 ,··· ,Nj )∈Tj q=1

Nq !

(17)
.

Note that if j = 0, then Wi,j (t) = 1. Substituting the
PDF fR (r) = 2πλr exp −πλr2 and (17) into the coverage
probability pc (θ) = E [pc (θ, r)] and by some algebraic
manipulations, we get
Z ∞
GHQ
M m
i −1 nX
X
X
Ci al
Wi,j (r/2πλ)
pc (θ) =
dr, (18)
j
exp (rGi )
0
i=0 j=0 l=1 (−1)
2/α

where Gi = 1 + (θ/Ωi bl ) PEh0 . Because of the fact that
j
P
Wi,j (t) =
Aj t q=1 Nq , where Aj is given in
(N1 ,··· ,Nj )∈Tj

Theorem 1, the integral term of (18) can be rewritten as
R∞
exp (−rGi ) Wi,j (r/2πλ) dr
0
Pj
R∞
P
(19)
=
Aj 0 exp (−rGi ) r q=1 Nq dr.
(N1 ,··· ,Nj )∈Tj

Finally, by evaluating the integral in (19) based on [14,
Eq. (3.35.3)] and with some algebraic manipulations in (18)
yields (9). This completes the proof.
Obviously, (9) avoids complicated integrals or high order
differential calculations. As the popular fading distributions
such as small-scale fading models, line-of-sight shadowing
models and traditional composite fading/shadowing models are
special cases of double shadowed fading [8, Fig. 1(b)], the
closed-form expression of the coverage probability obtained

4

1

1

0.9

0.9

0.7

0.95

0.8

0.9

0.7

Coverage Probability

Coverage Probability

0.8

0.85
0.6

0.8

0.5

Ana. (DS:κ=1,µ=2,m=1,σS=8)

0.4

Ana. (DS:κ=1,µ=2,m=1,σS=0)

0.75
−15

−10

−5

Ana. (DS:κ=1,µ=3,m=2,σS=8)
0.3

Ana. (DS:κ=1,µ=3,m=2,σS=0)

0.2

Sim. DS & σS=8dB

0.1

Sim. (LT:κ=1,µ=2,m=1,σS=8)
Sim. (LT:κ=1,µ=3,m=2,σS=8)

0
−30

−20

−10

0.6
0.5

Ana. (DS:κ=1,µ=1,m=1,σS=8)

0.4

Ana. (DS:κ=1,µ=1,m=1,σS=0)

0.3

Ana. (DS:κ=1,µ=3,m=3,σS=8)
Ana. (DS:κ=1,µ=3,m=3,σS=0)

0.2

Sim. DS & σS=8dB

0.1

Sim. (LT:κ=1,µ=1,m=1,σS=8)
Sim. (LT:κ=1,µ=3,m=3,σS=8)

0

SIR Threshold in dB

10

20

30

0
−30

−20

−10

0

SIR Threshold in dB

10

20

30

Fig. 2. Coverage probabilities with DS and LT for M = µ.

Fig. 3. Coverage probabilities with DS and LT for M = m − µ.

is generic. When the instantaneous shadowing of lognormal
distribution is ignored, i.e. σS = 0dB, Theorem 1 is reduced to
the κ-µ shadowed fading scenario experienced by the desired
signal. Most importantly, the theoretical analysis of coverage
probability concerning the interfering fading is avoided and it
is replaced by simple calculations Eh0 and Ehq .

distribution and κ-µ shadowed fading with integer fading
parameters are provided. Based on the resulting expression, the
generic and closed-form expression for the downlink coverage
probability is successfully derived and verified through simulation. The obtained expression of the coverage probability
indicates that the theoretical analysis of coverage probability
about the interfering fading is avoided. In addition, compared
to the long-term shadowed environment, numerical simulations
show that the instantaneous shadowing of large scale fading
is bad for the coverage of a Poisson network.

V. N UMERICAL RESULTS
In this section, we validate the accuracy of derived coverage
probability through Monte Carlo and make a comparison
between the double shadowed fading and long-term shadowed scenario. For notational simplicity, we use DS and LT
to stand for the double shadowed fading and comparative
term of long-term fading respectively, which are experienced
by the desired signal. The interfering signals experience
Rayleigh/Lognormal fading and Rayleigh fading separately.
The simulation parameters of the Poisson network are λ =
10−7 , α = 4 and P = 1. For numeral calculations, we set
nGHQ = 32. In Fig. 2, the coverage probability is plotted
for the DS model when M = µ (i.e. µ > m). It can be
observed that the analytically obtained results closely match
the simulation results. Intuitively, the coverage probabilities
increase as parameters of m and µ increase and decrease as
the values of σS increase. The coverage probability of the LT
model is also considered in Fig. 2. Comparing the coverage
probabilities of the DS model and LT model, we can see that
the instantaneous property of large scale fading makes a great
difference in coverage. The DS model provides a pessimistic
coverage. This is due to the fact that choosing the serving
BS depends solely on the distance. Hence, the instantaneous
property of large scale fading is presented. When M = m − µ
(i.e. µ ≤ m), the coverage probabilities of the DS model and
LT model are shown in Fig. 3 and exact matches can be seen
for various fading parameters. Similarly, it reveals the same
conclusions as the Fig. 2.
VI. C ONCLUSION
This paper analyzes the downlink coverage probability
for a Poisson network when the desired signal experiences
double shadowed fading and the interfering signals experience
arbitrary fading. The nearest BS association policy without
shadowing is considered. The exact and approximate expressions for the double shadowed fading composed of a lognormal

R EFERENCES
[1] H. Elsawy, E. Hossain, and M. Haenggi, “Stochastic geometry for
modeling, analysis, and design of multi-tier and cognitive cellular
wireless networks: A survey,” IEEE Commun. Surveys & Tuts., vol. 13,
no. 3, pp. 996–1019, 3rd Quart. 2013.
[2] J. G. Andrews, A. K. Gupta, and H. S. Dhillon, “A primer on cellular
network analysis using stochastic geometry,” arXiv:1604.03183, Oct.
2016.
[3] J. G. Andrews, F. Baccelli, and R. Ganti, “A tractable approach to
coverage and rate in cellular networks,” IEEE Trans. Commun., vol. 59,
no. 11, pp. 3122–3134, Nov. 2011.
[4] H. S. Dhillon, R. K. Ganti, F. Baccelli, and J. G. Andrews, “Modeling
and analysis of K-tier downlink heterogeneous cellular networks,” IEEE
J. Sel. Areas Commun., vol. 30, no. 3, pp. 550–560, Mar. 2012.
[5] J. F. Paris, “Statistical characterization of κ-µ shadowed fading,” IEEE
Trans. Veh. Technol., vol. 63, no. 2, pp. 518–526, Feb. 2014.
[6] S. Kumar, “Approximate out probability and capacity for κ-µ shadowed
fading,” IEEE Wireless Commun. Lett., vol. 4, no. 3, pp. 10–13, Jun.
2015.
[7] S. Parthasarathy and R. K. Ganti, ”Coverage analysis in downlink
Poisson cellular network with κ-µ shadowed fading,” IEEE Wireless
Commun. Lett., vol. 6, no. 1, pp. 10–13, Feb. 2017.
[8] Y. J. Chun, S. L. Cotton, H. S. Dhillon, F. J. Lopez-Martinez, J. F.
Paris, and S. K. Yoo, ”A comprehensive analysis of 5G heterogeneous
cellular systems operating over κ-µ shadowed fading channels,” IEEE
Trans. Wireless Commun., vol. 16, no. 11, pp. 6995–7010, Nov. 2017.
[9] H. S. Dhillon and J. G. Andrews, ”Downlink rate distribution in
heterogeneous cellular networks under generalized cell selection,” IEEE
Wireless Commun. Lett., vol. 3, no. 1, pp. 42–45, Feb. 2014.
[10] C. Saha, M. Afshang, and H. S. Dhillon, “Enriched K-tier hetnet model
to enable the analysis of user-centric small cell deployments,” IEEE
Trans. Wireless Commun., vol. 16, no. 3, pp. 1953–1608, Mar. 2017.
[11] F. J. Lopez-Martinez, J. F. Paris, and J. M. Romero-Jerez, “The κ-µ
shadowed fading model with integer fading parameters,” IEEE Trans.
Veh. Technol., vol. 66, no. 9, pp. 7653–7662, Sept. 2017.
[12] M. Abramowitz and I. A. Stegun, Handbook of mathematical functions
with formulas, graphs, and mathematical tables, 9nd ed. New York,
USA: Dover Publications, 1972.
[13] A. Papoulis, Probability, Random Variables, and Stochastic Processes.
New York, NY, USA: McGraw-Hill, 1991.
[14] I. S. Gradshteyn and I. M. Ryzhik, Tables of integrals, series, products,
7nd ed. San Diego, CA, USA: Academic, 2007.

