Next Article in Journal
Modified SHARK Cipher and Duffing Map-Based Cryptosystem
Previous Article in Journal
Quadrotor Real-Time Simulation: A Temporary Computational Complexity-Based Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Median Bernoulli Numbers and Ramanujan’s Harmonic Number Expansion

Department of Mathematics, University of Taipei, No. 1, Ai-Guo West Road, Taipei 100234, Taiwan
Mathematics 2022, 10(12), 2033; https://doi.org/10.3390/math10122033
Submission received: 20 March 2022 / Revised: 10 May 2022 / Accepted: 10 June 2022 / Published: 12 June 2022

Abstract

:
Ramanujan-type harmonic number expansion was given by many authors. Some of the most well-known are: H n γ + log n k = 1 B k k · n k , where B k is the Bernoulli numbers. In this paper, we rewrite Ramanujan’s harmonic number expansion into a similar form of Euler’s asymptotic expansion as n approaches infinity: H n γ + c 0 ( h ) log ( q + h ) k = 1 c k ( h ) k · ( q + h ) k , where q = n ( n + 1 ) is the nth pronic number, twice the nth triangular number, γ is the Euler–Mascheroni constant, and c k ( x ) = j = 0 k k j c j x k j , with c k is the negative of the median Bernoulli numbers. Then, 2 c n = k = 0 n n k B n + k , where B n is the Bernoulli number. By using the result obtained, we present two general Ramanujan’s asymptotic expansions for the nth harmonic number. For example, H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 j = 0 b j ( r ) ( q + 1 3 ) j 1 / r as n approaches infinity, where b j ( r ) can be determined.
MSC:
41A60; 11B83; 05A19

1. Introduction

