The generalization of elementary functions such as exponential function in q-exponential has been proposed since long time ago . However, more recently, another type of q-exponential has been introduced in the environment of the nonadditive statistical mechanics , where the theory has been proposed as a generalization of the standard Boltzmann-Gibbs statistical mechanics for the class of system that is non-ergotic. Many other q-representations were proposed such as the q-algebra , the q-Fourier transform  , the q-Dirac’s delta function  , the q-Gaussian    and so on. Moreover, combinations of complex q-exponentials were used to make a new representation of the classical nonlinear oscillator , where a combination of complex q-exponentials have been of interest also for recent q-generalized Schrödinger and Dirac equations and other field equations . In the non-additive statistical mechanics the entropy is described by the so called Tsallis entropy , previously introduced by Havrdat-Charvát  and later on by Daróczy , parameterized by an exponent q. The distribution that maximizes the Tsallis entropy under the constraint of the fixed energy generalizes the Boltzmann distribution through a generalized q-exponential function. Furthermore, in the nonadditive statistical mechanics, when the limit of the parameter of non-extensiveness q is , one recovers the standard exponential (and thus, the Gaussian or Boltzmann distribution). Many other classical equations (e.g. the Schrödinger equation) were generalized, replacing the classical exponential with the q-exponential. One example of non-ergodic systems that have been studied is system with long range interactions   .
In addition, there is another definition of q-exponential given in combinatorial mathematics, where q-exponential is a q-analog of the standard exponential function namely, the eigenfunction of a q-derivative. There are many q-derivatives, for example, the classical q-derivative, the Askey-Wilson operator and so on . Therefore, unlike the classical exponentials, q-exponentials are not unique. Moreover, it has been constructed rational expansions to by Németh and Newman in Ref. .
A special function very important is the gamma function , , . This is an extension for the complex plane of the factorial function , for , with its argument shifted down by 1, to real and complex numbers. Thus, the gamma function is defined for all complex numbers except for non positive integers. It raises as component in various probability distributions, and as such it is applicable in the fields of probability, statistics and statistical physics, as well as in combinatorics     . For instance, the well known Stirling formula depends on the factorial function as 
Furthermore, this function if relate with the Shannon’s entropy that gives an approach to the information measure
The gamma function is extended by analytic continuation to all complex numbers except the non-positive integers (where the function has simple poles), yielding a meromorphic function which we call the gamma function. It has no zeroes and hence, the reciprocal gamma function is a holomorphic function. In fact, gamma function corresponds to Mellin transform  of the negative exponential function: .
The generalized gamma function can be applied to nonlinear systems (NL) class obeyed by the nonlinear equations. This has become an important subject in the mathematical physics in recent years due to their ability to explain several complex behaviors in the nature  . Many areas of physics such as plasma physics and non-equilibrium have benefited from the study of NL equations . However, it is well known that, in general, nonlinear problems are impossible to solve analytically. The essential difference is that linear systems can be broken down into parts while the nonlinear not. When the parts of a system interfere, either cooperate or compete. The interactions are nonlinear and hence, the principle of superposition fails . Within the realm of physics, the nonlinearity is vital for many areas, for instance, for the superconductivity and Josephson junctions arrays .
In this work, we propose the representation of some special functions such as the gamma function in terms of q-exponentials functions of nonadditive statistical mechanics. The plan of this paper is the following. In Section 2, we represent the generalized gamma function expressed in terms of q-exponential functions. In Section 3, we use another representation for the q-exponential in the literature to represent the generalized gamma function. In the last section, Section 4, we present our conclusions and final remarks.
2. Generalized Gamma Function
The generalized gamma function, defined using the q-exponential distribution, is defined by the integral
where and and . In the limit , we have and for and for .
Since , when x is positive and , we can write
and, for the integral is limited for .
By making x fixed and decreasing the value of integral increases monotonically, i.e.
presents the properties . The q-exponential function is a deformation of the standard exponential function using the real parameter q  
The inverse of the q-exponential function is the q-logarithm function, , defined as
The formulae in Equation (6) hold true only for . In this interval x and q are mathematically independent. The complete definitions can be given in two equivalent ways: one either changes the full expression of the q-deformed functions appropriately and uses the interval only as in the equation above or one can consider a unique expression and change the deformation parameter interval as 
The parameter q represents the degree of non-additivity. Thus, solving the integral, Equation (3) and using the definition of q-exponential, we obtain the generalized gamma function given as
where for and , . Therefore, using the standard factorial function, we have the recurrence relation obeyed for the generalized gamma function given by
Consequently, we can obtain the expression for the q-factorial, as
We also obtain the generalized incomplete gamma functions
with , where
In addition, we have the following generalized functions that if relate with the incomplete generalized gamma function
that is the generalized complementary error function, besides
where we define the generalized exponential integral function as
In Figure 1 we plot the generalized gamma function for q value and for standard gamma function (that corresponds to case. The graphic of changes a lot with q value as showed in Figure 2. This is due to the fact the q-exponential function to represent a family of functions (one for each value of q within interval where the case (the standard exponential function, ) corresponds to only one function of these class of q-exponential functions. We get that for the particular case , the q-gamma function exhibits a behavior nearest to standard exponential than q-exponentials defined for other q values.
By using the definition of q-exponential, we obtain the generalized incomplete q-gamma function given as
Figure 1. Plot of the generalized gamma function for q value and for standard gamma function (that corresponds to case.
Figure 2. Plot of the generalized gamma function for q values (red-dashed-line) and (black-solid-line). The standard gamma function case corresponds to case.
In addition, we also obtain the generalized exponential integral function as
Finally, we have the generalized logarithm integral function , given by
Consequently, we obtain given as
3. Another Representation
Another type of q-exponential is defined using the following series expansion 
where is the q-factorial 
where we use the q-bracket given by
In addition, we have the q-Pochhammer symbol given as  
The graphic of , given by Equation (23), for q value of is given in Figure 3. The black-solid-line corresponds to series in Equation (23) until . The red-dashed-line correspond to sum until . So, for , the behavior of graphic does not change a lot (for ).
The generalized gamma function given by the integral
In Figure 4, we plot the generalized gamma function for and q value as . The behavior of the graphic changes with q as in the before case. However, we have here, and and . Remembering that here, this q-exponential is not that defined in the nonadditive statistical mechanics. Therefore, we obtain the integral Equation (27) as (Figure 4)
Figure 3. Plot of , given by Equation (23), for q value of . The black-solid-line corresponds to series in Equation (23) until . The red-dashed-line correspond to sum until .
Figure 4. Graphic of generalized gamma function for integer z, and . The sum in series given by Equation (23) is for . The axis has been rescaled as .
where the integration contour C in the integral above is represented in Figure 5; and . The integral above can be solved using the residue technique as
where , , where we have taken the limit and thus, from the above integral
which cancels for . For other z values the integral of the generalized gamma function diverges. However, for the incomplete gamma function
The integral of the incomplete generalized gamma function converges values. Because of the complicated behavior of the derivative of the q-exponential function, the generalized gamma function here does not obey the property obeyed by the standard gamma function.
For the generalized exponential integral function, we obtain
where . However for the integral above diverges. The same behavior we have for the generalized integrals Equation (16).
4. Final Remarks
In summary, we have proposed a generalization of the complete and uncomplete gamma functions using q-exponential functions. This subject can be of interest
Figure 5. Contour for the calculation of the integral in the Equation (29).
in many branches of physics due to the many phenomena associated with the nonlinearity. Recently, a connection between generalized functions of the nonadditive statistical mechanics framework and nonlinear equations has been investigated. In literature, were proposed other deformations of the -function in terms of deformed factorials , where the q-factorial is . These definitions are different from the proposed here.
This work was partially supported by CNPq (National Council for Scientific and Technological Development.