gamma.f90

Path: src/gamma.f90
Last Update: Fri Feb 27 14:51:08 +0900 2009

gamma

Authors:Shin-ichi Takehiro, Youhei SASAKI
Version:$Id: gamma.f90,v 1.5 2009-02-27 05:51:08 uwabami Exp $

概要

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

関数

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

Methods

factrl   gammaln  

Public Instance methods

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
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

[Validate]