Leonhard Euler in 1755 applied the Euler–Maclaurin sum formula to find the famous standard Euler asymptotic expansion for H n as n :
H n γ + log n k = 1 B k k · n k ,
where B k is the Bernoulli number defined by t e t 1 = k = 0 B k k ! t k , and γ = 0.57721 is the Euler–Mascheroni constant.
Ramanujan [1] proposed the following asymptotic expansion for H n :
H n γ + 1 2 log ( 2 m ) + 1 12 m 1 120 m 2 + 1 630 m 3 1 1680 m 4 + 1 2310 m 5 191 360360 m 6 + 29 30030 m 7 2833 1166880 m 8 + 140051 17459442 m 9 ,
where m = n ( n + 1 ) / 2 is the n-th triangular number. However, Ramanujan did not give any formulas for the general terms and also without any proof. Rewrite the above formula as the following notation:
H n γ + 1 2 log ( 2 m ) + k = 1 R k m k .
In 2008, Villarino [2] established an explicit expression for the coefficient sequence ( R k ) :
R k = ( 1 ) k 1 2 k · 8 k j = 0 k k j ( 4 ) j B 2 j ( 1 / 2 ) ,
where B k ( x ) are the Bernoulli polynomials defined by t e x t e t 1 = n = 0 B n ( x ) t n n ! . In 2015, Chen and Cheng [3] reconsidered Ramanujan’s formula and gave the following recurrence relation for ( R k ) :
R 1 = 1 12 , R k = 1 2 k 1 4 k B 2 k 2 k j = 1 k 1 2 j R j 2 k j 1 j 1 , k 2 .
In 2019, Chen [4] improved the recurrence relation as
R k = 1 2 k + 1 k 1 2 k + 1 j = 1 k 1 2 j + 1 R j 2 k j 2 k 2 j + 1 , for   k 2 .
Another Ramanujan-type harmonic number expansion was given by Wang [5] in 2018,
H n γ + 1 2 log ( 2 m + h ) + k = 1 α k ( h ) ( 2 m + h ) k ,
where h is a parameter and ( α k ( h ) ) is a coefficient sequence
α k ( h ) = h k 2 k + j = 1 k k 1 j 1 R j 2 j h k j .
In this paper, we rewrite Ramanujan’s harmonic number expansion into a similar form of Euler’s asymptotic expansion:
H n γ + c 0 log q k = 1 c k k · q k ,
where q = n ( n + 1 ) = 2 m is the nth pronic number, twice the nth triangular number. In fact, we prove that the number c k is the negative of the median Bernoulli number. The median Bernoulli number is studied by the author [6] in 2005. Then, we have for n 0 ,
c n = 1 2 k = 0 n n k B n + k .
Moreover, let
c k ( x ) = j = 0 k k j c j x k j .
Then, we could rewrite Wang’s expansion Equation (7) as follows:
H n γ + c 0 ( h ) log ( q + h ) k = 1 c k ( h ) k · ( q + h ) k .
We give simpler asymptotic expansion representations for H n using Equations (10) and (11), which effectively integrate the results of Villarino, Chen and Cheng, Chen, and Wang (see Equations (4)–(7)) and make their representations more meaningful. We discuss some properties of the numbers c n and the polynomials c n ( x ) in Section 3 and Section 4, respectively. For example, the Hankel determinant of c n ( x ) for any x can be evaluated as
2 n + 1 det 0 i , j n ( c i + j ( x ) ) = det 0 i , j n ( B 2 i + 2 j ( 1 / 2 ) ) .
Furthermore, Chen [7] gave a new asymptotic expansion. For any nonzero real number r, the n-th harmonic number H n may have an asymptotic expansion as n approaches infinity:
1 2 log ( 2 m ) + γ + 1 12 m j = 0 a j ( r ) m j 1 / r ,
where the parameters a j ( r ) satisfy the following recurrence relation
a 0 ( r ) = 1 , a j ( r ) = 1 j k = 1 j k ( 1 + r ) j ( 12 R k + 1 ) a j k ( r ) , j N .
Inspired by this, we give a more general asymptotic expansion in Section 5 using Equation (12). Given r , h real numbers with r 0 , h 1 / 3 , we get
H n γ + c 0 ( h ) log ( q + h ) 3 h 1 6 ( q + h ) 1 + j = 1 a j ( r , h ) ( q + h ) j 1 / r , n ,
We know that the formula with h = 0 is Equation (14) (see ([7], Theorem 2.3)). Since c 1 ( h ) = 3 h 1 6 , h = 1 / 3 will remove the ( q + h ) 1 term. This will improve the approximation. Thus, it can be seen that there are a lot of investigations for the h = 1 / 3 case, see [4,8,9,10].
If h = 1 / 3 , then the asymptotic expansion will become
H n γ + 1 2 log q + 1 3 1 180 q + 1 3 2 1 + j = 1 b j ( r ) q + 1 3 j 1 / r .
The parameters a j ( r , h ) and b j ( r ) in Equations (15) and (16) are determined by some recurrence relations, which will be illustrated in Theorems 2 and 3, respectively. At the end of this paper, we will compare how close these asymptotic formulas are to H n .

2. Median Bernoulli Numbers and R k

