arXiv:1811.03957v3 [math.GM] 11 Jun 2020 ∗∗ h amncoclao n nsaitclmcais[7–9] mechanics osci statistical harmonic in Gaussia the and of oscillator the state harmonic physics the ground possible In the is of it 4]. density variables [3, probability of statistics change and slight probability a in integr do Gaussian we The when gam 2]. Indeed, so-called [1, the Weierstrass involve and integral Gauss of Legendre, type this of Solutions n ifrnito ne h nerlsg 1] e snow us Let method [13]. the sign used integral he the under them methods differentiation among Others and integral, [10]. Gaussian method the integral solving double the being books where and o ipiiy eaecnieigta h ubrzr sin is zero number the that considering are we simplicity, For § ‡ † ∗ ¶ [email protected] [email protected] [email protected] [email protected] [email protected] [email protected] is ercl h ouino olwn asinintegral Gaussian following of solution the recall we First asinitga lokona rbblt nerli the is integral probability as known also integral Gaussian Γ 4 eatmnod ieca xtseNtri,Universidade Naturais, e Ciˆencias Exatas de Departamento ( x α ) 2 ∈ eatmnod ieca enlga,Uiesdd Fed Universidade Tecnologias, Ciˆencias e de Departamento stegmafnto.Tegmafnto a ra relevanc great has function gamma The function. gamma the is e ehdt ov h asinItga:ASlto Inspi Solution A Integral: Gaussian the Solve to Method New A R , rsnsanwfr fsligteGusa nerlhvn t having integral Gaussian the solving conclusion func of In form factorial new calculus. and a fractional gamma presents in the application of eva simple inve to class a have us certain we enabled a method, pr which for the framework, more solutions for Tsallis shown test in has a density As work number this literature. the in in method found proposed introduced also The have We paramet induction. newsworthy. this very that result demonstrated a have function, We function. simple a of htdpnso a on depends that 1 eatmnod ´sc,Uiesdd oEtd oRoGran Rio do Estado do F´ısica, Universidade de Departamento α Keywords: nti ae,w aepooe e ehdfrsligteGa the solving for method new a proposed have we paper, this In 3 6= eatmnod ´sc,Uiesdd eea aPara´ıba da Federal F´ısica, Universidade de Departamento .Teeaesvrlmtost ov hstp fitga,t integral, of type this solve to methods several are There 0. asinitga,Tyo eis pca ucin,frac functions, special series, Taylor integral, Gaussian ´ zr .Sales L. L´azaro .C .Souza M. C. T. n ne,w aefudagnrlslto o hstp finteg of type this for solution general a found have we index, Z − + ∞ ∞ Z − x + 2 ∞ 1 n ∞ , 4 ∗ , e x § oaa .Silva A. Jonatas , − noi .S Farias S. D. Antonio , 2 n α e x 2 − Z dx α .INTRODUCTION I. Dtd ue1,2020) 15, June (Dated: − x + 2 ∞ ∞ dx = ofidispriinfunction. partition its find to , e − nerlapasfeunl nqatmmcais[] ofi to [5], mechanics quantum in frequently appears integral n osdramr eea asinitga ie[,1,15] 14, [9, like integral Gaussian general more a consider afntositoue yElri 8hcnuyadimprove and century 18th in Euler by introduced functions ma lddi h e fntrlnmes Here, numbers. natural of set the in cluded ltr ntept nerlfruain[] ofidteprop the find to [6], formulation integral path the in llator, ( α sbellow: as , eea ua oSemi- do Rural Federal = 2 lhsawd ag fapiain nsvrlaeso knowle of areas several in applications of range wide a has al x n 2 nerlo h ucinexp function the of integral ocmuetenraiigcntn ftenra distribut normal the of constant normalizing the compute to dx 2 − a efudi 1,1] ordsnhszdeee asfor ways eleven synthesized Conrad 12]. [11, in found be can rlRrld Semi- do Rural eral Γ 2 fFuirtasom,Siln’ oml,cnorintegr contour formula, Stirling’s transforms, Fourier of n , α ebleei h eeac fti okbcueit because work this of relevance the in believe we , † 1 = lˆneaP Bento P. Eliˆangela , ) 2 2 !! n n edfeeta aclsa tool. a as calculus differential he 2 2 r rrpeet h alrsre ofcet fthis of coefficients series Taylor the represents er oeTerm htaepoe ymathematical by proved are that Theorems some + in r eie.Mroe,w aepresented have we Moreover, derived. are tions + 4 tgtdannetnievrinfrteparticle the for version non-extensive a stigated α 1 ut h ucinlt ftemto.Besides, method. the of functionality the luate 1 , √ ¶
α 85-7,J˜oPso-B rzland Brazil Jo˜ao Pessoa-PB, 58059-970, , ed ot,56020 Mossor´o-RN, Brazil 59610-210, Norte, do de π 2 n tai .Lavor P. Ot´avio and , n , π 2 + , cia n cesbeta oemethods some than accessible and actical o h eeomn fnwfntosta a be can that functions new of development the for e 1 ( inlderivative. tional ∀ sinitga.Itouigaparameter a Introducing integral. ussian , n ( ∈ ∀ ems elkonadwdl sdi text- in used widely and known well most he rd,59000 a o ersR,Brazil Ferros-RN, dos Pau 59900-000, Arido, ´ rd,57000 Cara´ubas-RN, Brazil 59780-000, Arido, n N ´ ∈ ) . N 3 a nprdb alrseries Taylor by inspired ral ) , ‡ , Hidalyn , ( e yTyo Series Taylor by red − x 4 2 , ) ∗∗ vrteetr ie( line entire the over x !i h obefactorial double the is !! gtrof agator − dthe nd ∞ ation, by d dge. , ∞ ion (2) (1) (3) ). 2 applied directly in physics. Normally this function is present in problems of physics such as, for example, in the normalization of Coulomb wave functions and the calculation of probabilities in statistical mechanics [16]. Notice that Eqs. (2) and (3) admits a solution for x2n+1, however, we present only the solutions for x2n because the method which will be presented in this work requires it. Inspired by the ways to solve the Gaussian integral as presented by Conrad, this work aims to present a new method for solving the Gaussian integral based on the numerical evolution of the expansion coefficients in Taylor Series of a simple function. This is relevant because enable us to treat certain mathematical and physical problems from another perspective having the differential calculus as a tool. The method consists of introducing a parameter that has the role of describing the evolution of the numerical sequence of Gaussian integral. Next, we show that this parameter can be thought as the coefficients of the expansion of function 1/2 f (x) = (1 x)− . Having done this, we were able to introduce some Theorems that are proved by mathematical induction. As a test for the− method, we applied the results due to Theoremsto obtain a non-extensiveversion for the particle number density via Tsallis statistics. Also, solutions for the gamma function of the form Γ(1/2 n), and the factorial function of a kind (n + 1/2)!, where n N, are derived. Besides, we present an application possibility in fractional± calculus using the definitions of fractional derivative∈ according to Riemann-Liouville and Caputo. The paper is structured as follows. In Section II, we show the entire construction process of the new method for solving the Gaussian integral. Simple applications in Tsallis statistical framework, special functions (gamma and factorial), and fractional calculus are presented in Section III. Finally, the conclusions are shown in Section IV.

II. THE METHOD

Consider the Gaussian integral as follows:

+∞ 2n αx2 I2n = x e− dx , ( n N;α > 0) . (4) ∞ ∀ ∈ Z− 1/2 3/2 5/2 Solving I2n by conventional methods, we arrive at the following results: I0 = √π/α , I2 = √π/2α , I4 = 3√π/4α , 7/2 I6 = 15√π/8α and so on. We can generalize these results by introducing a parameter that depends on n. The new method proposed here is based on the following definition:

√π I2n γ2n 2n+1 , (5) ≡ α 2 where γ2n is a parameter to be determined. Note that γ0 = 1, γ2 = 1/2, γ4 = 3/4, γ6 = 15/8 and so forth. The parameter γ2n has 1/2 an apparently unpredictable sequence, but it becomes evident when we perform the expansion of the function f (x) = (1 x)− in a Taylor series around x = 0, as bellow: −

∞ n 1/2 n d (1 x)− x f (x) = ∑ − n . (6) n 0 dx n! = x=0

Analyzing the first coefficients of the above expansion, we arrive at the following results: f (0)= 1, f ′(0)= 1/2, f ′′(0)= 3/4, f ′′′(0)= 15/8, f ′′′′(0)= 105/16 and so on. We now intend to find an expression for the parameter γ2n from the results shown above. Note that the coefficients of the expansion generate identical numbers those generated by the evolution of the parameter γ2n. Thus, it is reasonable to establish the following definition:

Definition 1. For n N, γ2n can be defined as ∈ n 1/2 d (1 x)− γ2n − n . (7) ≡ dx x=0

Lemma 1. By Definition 1 for all n N, we have ∈ n 1 1 γ = − + k . (8) 2n ∏ 2 k=0   Proof. The proof of this result is by induction on n. More generally, the product notation is defined as

n ∏ xi = xm xm+1 xn 1 xn . (9) i=m · ··· − · 3

If m > n, the product is an empty product, whose value is 1. The base case (n = 0) is an empty product, hence γ0 = 1. From Definition 1, for now with x = 0, it is easy to see that the inductive hypothesis is 6 n 1/2 n 1 d (1 x)− − 1 1 n − = + k (1 x)− 2 − . (10) dxn ∏ 2 − k=0   Then, we shall show that statement (10) it is true for n + 1. Therefore

n+1 1/2 n 1 d (1 x)− d − 1 1 n − = + k (1 x)− 2 − , dxn+1 dx ∏ 2 − "k=0   # n 1 1 (n+1) = + k (1 x)− 2 − , (11) ∏ 2 − k=0   where we used the inductive hypothesis to obtain the above result. Thus, at x = 0, Eq. (11) may take the form

n 1 − 1 γ2n = + k , n N . (12) ∏ 2 ∀ ∈ k=0  

Definition 2. From Eq. (2) we can write

(2n 1)!! γ¯2n − , n N . (13) ≡ 2n ∀ ∈ The double factorial m!! is defined as follow [16]:

m (m 2)...5 3 1 m > 0 odd m!! = m · (m − 2)...6 · 4 · 2 m > 0 even (14)  1 · − · · m = 1, 0 .  −

Theorem 1. Let I2n given by Eq. (2), then 

√π N I2n = γ2n 2n+1 , n , (15) α 2 ∀ ∈ where γ2n is given by Definition 1.

To prove the Theorem 1, we will take into account the Definition 2. In other words, we just need to show that γ2n = γ¯2n, for all n N. ∈ Proof. The base case (n = 0) shows that γ0 = γ¯0, since from Definition 2, ( 1)!! = 1. Suppose that γ2n = γ¯2n it is true for n = q with q N, then the inductive hypothesis is given by − ∈

γ2q = γ¯2q . (16)

We now must show that γ2n = γ¯2n it is also true for n = q + 1. So we have

γ2(q+1) = γ¯2(q+1) . (17)

By Lemma 1, with n = q + 1 and k = i, the left-hand side of relation (17) becomes

q 1 γ = + i , 2(q+1) ∏ 2 i=0   1 1 1 1 = + 1 + (q 1) + q , 2 2 ··· 2 − 2      1 = γ + q . (18) 2q 2   4

Applying the inductive hypothesis, we get

1 γ = γ¯ + q . (19) 2(q+1) 2q 2   Using the Definition 2 in Eq. (19), we obtain

(2q 1)!! 1 γ = − + q , 2(q+1) 2q 2   1 = [1 3 5 (2q 1)(2q + 1)] , 2q+1 · · ··· − = γ¯2(q+1) . (20) This completes the proof.

An immediate consequence of the Theorem 1 is that we can establish a relationship between the double factorial and the parameter γ2n, as bellow:

n 1/2 n n d (1 x)− N (2n 1)!! = 2 γ2n = 2 − n , n . (21) − dx ∀ ∈ x=0

For example, when n = 1, we have 1!! = 1, since γ2 = 1/2. For n = 2, we have 3!! = 3, since γ4 = 3/4 and so on. The method developed here is relevant because it presents a new form of solving the Gaussian integral of the type I2n. Since there are other methods, this has shown more practical than some methods found in the literature. For example, in Eq. (3), we need to solve a gamma function Γ(t) as shows the Definition 3, whereas our method just needs to derive a simple function. 1/2 Besides, we demonstrated that parameter γ2n represents the Taylor series coefficients of the function f (x) = (1 x) . This − − means that the evolution of Gaussian integral of a kind I2n has a strong relationship to the expansion terms of f (x), a result very newsworthy. One aspect that makes this method accessible is the fact that, for example, in physical problems n is not usually very large, which facilitates its application. Hereafter, we briefly present an application of our method in physics, specifically in the determination of the particle number density in Tsallis framework. Moreover, we present solutions for the gamma and factorial functions (special functions) in terms of parameter γ2n, and we also show a simple application in fractional calculus.

III. SIMPLE APPLICATIONS

A. Particle number density in Tsallis framework

We aim to present a simple application of the method proposed in this work. We chose to determinethe particle numberdensity in the Tsallis framework because an ideal scenario emerges in which it is possible to notice the functionality of the method. In 1988, Constantino Tsallis proposed a possible generalization of the Boltzmann-Gibbs (BG) entropy [17]. The proposed new entropy is expressed by

Ω kB q Sq = 1 p , (22) q 1 − ∑ i − i=1 ! where kB is the Boltzmann constant, pi is the probability of the system to be found in the microstate i and q is the parameter that characterizes the degree of nonextensivity of the system. The classical entropy is recovered in the limit q 1. In Tsallis’ statistics, the particle number of species i per volume can be written as follows: →

+∞ q gi 3 N q ni = 3 d p i , (23) (2πh¯) ∞ Z− N q where i is the generalized occupation number, gi is the degeneracy of species i and h¯ is the Planck reduced constant. The generalized occupation number for fermions in Tsallis framework is given by [18] 1 N q = , (24) i β (Ei µi) e2 q − + 1 − 5

x 1/(q 1) where e2 q [1 + (q 1)x] − . Further, β = 1/kBT, µi and Ei are the chemical potential and particle energy of species i, respectively.− ≡ − Expanding Eq. (24) up to the first order of (q 1), we obtain [19] − 1 (q 1) (β(E µ ))2eβ (Ei µi) N q = + i i − . (25) i β (E µ ) − − 2 e i i + 1 2 β (Ei µi) − e − + 1

β (Ei µi)
Considering the case in which kBT < (Ei µi), then e 1. Thus, Eq. (25) becomes − − ≫ q 1 q β (Ei µi) ( ) 2 β (Ei µi) N = e− − + − (β(Ei µi)) e− − . (26) i 2 − 2 2 In this way, assuming the energy for non-relativistic particles as Ei = mic + p /2mi and using Eq. (26) in Eq. (23), the generalized particle number density take the form

∞ ∞ 2πg 2 + 2 2πg q 1 2 + 2 2πg q 1 q i β (µi mic ) 2 α p i ( ) 2 β (µi mic ) 2 2 2 α p i ( ) 2 ni = 3 e − dpp e− + 3 − β e − (mic µi) dpp e− + 3 − β (2πh¯) ∞ (2πh¯) 2 − ∞ (2πh¯) 2 Z− Z− 2 +∞ 2 g q 1 2 π +∞ 2 2 β (µi mic ) 4 α p i ( ) 2 β (µi mic ) 6 α p mi(mic µi)e − dpp e− + 3 − β e − 2 dpp e− , (27) × − ∞ (2πh¯) 2 2m ∞ Z− i Z− where α β/2mi and c is the speed of light in vacuum. ≡ Note that in the above expression we have four Gaussian integrals of type I2n. It is at this point that we will apply our method. To solve these Gaussian integrals, we will use the Theorem 1. The results are:

+∞ 2 α p2 3/2 dpp e− = γ2(2mikBT ) √π , (28) Z ∞ −+∞ 4 α p2 5/2 dpp e− = γ4(2mikBT ) √π , (29) Z ∞ −+∞ 6 α p2 7/2 dpp e− = γ6(2mikBT ) √π . (30) ∞ Z− We will use the Definition 1 to compute the coefficients γ2n, thus we obtain γ2 = 1/2, γ4 = 3/4 and γ6 = 15/8. Finally, we can find the generalized particle number density. After a little algebra, we get [19]

2 3/2 µ m c2 2 2 q m k T i− i (q 1) mic µi mic µi 15 i B kBT ni = gi 2 e 1 + −2 − + 3 − + . (31) 2πh¯ ( " kBT kBT 4 #)       Notice that the usual particle number density is recovered when q = 1. We chose to apply our method to determine the particle number density in Tsallis framework because it arises integrals in which the n index assumes three different values that facilitated the application. This enables us to visualize from another perspective how the Gaussian integral evolve as n grows. As mentioned before, this method is viable to apply in mathematical and physical problems where the n index is small, since we would need to perform derivatives successive to obtain the result. On the other hand, the derivative order is half of the value of 2n index, for example, in Eq. (27) the integral containing the term 6 p , we need to derive only three times, making the practical and efficient method in the resolution of integrals of type I2n. Even that n is large, it was demonstrated that the method also works for all n N. ∈

B. Special functions

An immediate application of Theorem 1 is verified in the gamma and factorial functions. Below is the definition of the special functions we will consider from now on. Definition 3. (Euler, 1730) Let t R, and t > 0, the gamma function is defined by [20, 21] ∈ +∞ t 1 x Γ(t) = x − e− dx . (32) Z0 Using t = p + 1, and integrating by parts, we obtain the following recurrence relation: Γ(p + 1) = pΓ(p) . (33) 6

Definition 4. For m > 1, the factorial function is defined by [20, 21] − +∞ m x m! = x e− dx . (34) Z0 We begin with the gamma function given by Definition 3, setting t = n + 1/2, that is

2n + 1 +∞ 2n 1 Γ = x 2− e x dx . (35) 2 −   Z0 Taking x = r2,

+∞ 2n + 1 2n r2 Γ = r e− dr . (36) 2 ∞   Z− Note that we can apply Theorem 1, with α = 1, in the above expression to obtain

2n + 1 Γ = γ √π . (37) 2 2n  

Having presented (37), it is convenient to introduce a recurrence relation for the parameter γ2n.

Lemma 2. If γ2n is given by Definition 1, then 2n + 1 γ = γ . (38) 2(n+1) 2 2n Proof. Using (37), and the succeeding term, given by

2n + 1 Γ + 1 = γ √π , (39) 2 2(n+1)   we can insert these expressions in Eq. (33), defining p = n + 1/2. Therefore

2n + 1 γ = γ . (40) 2(n+1) 2 2n

The Theorem 1 ensure us that the parameter γ2n is given by Definition 1. Based on this, it is straightforward to show that the following equality is valid:

∞ n ∞ 1 x n n 2 γ2n = ( 1) x − , (41) ∑ n! ∑ − n n=0 n=0   x where y is the binomial coefficient. From the above expression, we can conclude that
n √π γ2n = ( 1) , (42) − 1 n ! − 2 − which yields the following result:

1 √π Γ n = ( 1)n , (43) 2 − − γ   2n where we use the identity m! = Γ(m + 1). We now let us apply Theorem 1 in the factorial function given by Definition 4, putting m = n + 1/2. Then we have

+∞ 2n + 1 2n+1 x ! = x 2 e dx . (44) 2 −   Z0 7

Replacing again x = r2, we find

+∞ 2n + 1 2(n+1) r2 ! = r e− dr . (45) 2 ∞   Z− Hence using Theorem 1, we obtain

2n + 1 ! = γ √π , (46) 2 2(n+1)   and making use of Lemma 2, we get

2n + 1 2n + 1 ! = γ √π . (47) 2 2 2n  

C. Fractional derivative: a possibility

Here we intend to present an interesting result that may be useful in fractional calculus. Initially let us consider the following definitions: Definition 5. Let Re(c) > Re(b) > 0, the hypergeometric function is defined by [22]

1 Γ(c) b 1 c b 1 a 2F1 (a,b;c;z) = τ − (1 τ) − − (1 zτ)− dτ , (48) Γ(b)Γ(c b) 0 − − − Z where Re(x) is real part of x. Definition 6. The fractional derivative according to Riemann-Liouville is defined by [23]

m t α 1 d f (τ) D f (t) = m α m 1 dτ , (49) Γ(m α) dt 0 (t τ) + − Z − − where α is a complex number such that Re(α) > 0 and yet shall is in interval m 1 < Re(α) m with m N. − ≤ ∈ Definition 7. The fractional derivative according to Caputo is given by [23]

t m α 1 m α 1 d D f (t) = (t τ) − − m f (τ)dτ , (50) Γ(m α) 0 − dτ − Z being Re(α) > 0 such that m 1 < Re(α) m with m N. − ≤ ∈ By using f (t)= tθ with θ > 1 in Definition 6, we obtain −

α θ Γ(θ + 1) θ α D t = t − . (51) Γ(θ α + 1) − For the case where α = θ, we have dαtα = α! . (52) dtα Let us now show that the above result can be represented by the k-th derivative of a simple function applied to a point. Lemma 3. Let k N, then ∈ k 1 d (1 x)− −k = k! . (53) dx x=0

8

Proof. In fact,

dk(1 x) 1 ∞ n − − = Γ(k + 1) xn k , dxk ∑ k − n=k   k! = , (54) (1 x)k+1 − where ∞ n n k 1 x − = . (55) ∑ k (1 x)k+1 n=k   − Hence, at x = 0, we obtain

k 1 d (1 x)− −k = k! . (56) dx x=0

Suppose that k can take the form k = n + 1/2, the Lemma 3 may be written as follow:

2n+1 1 d 2 (1 x)− 2n + 1 2n−+1 = ! . (57) dx 2 2 x=0  

Let us now analyze the left-hand side of the above equation. To accomplish this, we will use the Definitions 6 and 7 to evaluate the fractional derivatives. Thus, using the Definition 6 with α = n + 1/2, we find the following result:

n 1 2n+1 1 x− − 2 1 D 2 (1 x)− = 2F1 1,1;m n + ;x , (58) − Γ 1 n − 2 2 −  
where 2F1 (a,b;c;z) is the hypergeometric function given by the Definition 5. This result applied at the point x = 0 diverge. Now using the Definition 7 with α = n + 1/2 on the left-hand side of Eq. (57), we obtain

2n+1 1 Γ(m + 1) m n 1 1 D 2 (1 x)− = x − − 2 2F1 m + 1,1;m n + ;x . (59) − Γ m n + 1 − 2 − 2   Note that at the point x = 0, the above result goes to zero, since
m > n + 1/2. We show that both Riemann-Liouville and Caputo 1 definitions are flimsy when computed the fractional derivatives of the function f (x) = (1 x)− at the point x = 0. On the other hand, this problem can be avoided considering the result presented in Eq. (47). Hence, we− may write

2n+1 1 2n + 1 D 2 (1 x)− = γ2n√π . (60) − x=0 2

IV. CONCLUSIONS

It was presented a new method for solving the Gaussian integral inspired by expansion in Taylor series of a simple function, namely f (x) = (1 x) 1/2. Introducing a parameter with a n index dependence, we have found a general solution for this type − − of integral being able to use in any situation, since the Gaussian integral is of type I2n. We have demonstrated that the parameter γ2n represents the Taylor series coefficients of the function f (x), a result newsworthy. The reliability of the method is guaranteed through the proof of some Theorems, which were proved by mathematical induction. To check the functionality of the method, we have obtained the particle number density in the Tsallis framework. Besides, we have presented solutions for the gamma function of the form Γ(1/2 n), and the factorial function of a kind (n + 1/2)!, in terms of the parameter γ2n. We have also shown that the method is useful± in fractional calculus, for example, the definitions of fractional derivative according to Riemann- 1 Liouville and Caputo are flimsy when evaluated the fractional derivatives of the function f (x) = (1 x)− at the point x = 0, whereas using our method we find the result as shows Eq. (60), showing again the functionality and efficiency− of the method. In conclusion, we believe in the relevance of the present work because it reveals a new way of solving the Gaussian integral: one of the most famous integrals of the exact sciences. 9

Acknowledgements

The authors are very grateful to R. C. Duarte, N. S. Almeida, T. Dumelow and A. M. Filho by the helpful discussions and the Brazilian agency CAPES for financial support.

[1] P. J. Davis, The American Mathematical Monthly 66, 847-869 (1959). [2] D. Gronau, Teaching Mathematics and Cumputer Science 1, 43-53 (2003). [3] S. Stahl, Mathematics Magazine 79, 2, 96-113 (2006). [4] M. R. Spiegel, J. Schiller, and R. A. Srinivasan, Probability And Statistics (Mc Graw Hill, 2001). [5] W. Greiner, Quantum Mechanics An Introduction (Springer, Berlin, 1990). [6] J.J. Sakurai, Modern Quantum Mechanics (AddisonWesley, Redwood City 1985). [7] W. Greiner, Thermodynamics and Statistical Mechanics (Springer Verlag, New York, 1995). [8] R. K. Pathria, Statistical Mechanics (Butterworth-Heinemann, Oxford, 1996). [9] S. R. A. Salinas, Introduction to Statistical Physics (Springer, New York, 2001). [10] J. K. F. Sturm, Cours d’Analyse de l´ecole polytechnique (Paris, Mallet-Bachelier, 1857). [11] S. M. Stigler, Statistical Science 1, 359378, (1986). [12] P. S. Laplace, Th´eorie Analytiques des Probabilit´es. Paris, Courcier 1812, reprinted (Bruxelles, Culture et Civilisation, 1967). [13] K. T. Conrad, The Gaussian Integral. (2013). Online at https://www.semanticscholar.org/paper/THE-GAUSSIAN-INTEGRAL-Conrad/4687538f80e333c175691d627dc1254eef3605f8 [14] S. M. Hernandez, Termodin`amica i Mec`anica estad´ıstica (Libre Lulu, Breda, 2015). [15] E. W. Weisstein, Gaussian integral. From MathWorld-A Wolfram Web Resource. http://mathworld.wolfram.com/GaussianIntegral.html [16] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, New York, 2005). [17] C. Tsallis, Possible generalization of Boltzmann-Gibbs statistics 52, 479-487 (1988). [18] K.-M. Shen, B.-W. Zhang, and E.-K. Wang, Physica A: Statistical Mechanics and its Applications 487, 215 (2017) [19] M. E. Pessah, D. F. Torres and H. Vucetich, Statistical mechanics and the description of the early universe. (I). Foundations for a slightly non-extensive cosmology 297, 164-200 (2001). [20] M. L. Boas, Mathematical methods in the physical sciences. (John Wiley & Sons, New Jersey, 2006). [21] K. F. Riley, M. P. Hobson, and S. J. Bence, Mathematical Methods for Physics and Engineering (Cambridge University Press, Cambridge, 2006). [22] E. W. Weisstein, Hypergeometric Function. From MathWorld-A Wolfram Web Resource. https://mathworld.wolfram.com/HypergeometricFunction.html [23] E. C. Oliveira and J. A. T. Machado, A review of definitions for fractional derivatives and integral, Mathematical Problems in Engineering 2014, (2014).