sidef-scripts/Math/gamma_function.sf

#!/usr/bin/ruby

#
## https://rosettacode.org/wiki/Gamma_function
#

define ℯ = Num.e
define τ = Num.tau
 
func Γ(t, r=50) {
    t < r ? (__FUNC__(t + 1) / t)
           : (sqrt(τ*t) * pow(t/ℯ + 1/(12*ℯ*t), t) / t)
}
 
for i in (1..10) {
    say ("%.14f" % Γ(i/3))
}