Set a 0 , n = B n , for n 0 . And for n 1 , k 0 ,
a n , k = a n 1 , k + a n 1 , k + 1 ,
or equivalently,
a n , k = j = 0 n n j a 0 , k + j .
The corresponding matrix is represented as follows.
1 1 / 2 1 / 6 0 1 / 30 0 1 / 42 0 1 / 2 1 / 3 1 / 6 1 / 30 1 / 30 1 / 42 1 / 42 1 / 6 1 / 6 2 / 15 1 / 15 1 / 105 1 / 21 0 1 / 30 1 / 15 8 / 105 4 / 105 1 / 30 1 / 30 1 / 105 4 / 105 0 1 / 42 1 / 21 1 / 42 1 / 42 0
This matrix is called the “ B S -matrix” in [6], which is a special Euler–Seidel matrix. Let
c n = a n + 1 , n
be the lower diagonal sequence of the B S -matrix. The number c n is the negative of the median Bernoulli number K n , which is the upper diagonal sequence of the B S -matrix, i.e., K n = a n , n + 1 = c n (ref. [6]). Therefore, by [6], and Equations (8), (15) and (16), we have
c n = 1 2 k = 0 n n k B n + k = k = 0 n + 1 n + 1 k B n + k = k = 0 n n k B n + 1 + k .
Let the ordinary generating function of c n as follows.
m ( x ) = n = 0 c n x n + 1 .
Let ψ ( x ) be the formal Laplace transform of t / sinh t . Then, the following relation was obtained ([6], Theorem 4.2, Equation (29))
2 x · ψ ( x ) = m 4 x 2 1 x 2 .
Since t / sinh t = n = 0 4 n B 2 n ( 1 / 2 ) t 2 n / ( 2 n ) ! , we have that for n 0 ([6], Equation (32)),
2 2 n + 1 ( 1 ) n c n = j = 0 n n j ( 1 ) j 2 2 j B 2 j ( 1 / 2 ) .
Using Villarino’s explicit formula for R k , Equation (4), we have for k 1 ,
2 k · k · R k = j = 0 k k j ( 1 ) k + j 2 2 j 2 k 1 B 2 j ( 1 / 2 ) = c k .
This implies that
H n γ + c 0 log q k = 1 c k k q k , n .
On the other hand, we substitute R k as Equation (20) in Wang’s formula for α k ( h ) , (see Equation (8)), we have for k 1 ,
k · α k ( h ) = h k 2 j = 1 k k 1 j 1 R j 2 k h k j k = h k 2 + j = 1 k k 1 j 1 h k j c j j · k = j = 0 k k j h k j c j = c k ( h ) .
Therefore, we conclude our result in the following.
Theorem 1.
For n , we have
H n γ + c 0 ( h ) log ( q + h ) k = 1 c k ( h ) k · ( q + h ) k ,
where q = n ( n + 1 ) is the nth pronic number,
c k ( x ) = j = 0 k k j c j x k j , a n d c n = 1 2 k = 0 n n k B n + k .

3. Some Properties of c n

Let the ordinary generating function of B n as follows.
b ( x ) = n = 0 B n x n + 1 .
Using the relation between the ordinary generating functions of a 0 , n , a n , n , and a n , n + 1 of the B S -matrix, we have the following relation ([6], Theorem 4.2, Equation (29))
b ( x ) = 1 + 2 x m x 2 1 + x .
Then, the following identity is obtained ([6], Equation (27)).
j = 0 n / 2 n j j n n j c j = B 1 , if n = 1 , B n , if n 2 .
In the above formula, the formula obtained by substituting n = 2 k appears in the recurrence relation of R k given by Chen and Cheng [3] in 2015 (see Equation (5)). Furthermore, if we substitute n = 2 k + 1 into the above identity, we obtain Equation (6) given by Chen [4] in 2019.
There are a lot of properties of c n obtained from [6]. For example, let the denominators and the numerators of the rational number c n be D n , N n , respectively. We have the following properties ([6], Theorem 1.1):
  • The denominator D n is a square-free integer.
  • The set of the all odd prime divisors of D n is
    p : odd prime n m p 1 2 n 2 m 1 , m N .
  • The denominator D n is an odd integer, for n 2 .
  • The largest power of 2 that divides the numerator N n is 2 n 1 2 .
The ordinary generating function m ( x ) has the following continued fraction representation ([6], Theorem 5.5)
m ( x ) = c 0 x 1 + a 0 x b 1 x 2 1 + a 1 x b 2 x 2 1 + a 2 x + ,
where for n 0 ,
a n = 8 n 4 + 8 n 3 + 6 n 2 + 2 n 1 ( 4 n + 3 ) ( 4 n 1 ) , b n + 1 = ( 2 n + 1 ) 4 ( n + 1 ) 4 ( 4 n + 1 ) ( 4 n + 3 ) 2 ( 4 n + 5 ) .
Using this representation, we have the Hankel determinant of c n (ref. [6], Theorem 5.5)
det 0 i , j n ( c i + j ) = 1 2 n + 1 j = 1 n ( 2 j 1 ) 4 j 4 ( 4 j 3 ) ( 4 j 1 ) 2 ( 4 j + 1 ) n j + 1 .
Since the finite product in Equation (28) is the Hankel determinant of B 2 n ( 1 / 2 ) (see [6], Equation (41)), we have
det 0 i , j n B 2 i + 2 j ( 1 / 2 ) = 2 n + 1 det 0 i , j n ( c i + j ) .
By Equation (19) and an integral representation of B 2 n ( 1 / 2 ) ([11], Equation (28))
B 2 n ( 1 / 2 ) = ( 1 ) n π 0 t 2 n sech 2 ( π t ) d t ,
we have an integral representation of c n , for n 0 ,
c n = ( 1 ) n π 2 2 n + 1 0 ( 4 t 2 + 1 ) n sech 2 ( π t ) d t .

