RosettaCodeData/Task/Gamma-function/Julia/gamma-function-2.julia

using QuadGK
gammaquad(t::Float64) = first(quadgk(x -> x ^ (t - 1) * exp(-x), zero(t), Inf, reltol = 100eps(t)))
@show gammaquad(1.0)
June 2018 Update 2018-06-22 20:57:24 +00:00			`using QuadGK`
			`gammaquad(t::Float64) = first(quadgk(x -> x ^ (t - 1) * exp(-x), zero(t), Inf, reltol = 100eps(t)))`
			`@show gammaquad(1.0)`