gamma

Version:

$Id: gamma.f90 590 2013-08-19 08:48:21Z uwabami $

概要

このファイルは gamma 関数に関する計算用の関数を提供する

関数

factrl :	与えられた整数に対する階乗(ガンマ関数)の値を返す
gammaln :	与えられた正の実数に対する階乗の対数(ガンマ関数の対数)の値を返す

Methods

factrl gammaln

Public Instance methods

factrl( n ) result(factrl)

*Function* :
`factrl` :	`real(8)`
`n` :	`integer, intent(IN)`

Returns the value n! as a floating point number

[Source]

function factrl(n)
  !
  ! Returns the value n! as a floating point number
  !
  integer, intent(IN) :: n
  real(8)             :: factrl
  real(8)             :: gammaln
  external gammaln

  factrl=exp(gammaln(n+1.0d0))

  return

end function factrl

gammaln( xx ) result(gammaln)

*Function* :
`gammaln` :	`real(8)`
`xx` :	`real(8)`

Returns the vaule ln[Γ(x)] for x >= 0

[Source]

function gammaln(xx)
  !
  ! Returns the vaule ln[Γ(x)] for x >= 0
  !
  real(8) ::  gammaln, xx
  real(8) ::  dlgamma
  external dlgamma

  gammaln = dlgamma(xx)

  return
end function gammaln