Comparison between two common collocation approaches based on
radial basis functions for the case of heat transfer equations arising
in porous medium
K. Parand1,a , S. Abbasbandyb,∗, S. Kazemb , A.R. Rezaeia
a

arXiv:1008.2348v2 [math-ph] 16 Aug 2010

b

Department of Computer Sciences, Shahid Beheshti University, G.C., Tehran, Iran
Department of Mathematics, Imam Khomeini International University, Ghazvin 34149-16818, Iran

Abstract
In this paper two common collocation approaches based on radial basis functions have been
considered; one be computed through the integration process (IRBF) and one be computed
through the differentiation process (DRBF). We investigated the two approaches on natural
convection heat transfer equations embedded in porous medium which are of great importance in the design of canisters for nuclear wastes disposal. Numerical results show that the
IRBF be performed much better than the common DRBF, and show good accuracy and
high rate of convergence of IRBF process.
Keywords: Collocation method; Nonlinear ODE; Radial Basis Functions; Direct Inverse
Multiquadric; Indirect Multiquadric; Porous media.
PACS: 47.56.+r, 02.70.Hm

1. Introduction
Natural convective heat transfer in porous media has received considerable attention
during the past few decades. This interest can be attributed due to its wide range of
applications in ceramic processing, nuclear reactor cooling system, crude oil drilling, chemical
reactor design, ground water pollution and filtration processes. External natural convection
in a porous medium adjacent to heated bodies was analyzed by Nield and Bejan [1], Merkin
[2, 3], Minkowycz and Cheng [4, 5], Pop and Cheng [6, 7], Ingham and Pop [8]. In all of
these analysis, it was assumed that boundary layer approximations are applicable and the
coupled set of governing equations were solved by numerical methods.
In this paper, the same approximations are applied to the problem of natural convection
about an inverted heated cone embedded in a porous medium of infinite extent. No similarity
∗

Corresponding author. Tel:(+98912) 1305326 Fax:(+98281) 3780040
Email addresses: k_parand@sbu.ac.ir (K. Parand), abbasbandy@yahoo.com (S. Abbasbandy),
saeedkazem@gmail.com (S. Kazem), alireza.rz@gmail.com (A.R. Rezaei)
1
Member of research group of Scientific Computing.
August 17, 2010

solution exists for the truncated cone, but for the case of full cone, if the prescribed wall
temperature or surface heat flux is a power function of distance from the vertex of the
inverted cone similarity solutions exist [1, 6], a great deal of information is available on heat
and fluid flow about such cones as reviewed by Refs. [9, 10].
Bejan and Khair [11] used Darcy’s law to study the vertical natural convective flows
driven by temperature and concentration gradients. Nakayama and Hossain [12] applied
the integral method to obtain the heat and mass transfer by free convection from a vertical
surface with constant wall temperature and concentration. Yih [13] examined the coupled
heat and mass transfer by free convection over a truncated cone in porous media for variable
wall temperature and concentration or variable heat and mass fluxes and [14] applied the
uniform transpiration effect on coupled heat and mass transfer in mixed convection about
inclined surfaces in porous media for the entire regime. Cheng [15] used an integral approach
to study the heat and mass transfer by natural convection from truncated cones in porous
media with variable wall temperature and [16] studied the Soret and Dufour effects on the
boundary layer flow due to natural convection heat and mass transfer over a vertical cone in
a porous medium saturated with Newtonian fluids with constant wall temperature. Natural
convective mass transfer from upward-pointing vertical cones, embedded in saturated porous
media, has been studied using the limiting diffusion [17]. The natural convection along
an isothermal wavy cone embedded in a fluid-saturated porous medium are presented in
[18, 19]. Lai and Kulacki [20] studied the natural convection boundary layer flow along a
vertical surface with constant heat and mass flux including the effect of wall injection. In
[21] fluid flow and heat transfer of vertical full cone embedded in porous media have been
solved by Homotopy analysis method [22, 23].
Mathematical modeling of many problems in science and engineering leads to ordinary
differential equations (ODEs) [24–28]. The methods based on radial basis functions (RBF)
which are part of an emerging field of mathematics are famous ways to solve these kinds
of problems . First studied by Roland Hardy, an Iowa State geodesist, in 1968, these
methods allow for scattered data to be easily used in computations [29]. The concept of
solving DEs by using RBFs was first introduced by Kansa [30]. Since then, it has received
a great deal of attention from researchers. And consequently, many further interesting
developments and applications have been reported (e.g. Zerroukat et al.[31], Mai-Duy and
Tran-Cong[32, 33]). Essentially, in a typical RBF collocation method, each variable and its
derivatives are all expressed as weighted linear combinations of basis functions, where the
sets of network weights are identical. These closed forms of representations are substituted
with the governing equations as well as boundary conditions, and the point collocation
technique is then employed to discretize the system [27]. If all basis functions in networks
are available in analytic forms, the RBF collocation methods can be regarded as truly
meshless methods [34]. There are two basic approaches for obtaining new basis functions
from RBFs, namely direct approach (DRBF) based on a differential process (Kansa [30])
and indirect approach (IRBF) based on an integration process (Mai-Duy and Tran-Cong
[29, 32, 35]). Both approaches were tested on the solution of second order DEs and the
indirect approach was found to be superior to the direct approach (Mai-Duy and Tran-Cong
[32]).
2