4. Some Properties of c n ( x )

We first list c n ( x ) for n = 0 , 1 , 2 , , 5 .
c 0 ( x ) = 1 2 , c 1 ( x ) = 1 6 + x 2 , c 2 ( x ) = 1 15 x 3 + x 2 2 , c 3 ( x ) = 4 105 + x 5 x 2 2 + x 3 2 , c 4 ( x ) = 4 105 16 x 105 + 2 x 2 5 2 x 3 3 + x 4 2 , c 5 ( x ) = 16 231 + 4 x 21 8 x 2 21 + 2 x 3 3 5 x 4 6 + x 5 2 .
Differentiating Equation (11) with respect to x we obtain
d d x c n ( x ) = k = 0 n n k c n k k x k 1 = n k = 0 n 1 n 1 k c n 1 k x k = n c n 1 ( x ) .
Therefore,
x y c n ( t ) d t = c n + 1 ( y ) c n + 1 ( x ) n + 1 .
On the other hand, we use Equation (30) to get an integral representation of c n ( x ) :
c n ( x ) = ( 1 ) n π 2 2 n + 1 0 ( 4 t 2 4 x + 1 ) n sech 2 ( π t ) d t .
Let us consider the function c n ( x + y ) . We express ( x + y ) k as its binomial expansion.
c n ( x + y ) = k = 0 n n k c n k ( x + y ) k = k = 0 n c n k = 0 k k x y k .
We interchange the order of summation and the inner sum becomes c n ( y ) :
c n ( x + y ) = = 0 n k = n n k k c n k x y k = = 0 n n x c n ( y ) .
Thus, we have
c n ( x + y ) = k = 0 n n k c k ( x ) y n k .
Using the inversion binomial theorem to Equation (19) we have
4 n B 2 n ( 1 / 2 ) = k = 0 n n k 2 2 k + 1 c k = 2 2 n + 1 c n ( 1 / 4 ) .
This implies that
2 c n ( 1 / 4 ) = B 2 n ( 1 / 2 ) .
Applying the above identity and Equation (29), we know that the Hankel determinant of c n ( 1 / 4 ) is the same as the Hankel determinant of c n . However, we use ([12], Proposition 1), indeed that for any value of x,
det 0 i , j n ( c i + j ) = det 0 i , j n ( c i + j ( x ) ) .

5. New Asymptotic Expansions

