The gamma function is defined as:
Use the LNGAMMA function to obtain the natural logarithm of the gamma function when there is a possibility of overflow.
Plot the gamma function over the range 0 to 0.5:
; Generate the numbers to use in the GAMMA function.
Z = FINDGEN(50)/100.
; Plot the GAMMA of each of these numbers.
YTITLE='Gamma of Numbers', COLOR='red', $
TITLE="Plot of a Number vs Its GAMMA", $
NAME='Number to the Power of 1')
; Plot the GAMMA of the numbers vs the number
; raised to one of several exponents.
NAME='Number to the Power of 2', /OVERPLOT)
NAME='Number to the Power of 3', /OVERPLOT)
NAME='Number to the Power of 4', /OVERPLOT)
NAME='Number to the Power of 5', /OVERPLOT)
L = LEGEND(POSITION=[0.4,90], /DATA)
Result = GAMMA(Z)
If Z is double-precision, the result is double-precision (either double or double complex), otherwise the result is single-precision (either float or complex).
The expression for which the gamma function will be evaluated. Z may be complex.
Thread Pool Keywords
This routine is written to make use of IDL’s thread pool, which can increase execution speed on systems with multiple CPUs. The values stored in the !CPU system variable control whether IDL uses the thread pool for a given computation. In addition, you can use the thread pool keywords TPOOL_MAX_ELTS, TPOOL_MIN_ELTS, and TPOOL_NOTHREAD to override the defaults established by !CPU for a single invocation of this routine. See Thread Pool Keywords for details.
Z argument accepts complex input