In this paper we apply the DRBF and IRBF for solving natural convection of Darcian
fluid about a vertical full cone embedded in porous media prescribed surface heat flux which
is third order nonlinear ODE.
2. Problem formulation
Consider an inverted cone with semi-angle γ and take axes in the manner indicated in
Fig. 1(a). The boundary layer develops over the heated frustum x = x0 .
The boundary layer equations for natural convection of Darcian fluid about a cone are
[6]:
∂
∂
(ru) +
(rv) = 0,
∂x
∂y
ρ∞ βKg cos γ(T − T∞ )
u=
,
µ
∂T
∂2T
∂T
+v
=α 2.
u
∂x
∂y
∂y

(2.1)
(2.2)

For a thin boundary layer, r is obtained approximately xsin(γ). Suppose that a power law
of heat flux is prescribed on the frustum. Accordingly, the boundary conditions at infinity
are:
u = 0,

T = T∞ , if y → ∞

(2.3)
(2.4)

and at the wall are
v = 0 if y = 0.
If the surface heat flux qw [6] is prescribed, qw is obtained as
qw = −k(

∂T
)y=0 = A(x − x0 )λ ,
∂y

x0 ≤ x ≤ ∞.

For the case of a full cone (x0 = 0, F ig.1(b)) a similarity solution exists [6].
In the case of prescribed surface heat flux the similarity solution for the stream function
ψ and T where
u=

1 ∂ψ
,
r ∂y

v=−

1 ∂ψ
r ∂x

(2.5)

is of the form [6]:
ψ = αr(Rax )1/3 f (η),
1
qw x
T − T∞ =
(Rax )− 3 θ(η),
k
y
1/3
η = (Rax ) ,
x

(2.6)

3

where
Rax =

ρ∞ βgK cos(γ)qw x2
µαk

(2.7)

is the local Rayleigh number for the case of prescribed surface heat flux. The governing
equations become
f ′ = θ,
λ + 5 ′ 2λ + 1 ′
fθ −
f θ = 0,
θ′′ +
2
3
subjected to boundary conditions as:
f (0) = 0,

θ′ (0) = −1,

θ(∞) = 0.