To derive our new asymptotic expansions are inspired by ([7], Theorem 2.3). We need the following lemma.
Lemma 1
([7], Lemma 1). If j = 0 q j x j is an asymptotic expansion for g ( x ) as x approaches infinity. Given any real number r, the parameters Q j ( r ) are defined by Q 0 ( r ) = 1 and for j N ,
Q j ( r ) = 1 j k = 1 j k ( 1 + r ) j q k Q j k ( r ) .
Then j = 0 Q j ( r ) x j is an asymptotic expansion for g ( x ) r .
Our new asymptotic expansions are derived from Equation (23). It is note that c 1 ( h ) = 3 h 1 6 . Therefore, we divide into two cases depending on whether h is 1 / 3 or not.
Theorem 2.
Let r and h be any given real numbers with r 0 and h 1 / 3 . The n-th harmonic number H n has the following asymptotic expansion as n approaches infinity:
γ + c 0 ( h ) log ( q + h ) 3 h 1 6 ( q + h ) 1 + j = 1 a j ( r , h ) ( q + h ) j 1 / r
where the parameters a j ( r , h ) given by the recurrence relation
a 0 ( r , h ) = 1 , a j ( r , h ) = 1 j k = 1 j k ( 1 + r ) j 6 c k + 1 ( h ) ( 3 h 1 ) ( k + 1 ) a j k ( r , h ) , j 1 .
Proof. 
Rewrite Equation (37) as the following representation:
6 ( q + h ) 3 h 1 H n γ c 0 ( h ) log ( q + h ) 1 + j = 1 a j ( r , h ) ( q + h ) j 1 / r .
In view of Equation (23), we have
6 ( q + h ) 3 h 1 H n γ c 0 ( h ) log ( q + h ) 1 + 6 3 h 1 k = 1 c k + 1 ( h ) ( k + 1 ) ( q + h ) k .
Comparing the above two expressions, we know that
1 + 6 3 h 1 k = 1 c k + 1 ( h ) ( k + 1 ) ( q + h ) k r 1 + j = 1 a j ( r , h ) ( q + h ) j .
We apply Lemma 1 and get the result we want. □
Using a similar approach, we can easily derive the following theorem for the situation h = 1 / 3 .
Theorem 3.
Given a real number r with r 0 . The n-th harmonic number H n has the asymptotic expansion as n approaches infinity:
γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 1 + j = 1 b j ( r ) ( q + 1 3 ) j 1 / r
where the parameters b j ( r ) are defined by the following relation
b 0 ( r ) = 1 , b j ( r ) = 1 j k = 1 j k ( 1 + r ) j 180 c k + 2 ( 1 3 ) k + 2 b j k ( r ) , j 1 .
Chen [7] discussed many properties of the h = 0 case. Therefore, we mainly deal with the case of h = 1 / 3 here.
The first few parameters b j ( r ) are:
b 0 ( r ) = 1 , b 1 ( r ) = 32 63 r , b 2 ( r ) = 3701 7938 r + 512 3969 r 2 , b 3 ( r ) = 7264240 8251551 r 59216 250047 r 2 16384 750141 r 3 , b 4 ( r ) = 47882328785 18021387384 r + 2311659673 4158781704 r 2 + 947456 15752961 r 3 + 131072 47258883 r 4 , b 5 ( r ) = 8014919889976 709592128245 r 749340134980 425755276947 r 2 789621116 4678629417 r 3 30318592 2977309629 r 4 4194304 14886548145 r 5 .
For r = 1 in Equation (39), the resulting asymptotic expansion is as follows ([10], Equation (3.24)):
H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 + 8 2835 ( q + 1 3 ) 3 5 1512 ( q + 1 3 ) 4 + 592 93555 ( q + 1 3 ) 5 796801 43783740 ( q + 1 3 ) 6 + 268264 3648645 ( q + 1 3 ) 7
as n .
For r = 1 in Equation (39), we obtain a new asymptotic expansion:
H n γ + 1 2 log ( q + 1 3 ) 180 ( q + 1 3 ) 2 + 640 7 ( q + 1 3 ) 26770 441 + 36602240 305613 ( q + 1 3 ) 97247611025 250297047 ( q + 1 3 ) 2 + 27515011460000 15768713961 ( q + 1 3 ) 3 1
as n .
For r = 23 / 2 in Equation (39), we obtain a new asymptotic expansion:
H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 1 + 368 63 ( q + 1 3 ) + 185725 15876 ( q + 1 3 ) 2 + 3674204 305613 ( q + 1 3 ) 3 5793677 728136864 ( q + 1 3 ) 4 + 1021070020123 31537427922 ( q + 1 3 ) 5 + 2 / 23
as n .
From a computational point of view, the formulas Equations (42) and (43) are better than Equation (41).
It follows from Equations (41)–(43) that for n ,
H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 + 8 2835 ( q + 1 3 ) 3 5 1512 ( q + 1 3 ) 4 + 592 93555 ( q + 1 3 ) 5 : = u n ,
H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 + 640 7 ( q + 1 3 ) 26770 441 + 36602240 305613 ( q + 1 3 ) : = v n ,
H n γ + 1 2 log ( q + 1 3 ) 1 180 ( q + 1 3 ) 2 1 + 368 63 ( q + 1 3 ) + 185725 15876 ( q + 1 3 ) 2 + 3674204 305613 ( q + 1 3 ) 3 2 / 23 : = w n .
From Table 1, we observe that, among approximation formulas Equations (44)–(46), for n 1 , the formula Equation (46) would be the best one. There seems to be an optimal real number r in Equation (39), and when we substitute it in this formula, the resulting approximation should be optimal. We guess that this real number r should be close to 11.502534

