calc_gamma_cf Function

private function calc_gamma_cf(a, x) result(gamma_cf)

solves the normalized upper gamma function $\gamma(a,x) / \Gamma(a)$ using a continued fractions method (modified Lentz Algorithm)

input: a = gamma function argument, x = lower limit of integration

output: normalized upper gamma function: $\gamma(a, x) / \Gamma(a)$

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	a
real(kind=wp),	intent(in)			::	x

Return Value real(kind=wp)