Finally from Eqs. (2.8) and (2.9) we have:
(
 ′′
 ′ 2
f f − 2λ+1
(f ) = 0,
ODE. f ′′′ + λ+5
2
3
′′
B.C. f (0) = 0, f (0) = −1, f ′ (∞) = 0.

(2.8)

(2.9)

(2.10)

It is of interest to obtain the value of the local Nusselt number which is defined as [6]:
qw x
Nux =
.
(2.11)
k(Tw − T∞ )

From Eqs. (2.11), (2.6) and (2.7) it follows that the local Nusselt number which is interest
to obtain given by:
Nux = Ra1/3
x [−θ(0)].

(2.12)

3. RBF Functions
Let R+ = {x ∈ R, x ≥ 0} be the non-negative half-line and let φ : R+ → R be a
continuous function with φ(0) ≥ 0. A radial basis function on Rd is a function of the form
φ(kX − Xi k)

(3.1)

where X, Xi ∈ Rd and k.k denotes the Euclidean distance between X, Xi . If one chooses
N points {Xi }N
i=1 in R then by custom
s(X) =

N
X
i=1

λi φ(kX − Xi k);

λi ∈ R,

(3.2)

is called a radial basis function as well [36].
In order to explain RBF methods briefly, suppose that the one-dimensional input data point
set or the center set xi in the given domain Ω ⊆ R is given. The center point is not necessarily
structured, that is, it can have an arbitrary distribution. The arbitrary grid structure is one
of the major differences between the RBF method and other global methods. Such a meshfree grid structure yields high flexibility especially when the domain is irregular. In this
work the uniform grid is used for RBF approximation.
4

3.1. Properties of RBF
With a radial function φ(r) and with data values ui given at the locations xi , for i =
1, 2, ..., N the function

s(x) =

N
X

λi φi (x)

(3.3)

i=1

where r = ri = kx − xi k and φi (x) = φ(kx − xi k), interpolates the data if we choose the
expansion coefficients λi in such a way that s(xj ) = uj , for j = 1, 2, ..., N [37, 38]. The
expansion coefficients λi can therefore be obtained by solving the linear system AΛ = U,
where:
Aij = φ(kxj − xi k),
Λ = [λ1 , λ2 , ..., λN ]T ,
U = [u1 , u2 , . . . , uN ]T .

(3.4)
(3.5)
(3.6)

All the infinitely smooth RBF choices listed in Table (1) will give coefficient matrices
A in (3.4) which are symmetric and nonsingular [39], i.e. there is a unique interpolant of
the form (3.3) no matter how the distinct data points are scattered in any number of space
dimensions. In the cases of inverse quadratic, inverse multiquadric and GA the matrix
A is positive definite and, for multiquadric (MQ), it has one positive eigenvalue and the
remaining ones are all negative [39].
Interpolation using Conical splines and thin-plate splines (TPSs) can become singular in
multidimensions [34]. However, low-degree polynomials can be added to the RBF interpolant
to guarantee that the interpolation matrix is positive definite (a stronger condition than
nonsingularity). For example, for the Conical RBF and
Pm the TPS inPdNdimensions this
becomes the case if we use as an interpolant s(x) =
i=1 λi φ(kx − xi k)
i=1 ai pi (x) +
PN
together with the constraints j=1 λj pi (xj ) = 0, for i = 1, 2, ..., m. Here pi (x) denotes a
basis for polynomials of Pqd in Rd (Pqd denotes the space of d-variate polynomials of order
not exceeding q) and m = (q − 1 + d)!/(d!(q − 1)!) [40].
3.2. RBF Interpolation
One dimensional function u(x) to be interpolated or approximated can be represented
by an RBF as:
N
X

λi φi (x) = ΦT (x)Λ

(3.7)

ΦT (x) = [φ1 (x), φ2 (x), ..., φN (x)],
Λ = [λ1 , λ2 , ..., λN ]T ,

(3.8)
(3.9)

u(x) ≈ s(x) =

i=1

where

5

x is the input and λi s are the set of coefficients to be determined. By choosing N interpolate
nodes (xj , j = 1, 2, ..., N) in Ω ∪ ∂Ω, the function u(x) can be approximated in Ω ∪ ∂Ω.
uj =

N
X

λi φi (xj ),

(j = 1, 2, ..., N).

(3.10)

i=1

To brief discussion on coefficient matrix we define:
AΛ = U

(3.11)

where
U =[u1 , u2 , ..., uN ]T ,

(3.12)

A =[ΦT (x1 ), ΦT (x2 ), ..., ΦT (xN )]T

φ1 (x1 ) φ2 (x1 ) . . . φN (x1 )
 φ1 (x2 ) φ2 (x2 ) . . . φN (x2 )

=
..
..
..
..

.
.
.
.
φ1 (xN ) φ2 (xN ) . . . φN (xN )





.


(3.13)

Note that φi (xj ) = φ(kxi − xj k) therefore φi (xj ) = φj (xi ) consequently A = AT .
The shape parameter c which is shown in Table (1) affects both the accuracy of the approximation and the conditioning of the interpolation matrix [44]. In general, for a fixed number
of N, smaller shape parameters produce the more accurate approximations, but also are
associated with a poorly conditioned A. The condition number also grows with N for fixed
values of the shape parameter c. Many researchers (e.g.[41, 42]) have attempted to develop
algorithms for selecting optimal values of the shape parameter. The optimal choice of the
shape parameter is still an open question. In practice it is most often selected by brute force.
Recently, Fornberg et. al.[43] have developed a Contour-Padé algorithm which is capable of
stably computing the RBF approximation for all c > 0 [44].
The following theorem about the convergence of RBF interpolation is discussed [45, 46].
Theorem 3.1. assume xi , (i = 1, 2, ..., N) are N nodes in Ω which is convex, let
h = max min kx − xi k2
x∈Ω 1≤i≤N

R
when φ̂(η) < c(1 + |η|)−(2l+d) for any u(x) satisfies (û(η))2 /φ̂(η)dη < ∞ we have
kuN (α) − u(α) k∞ ≤ chl−α

where φ(x) is RBF and the constant c depends on the RBF, d is space dimension, l and
α are nonnegative integer. It can be seen that not only RBF itself but also its any order
derivative has a good convergence.
6

3.3. Direct RBF for ODEs (DRBF)
In the direct method, the closed form DRBF approximating function (3.7) is first obtained from a set of training points, and its derivative of any order, e.g. pth order, can then
be calculated in a straightforward manner by differentiating such a closed form DRBF as
follows:
N
N
N
X
dk s(x)
dk X
dk φi (x) X
[k]
= k(
λi φi (x)) =
λi
=
λi Gi (x)
k
k
dx
dx i=1
dx
i=1
i=1

(3.14)

where
[k]

Gi (x) = dk φi (x)/dxk ,

k = 0, 1, ..., p.

Now we aim to apply the DRBF method for solving the ODEs in general form :
(
F (x, u, u′, ..., u(p−1) , u(p) ) = 0, a ≤ x ≤ b,
u(i) (ei ) = αi+1 ,
i = 0, 1, ..., p − 1,

(3.15)

where ei ∈ {a, b} and u(i) (x) = di u(x)/dxi , F is known function and {αi }pi=1 are known
constants. By substituting Eq. (3.7) in (3.15) and using Eq. (3.14) we have:
F (x, s, s′ , ..., s(p) ) = F (x,

N
X
i=1

[0]

λi Gi (x),

N
X
i=1

[1]

λiGi (x), ...,

N
X

[p]

λi Gi (x)).

i=1

Now, to obtain λi s (i = 1, 2, ..., N) we define the residual function:
Res(x) = F (x, s, s′ , ..., s(p) ).

(3.16)

The set of equations for obtaining the coefficients {λi }N
i=1 come from equalizing Eq. (3.16)
N −p
to zero at N − p interpolate nodes {xj }j=1 plus p boundary conditions:
(
Res(xj ) = 0,
j = 1, 2, ..., N − p,
(3.17)
PN
[k]
i=1 λi Gi (ei ) = αi+1 , k = 0, 1, ..., p − 1.
Since the direct approach is based on a differentiation process, all derivatives obtained here
are very sensitive to noise arising from the interpolation of DRBFs from a set of discrete
data points. Any noise here, even at the small level, will be badly magnified with an increase
in the order of derivative [29].

7

3.4. Indirect RBF for ODEs (IRBF)
In the indirect method, the formulation of the problem starts with the decomposition
of the highest order derivative under consideration into RBFs. The obtained derivative
expression is then integrated to yield expressions for lower order derivatives and finally
for the original function itself. In contrast, the integration process, where each integral
represents the area under the corresponding curve, is much less sensitive to noise. Based
on this observation, it is expected that through the integration process, the approximating
functions are much smoother and therefore have higher approximation power. Also To
numerically explore tile IRBF methods with shape parameters for which the interpolation
matrix is too poorly conditioned to use standard methods [47]. Let p be the highest order
of the derivative under consideration the boundary value ODEs in general form Eq. (3.15)
when (∃k s.t. ek = b) then we can define:
N

dp ŝ(x) X
=
λi φi (x),
dxp
i=1
Z X
Z
N
N
N
p−1
X
X
d ŝ(x)
[p−1]
=
λ
φ
(x)dx
=
λ
φ
(x)dx
=
λi hi (x) + d1 ,
i i
i
i
dxp−1
i=1
i=1
i=1

(3.18)

..
.

dŝ(x)
=
dx
=

Z X
N

[2]

λi hi (x)dx + d1

i=1
N
X

[1]
λi hi (x)

i=1

xp−3
xp−2
+ d2
+ ... + dp−1
(p − 2)!
(p − 3)!

xp−2
xp−3
+ d1
+ d2
+ ... + dp−1 ,
(p − 2)!
(p − 3)!

Z X
N
xp−3
xp−2
[1]
+ d2
+ ... + dp−1)dx
λi hi (x) + d1
ŝ(x) = (
(p − 2)!
(p − 3)!
i=1
=

N
X
i=1

[0]

λi hi (x) + d1

xp−2
xp−1
+ d2
+ ... + dp−1 x + dp ,
(p − 1)!
(p − 2)!

where
(R
φi (x)dx,
k = p − 1,
[k]
hi (x) = R [k+1]
hi (x)dx, k = 0, 1, ..., p − 2.

(3.19)

Substituting Eqs. (3.18) in (3.16) at N interpolate nodes {xj }N
j=1 plus p boundary conditions
p
the set of coefficients {λi }N
and
{d
}
is
obtained
as
follow
:
j j=1
i=1
(
Res(xj ) = 0,
ŝ(i) (ei ) = αi+1 ,

j = 1, 2, ..., N,
i = 1, 2, ..., p.
8

4. Solving the model
Consider governing equation of fluid flow and heat transfer of full cone embedded in
porous medium that is expressed by Eq. (2.10) for prescribed surface heat flux.
In the first step of our analysis, we approximate f (η) for solving the model by DRBF:
f (η) ≃ s(η) =

N
X

λi φi (η),

(4.1)

i=1

and f ′′′ (η) for solving the model by IRBF:
N

X
d3
f (η) ≃ 3 ŝ(η) =
λi φi (η).
dη
i=1
′′′

The general form of problem appear to:
(
Fλ (η, f, f ′ , f ′′ , f ′′′ ) = 0,
f (0) = 0, f ′′ (0) = −1, f ′ (∞) = 0.

(4.2)

(4.3)

To solve this problem we define residual function:
Res(η) = Fλ (η, s, s′, s′′ , s′′′ );

f or DRBF,

(4.4)

Res(η) = Fλ (η, ŝ, ŝ′ , ŝ′′ , ŝ′′′ );

f or IRBF.

(4.5)

The unknown coefficients {λi }N
i=1 come from equalizing Res(η) to zero at N interpolate
nodes ηi from Uniform distribution between 0 and η∞ which we set 9/2 for this problem.
4.1. Solving the model by DRBF
In the first step of solving, φi (η) is set by inverse multiquadric function which is shown in
Table (1). Now, the residual function is constructed by substituting Eq. (4.1) in Eq. (4.4):
Res(η) =

N
X
i=1

N
N
N
X
2λ + 1 X
λ+5 X
′′
′
)(
λi φi (η))(
λi φi (η)) − (
)(
λi φi (η))2 .
λi φi (η) + (
2
3
i=1
i=1
i=1
′′′

N −2
By using N − 2 interpolate nodes {ηj }j=1
plus two boundary conditions of Eq. (4.3)
′′
(f (0) = 0, f (0) = −1), the set of equations can be solved, consequently the coefficients
{λi }N
i=1 will be obtained:


Res(ηj ) = 0,
j = 1, 2, ..., N − 2,
P
N
i=1 λi φi (0) = 0,

PN
′′
i=1 λi φi (0) = −1.

9

Take into account φ′i (∞) = 0, for i = 1, 2, ..., N the infinity boundary condition (f ′ (∞) = 0)
is already satisfied.
Table (2) show the f ′ (η) for some λ in comparison with solutions of [21]. Also f ′ (η) for two
selected λ = 1/4 and 3/4 are showed in Table (3) in comparison with Runge-Kutta solution
is obtained by the MATLAB software command ODE45 which is used and applied by the
authors in ref. [21]. Absolute errors show that DRBF give us approximate solution with
a high degree of accuracy with a small N. The resulting graph of Eq. (2.10) is shown in
Figure (2).
4.2. Solving the model by IRBF
In the first of solving φi (η) is set by multiquadric function which is shown in Table (1).
Now, the residual function is constructed by substituting Eq. (4.2) in Eq. (4.5) and using
Eq. (3.18):
Z Z Z
Z
N
N
N
X
X
λ+5 X
)
λi
φi (η)dη
λi φi (η)dη
Res(η) =
λi φi (η) + (
2
i=1
i=1
i=1
Z
Z
N
2λ + 1 X
−(
)(
λi
φi (η)dη)2
3
i=1
=

N
X
i=1

N
N
X
d1 η 2
λ+5 X
[0]
[2]
)(
λi hi (η) +
+ d2 η + d3 )(
λi hi (η) + d1 )
λi φi (η) + (
2
2
i=1
i=1

N
2λ + 1 X
[1]
)(
λi hi (η) + d1 η + d2 )2 .
−(
3
i=1

By using N interpolate nodes {ηj }N
j=1 plus three boundary conditions of Eq. (4.3) the set of
3
equations can be solved, consequently the coefficients {λi }N
i=1 and {di }i=1 will be obtained.
In this method we put η∞ instead of infinity condition:

Res(ηj ) = 0,
j = 1, 2, ..., N,



PN λ h[0] (0) + d = 0,
i i
3
Pi=1
[2]
N


i=1 λi hi (0) + d1 = −1,
PN

[1]
i=1 λi hi (η∞ )d1 η∞ + d2 = 0.

Table (4) shows the f ′ (η) for some λ in comparison with solutions of [21]. Also f ′ (η)
for two selected λ = 1/4 and 3/4 are showed in Table (5) in comparison with Runge-Kutta
solution is obtained by the MATLAB software command ODE45 which is used and applied
by the authors in ref. [21]. Absolute errors show that IRBF give us approximate solution
with a high degree of accuracy with a small N. The resulting graph of Eq. (2.10) is shown
in Figure (3). A graph in figures (4) for λ = 2/3 show kResk2 for some N. Table (6) shows
kResk2 for some N and λ.
Comparison between DRBF solution in Table (3) with N = 12 and IRBF solution in Table
(5) with N = 10 for f ′ (η) show that the convergence of the IRBF method is faster, because
of using less numbers of collocation points.
10

5. Conclusion
In this paper we made a comparison between the two common collocation approaches
based on radial basis functions namely DRBF and IRBF methods on natural convection
equation about an inverted heated cone embedded in a porous medium of infinite extent
which are of great importance in the design of canisters for nuclear wastes disposal. These
functions are proposed to provide an effective but simple way to improve the convergence
of the solution by collocation method. The direct approach (DRBF) is based on a differentiation process, all derivatives are very sensitive to noise arising from the interpolation of
DRBFs from a set of discrete data points. Any noise, even at the small level, will be bad
magnified with an increase in the order of derivative. The indirect technique (IRBF) which
is based on integration process, each integral represents the area under the corresponding
curve, is much less sensitive to noise. Based on this observation, it is expected that through
the integration process, the approximating functions are much smoother and therefore have
higher approximation power. Additionally, through the comparison with other methods
such as HAM we show that the RBFs methods have good reliability and efficiency. Also
high convergence rates and good accuracy are obtained with the proposed method using
relatively low numbers of data points.
References
[1] D.A. Nield, A. Bejan, Convection in Porous Media, third ed., Springer-Verlag, New York, (2006).
[2] J.H. Merkin, Free convection boundary layers in a saturated porous medium with lateral mass flux, Int.
J. Heat Mass Transfer, 21 (1978) 1499-1504.
[3] J.H. Merkin, Free convection boundary layers on axisymmetric and two dimensional bodies of arbitrary
shape in a saturated porous medium, Int. J. Heat Mass Transfer, 22 (1979) 1461-1462.
[4] W.J. Minkowycz, P. Cheng, Local non-similar solutions for free convective flow with uniform lateral
mass flux in a porous medium, Lett. Heat Mass Transfer, 9 (1982) 159-168.
[5] W.J. Minkowycz, P. Cheng, F. Moalem, The effect of surface mass transfer on buoyancy induced Darcian
flow adjacent to a horizontal heated surface, Int. Commun. Heat Mass Transfer, 12 (1985) 55-65.
[6] P. Cheng, T.T. Le, I. Pop, Natural convection of a Darcian fluid about a cone, Int. Commun. Heat Mass
Transfer, 12 (1985) 705-717.
[7] I. Pop, P. Cheng, An integral solution for free convection of a Darcian fluid about a cone with curvature
effects, Int. Commun. Heat Mass Transfer, 13 (1986) 433-438.
[8] D.B. Ingham, I. Pop, Natural convection about a heated horizontal cylinder in a porous medium, J.
Fluid Mech. 184 (1987) 157-181.
[9] I. Pop, D.B. Ingham Convective heat transfer: mathematical and computational modeling of viscous
fluids and porous media, Pergamon Press, Oxford, 2001.
[10] K. Vafai, Handbook of porous media. Marcel Dekker, New York, 2000.
[11] A. Bejan, K.R. Khair, Heat and mass transfer by natural convection in a porous medium, Int. J. Heat
Mass Transfer, 28 (1985) 909-918.
[12] A. Nakayama, M.A. Hossain, An integral treatment for combined heat and mass transfer by natural
convection in a porous medium, Int. J. Heat Mass Transfer, 38 (1995) 761-765.
[13] K.A. Yih, Coupled heat and mass transfer by free convection over a truncated cone in porous media:
VWT/VWC or VHF/VMF, Acta Mech. 137 (1999) 83-97.
[14] K.A. Yih, Uniform transpiration effect on coupled heat and mass transfer in mixed convection about
inclined surfaces in porous media : the entire regime, Acta Mech. 132 (1999) 229-240.

11

[15] C.Y. Cheng, An integral approach for heat and mass transfer by natural convection from truncated cones
in porous media with variable wall temperature and concentration, Int. Commun. Heat Mass Transfer,
27 (2000) 437-548.
[16] C.Y. Cheng, Soret and Dufour effects on natural convection heat and mass transfer from a vertical cone
in a porous medium, Int. Commun. Heat Mass Transfer, 36 (2009) 1020-1024.
[17] S.U. Rahman, K. Mahgoub, A. Nafees, Natural Convective Mass Transfer from Upward Pointing Conical Surfaces in Porous Media, Chem. Eng. Commun. 194 (2007) 280-290.
[18] I. Pop, T.Y. Na, Naturnal convection of a Darcian fluid about a wavy cone, Int. Commun. Heat Mass
Transfer, 21 (1994) 891-899.
[19] I. Pop, T.Y. Na, Natural convection over a frustum of a wavy cone in a porous medium, Mech. Res.
Commun. 22 (1995) 181-190.
[20] F. C. Lai and F. A. Kulacki, Coupled heat and mass transfer by natural convection from vertical
surfaces in porous media, Int. J. Heat Mass Transfer, 34,4-5 (1991) 1189-1194.
[21] A.R. Sohouli, M. Famouri, A. Kimiaeifar, G. Domairry, Application of homotopy analysis method for
natural convection of Darcian fluid about a vertical full cone embedded in pours media prescribed surface
heat flux, Commun. Nonlinear Sci. Numer. Simul. 15 (7) (2010) 1691-1699.
[22] S. Abbasbandy, T. Hayat, Solution of the MHD Falkner-Skan flow by homotopy analysis method,
Commun. Nonlinear Sci. Numer. Simul. 14 (9-10) (2009) 3591-3598.
[23] S. Abbasbandy, Approximate solution for the nonlinear model of diffusion and reaction in porous
catalysts by means of the homotopy analysis method, Chem. Eng. J. 136 (2008) 144-150.
[24] K. Parand, M. Razzaghi, Rational Legendre approximation for solving some physical problems on
semi-infinite intervals, Phys. Scr. 69 (2004) 353-357.
[25] K. Parand, M. Dehghan, A. Pirkhedri, Sinc-collocation method for solving the Blasius equation, Phys.
Lett. A, 373 (44) (2009) 4060-4065.
[26] K. Parand, M. Dehghan, A.R. Rezaei, S.M. Ghaderi.An approximational algorithm for the solution
of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation
method, Comput. Phys. Commun. 2010; DOI:10.1016/j.cpc.2010.02.018.
[27] K. Parand, A.R. Rezaei, S.M. Ghaderi, An approximate solution of the MHD Falkner-Skan
flow by Hermite functions pseudospectral method, Commun. Nonlinear Sci. Numer. Simul. 2010;
Doi:10.1016/j.cnsns.2010.03.022.
[28] S. Abbasbandy, E. Shivanian, Exact analytical solution of a nonlinear equation arising in heat transfer,
Phys. Lett. A, 374 (4) (2010) 567-574.
[29] N. Mai-Duy, Solving high order ordinary differential equations with radial basis function networks, Int.
J. Numer. Meth. Engng. 62 (6) (2005) 824-852.
[30] E.J. Kansa, Multiquadrics-A scattered data approximation scheme with applications to computational
fluiddynamics II. Solutions to parabolic, hyperbolic and elliptic partial differential equations, Comput.
Math. Appl. 19 (8,9) (1990) 147-161.
[31] M. Zerroukat, H. Power, C.S. Chen, A numerical method for heat transfer problems using collocation
and radial basis functions, Int. J. Numer. Meth. Engng. 42 (1998) 1263-1278.
[32] N. Mai-Duy, T. Tran-Cong, Numerical solution of differential equations using multiquadric radial basis
function networks, Neural Netw. 14(2)(2001) 185-199.
[33] N. Mai-Duy, T. Tran-Cong, Numerical solution of NavierStokes equations using multiquadric radial
basis function networks, Int. J. Numer. Meth Fluids, 37 (2001) 65-86.
[34] B. Fornberg, Comparisons between pseudospectral and radial basis function derivative approximations,
IMA J. Numer. Anal. 30 (1) (2008) 149-172.
[35] N. Mai-Duy, T. Tran-Cong, Approximation of function and its derivatives using radial basis function
network methods, Appl. Math. Modelling, 27 (2003) 197-220.
[36] M.A. Golberg, C.S. Chen, H. Bowman, Some recent results and proposals for the use of radial basis
functions in the BEM, Eng. Anal. Bound. Elem. 23 (4) (1999) 285-296.
[37] M.D. Buhmann, Radial basis functions, Cambridge University Press, Cambridge, 2003.
[38] G.E. Fasshauer, Meshfree Approximation Methods with Matlab, World Scientific Publishing, Singapore,

12

(2007).
[39] M.J.D. Powell, The theory of radial basis function approximation in 1990, Advances in Numerical
Analysis, Clarendon, Oxford, (1992).
[40] M. Dehghan, A. Shokri, A meshless method for numerical solution of the one-dimensional wave equation
with an integral condition using radial basis functions, Numer. Algo. 52(3) (2009) 461-477.
[41] R.E. Carlson, T.A. Foley, The parameter r2 in multiquadric interpolation, Comput. Math. Appl. 21
(9) (1991) 29-42.
[42] S. Rippa, An algorithm for selecting a good parameter c in radial basis function interpolation, Adv.
Comput. Math. 11 (1999) 193-210.
[43] B. Fornberg, T. Dirscol, G. Wright, R. Charles, Observations on the behavior of radial basis function
approximations near boundaries, Comput. Math. Appl. 43 (2002) 473-490.
[44] S.A. Sarra, Adaptive radial basis function method for time dependent partial differential equations,
Appl. Numer. Math. 54 (2005) 79-94.
[45] Z.M. Wu, Radial basis function scattered data interpolation and the meshless method of numerical
solution of PDEs, J. Eng. Math. 19 (2) (2002), pp. 112 [In Chinese].
[46] Z.M. Wu, R. Schaback, Local error estimates for radial basis function interpolation of scattered data,
IMA J. Numer. Anal. 13 (1993) 13-27.
[47] B. Fornberg, G. Wright, Stable computation of multiquadric interpolants for all values of the shape
parameter, Computers Math. Applic. 48 (5/6) (2004) 853-867.

13

List of Figures
1
2
3
4

(a) Coordinate system for the boundary layer on a heated frustum of a cone,
(b) full cone, x0 = 0. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
DRBF approximation of f ′ (η) for different values λ = 0, 1/4, 1/3, 1/2, 3/4
and 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
IRBF approximation of f ′ (η) for different values λ = 0, 1/4, 1/3, 1/2, 3/4
and 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
kResk2 for λ = 2/3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

15
16
17
18

List of Tables
1
2
3
4
5
6

Some well–known functions that generate RBFs (r = kx − xi k = ri ), c > 0 .
A comparison between solutions obtained by [21] and the DRBF method for
f ′ (0) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Comparison between DRBF solution and Runge-Kutta solution for f ′ (η) for
λ = 1/4 and λ = 3/4 with N = 12 . . . . . . . . . . . . . . . . . . . . . . . .
A comparison between solutions obtained by [21] and the IRBF method for
f ′ (0) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Comparison between IRBF solution and Runge-Kutta solution for f ′ (η) for
λ = 1/4 and λ = 3/4 with N = 10 . . . . . . . . . . . . . . . . . . . . . . . .
kResk2 for different N and λ by IRBF . . . . . . . . . . . . . . . . . . . . .

14

20
21
21
22
22
23

xxxxxxxx
xxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxx r(x)xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx

xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxx

g

xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxx

u

xxxxx
xxxxx

X0

x

v

xxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

xxxxxxxxx
xxxxxxxxx
xxxxxxxxx
xxxxxxxxx

y
(a)

(b)

Figure 1: (a) Coordinate system for the boundary layer on a heated frustum of a cone, (b) full cone, x0 = 0.

15

Figure 2: DRBF approximation of f ′ (η) for different values λ = 0, 1/4, 1/3, 1/2, 3/4 and 1

16

Figure 3: IRBF approximation of f ′ (η) for different values λ = 0, 1/4, 1/3, 1/2, 3/4 and 1

17

Figure 4: kResk2 for λ = 2/3

18

Nomenclature
A
prescribed constant
f
similarity function for stream function temperature
g
acceleration due to gravity parameter
K
permeability of the fluid-saturated porous medium
N ux local Nusselt number
qw
surface heat flux fluid-saturated porous medium
r
local radius of the cone fluid
Rax local Raleigh number
T
temperature
T∞
ambient temperature
u, v
velocity vector along x,y axis
x, y
Cartesian coordinate system
x0
distance of start point of cone from the vertex
Greek symbols
α
thermal diffusivity the fluid-saturated porous medium
β
expansion coefficient of the fluid
η
independent dimensionless
θ
similarity function for
λ
prescribed constants
µ
viscosity of the fluid
ρ∞
density of the fluid at infinity
ψ
stream function

19

Table 1: Some well–known functions that generate RBFs (r = kx − xi k = ri ), c > 0

Name of functions
Multiquadrics (MQ)
Inverse multiquadrics (IMQ)
Thin plate (polyharmonic)Splines (TPS)
Conical splines
Gaussian (GA)
Exponential spline

Definition
√
r 2√+ c2
1/( r 2 + c2 )
(−1)k+1 r 2k log(r)
r 2k+1
exp(−cr 2 )
exp(−cr)

20

Table 2: A comparison between solutions obtained by [21] and the DRBF method for f ′ (0)

λ
0
1/4
1/3
1/2
3/4
1

Runge-Kutta
Solution[21]
0.94760
0.91130
0.90030
0.87980
0.85220
0.82760

N
10
12
10
10
12
10

DRBF method
c
DRBF
Error with RK
3.46543
0.94750
0.0001
3.943
0.91086
0.00044
4.9665
0.90038
0.0008
5.36
0.87981
0.00001
5.23
0.85227
0.00007
5.89
0.82737
0.00023

Other methods
HAM[21]
0.94783
0.91119
0.90103
0.87964
0.85242
0.82726

Table 3: Comparison between DRBF solution and Runge-Kutta solution for f ′ (η) for λ = 1/4 and λ = 3/4
with N = 12

η
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
1.5

λ = 1/4, f ′(η)
DRBF
Runge-Kutta
Solution
Solution[21]
0.910886
0.911295
0.813122
0.813604
0.720720
0.721351
0.634656
0.635531
0.555601
0.556661
0.483877
0.484997
0.419467
0.420587
0.362196
0.363276
0.311655
0.312677
0.267358
0.268264
0.228756
0.229508
0.195268
0.195878
0.166342
0.166847
0.141450
0.141837
0.120090
0.120362
0.101811
0.102025

Error
with RK
0.000409
0.000482
0.000631
0.000875
0.001061
0.001120
0.001120
0.001080
0.001022
0.000906
0.000752
0.000610
0.000505
0.000387
0.000272
0.000214

21

λ = 3/4, f ′(η)
DRBF
Runge-Kutta
Solution
Solution[21]
0.852268
0.852193
0.755678
0.755377
0.666229
0.665448
0.584176
0.582985
0.509628
0.508141
0.442519
0.440849
0.382620
0.380907
0.329539
0.327973
0.282878
0.281536
0.242076
0.241013
0.206652
0.205832
0.176015
0.175434
0.149653
0.149275
0.127026
0.126821
0.107681
0.107596
0.091167
0.091196

Error
with RK
0.000075
0.000301
0.000781
0.001191
0.001487
0.001670
0.001713
0.001566
0.001342
0.001043
0.000820
0.000581
0.000378
0.000205
0.000085
0.000029

Table 4: A comparison between solutions obtained by [21] and the IRBF method for f ′ (0)

λ
0
1/4
1/3
1/2
3/4
1

Runge-Kutta
Solution[21]
0.94760
0.91130
0.90030
0.87980
0.85220
0.82760

N
10
10
10
10
10
10

c
1.860
2.005
2.050
2.150
2.418
2.380

IRBF method
IRBF
Error with RK
0.94758
0.00002
0.91128
0.00002
0.90030
0.00000
0.87979
0.00001
0.85206
0.00014
0.82762
0.00002

Other methods
HAM[21]
0.94783
0.91119
0.90103
0.87964
0.85242
0.82726

Table 5: Comparison between IRBF solution and Runge-Kutta solution for f ′ (η) for λ = 1/4 and λ = 3/4
with N = 10

η
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
1.1
1.2
1.3
1.4
1.5
2
2.5
3
3.5
4
4.5

λ = 1/4, f ′ (η)
IRBF
Runge-Kutta
Solution
Solution[21]
0.911278
0.911295
0.813594
0.813604
0.721302
0.721351
0.635394
0.635531
0.556440
0.556661
0.484753
0.484997
0.420334
0.420587
0.362989
0.363276
0.312347
0.312677
0.267934
0.268264
0.229197
0.229508
0.195584
0.195878
0.166546
0.166847
0.141548
0.141837
0.120107
0.120362
0.101771
0.102025
0.043769
0.043951
0.018310
0.018546
0.007391
0.007610
0.002716
0.002953
0.000719
0.000962
0.000010
0.000123

Error
with RK
0.000017
0.000010
0.000049
0.000137
0.000221
0.000244
0.000253
0.000287
0.000330
0.000330
0.000311
0.000294
0.000301
0.000289
0.000255
0.000254
0.000182
0.000236
0.000219
0.000237
0.000243
0.000113

22

λ = 3/4, f ′(η)
IRBF
Runge-Kutta
Solution
Solution[21]
0.852059
0.852193
0.755254
0.755377
0.665278
0.665448
0.582736
0.582985
0.507827
0.508141
0.440510
0.440849
0.380599
0.380907
0.327647
0.327973
0.281184
0.281536
0.240696
0.241013
0.205482
0.205832
0.175053
0.175434
0.148891
0.149275
0.126408
0.126821
0.107158
0.107596
0.090757
0.091196
0.038887
0.039223
0.016225
0.016574
0.006433
0.006832
0.002255
0.002668
0.000445
0.000913
−0.00001
0.000237

Error
with RK
0.000134
0.000123
0.000170
0.000249
0.000314
0.000339
0.000308
0.000326
0.000352
0.000317
0.000350
0.000381
0.000384
0.000413
0.000438
0.000439
0.000336
0.000349
0.000399
0.000413
0.000468
0.000247

λ
0
1/4
1/3
1/2
3/4
1

N=5
1.161037
0.851976
0.770207
0.637068
0.469112
0.355474

Table 6: kResk2 for different N and λ by IRBF

N=6
0.306257
0.278068
0.269714
0.216095
0.230425
0.221336

N=8
0.016229
0.018364
0.018842
0.016813
0.020280
0.020463

23

N=10
0.000108
0.000363
0.000494
0.000587
0.001051
0.001077

N=12
0.77e − 5
0.52e − 4
0.24e − 4
0.83e − 5
0.15e − 4
0.31e − 4

N=15
0.11e − 7
0.27e − 6
0.13e − 5
0.47e − 6
0.83e − 6
0.39e − 6