Funding

The author was funded by the Ministry of Science and Technology, Taiwan, Republic of China, through grant MOST 110-2115-M-845-001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berndt, B. Ramanujan’s Notebooks; Springer: New York, NY, USA, 1998; Volume 5. [Google Scholar]
  2. Villarino, M.B. Ramanujan’s harmonic number expansion into negative powers of a triangular number. J. Inequal. Pure. Appl. Math. 2008, 9, 89. [Google Scholar]
  3. Chen, C.-P.; Cheng, J.X. Ramanujan’s asymptotic expansion for the harmonic numbers. Ramanujan J. 2015, 38, 123–128. [Google Scholar] [CrossRef]
  4. Chen, C.-P. Asymptotic series related to Ramanujan’s expansion for the harmonic number. Proc. Indian Aca. Sci. (Math. Sci.) 2019, 129, 77. [Google Scholar] [CrossRef]
  5. Wang, W. Harmonic number expansions of the Ramanujan type. Results Math. 2018, 73, 161. [Google Scholar] [CrossRef]
  6. Chen, K.-W. A summation on Bernoulli numbers. J. Number Theory 2005, 111, 372–391. [Google Scholar] [CrossRef]
  7. Chen, C.-P. Ramanujan’s formula for the harmonic number. Appl. Math. Comput. 2018, 317, 121–128. [Google Scholar] [CrossRef]
  8. Mortici, C.; Chen, C.-P. On the harmonic number expansion by Ramanujan. J. Inequal. Appl. 2013, 222, 1–10. [Google Scholar] [CrossRef]
  9. Mortici, C.; Villarino, M.B. On the Ramanujan-Lodge harmonic number expansion. Appl. Math. Comput. 2015, 251, 423–430. [Google Scholar] [CrossRef]
  10. Chen, C.-P. On the coefficients of asymptotic expansion for the harmonic number by Ramanujan. Ramanujan J. 2016, 40, 279–290. [Google Scholar] [CrossRef]
  11. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1953; Volume 1. [Google Scholar]
  12. Junod, A. Hankel determinants and orthogonal polynomials. Expo. Math. 2003, 21, 63–74. [Google Scholar] [CrossRef]
Table 1. Comparison of approximation Formulas (44)–(46).
Table 1. Comparison of approximation Formulas (44)–(46).
n u n H n v n H n w n H n
1 4.625 × 10 5 1.997 × 10 5 1.405 × 10 6
10 9.735 × 10 15 6.332 × 10 15 6.750 × 10 17
10 2 1.713 × 10 26 1.129 × 10 26 2.162 × 10 30
10 3 1.809 × 10 38 1.192 × 10 38 3.805 × 10 42
10 4 1.819 × 10 50 1.198 × 10 50 3.841 × 10 54
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, K.-W. Median Bernoulli Numbers and Ramanujan’s Harmonic Number Expansion. Mathematics 2022, 10, 2033. https://doi.org/10.3390/math10122033

AMA Style

Chen K-W. Median Bernoulli Numbers and Ramanujan’s Harmonic Number Expansion. Mathematics. 2022; 10(12):2033. https://doi.org/10.3390/math10122033

Chicago/Turabian Style

Chen, Kwang-Wu. 2022. "Median Bernoulli Numbers and Ramanujan’s Harmonic Number Expansion" Mathematics 10, no. 12: 2033. https://doi.org/10.3390/math10122033

APA Style

Chen, K. -W. (2022). Median Bernoulli Numbers and Ramanujan’s Harmonic Number Expansion. Mathematics, 10(12), 2033. https://doi.org/10.3390/math10122033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop