﻿功能:求伽马函数值

格式:BgGamma(x)

说明:x为数据或者符号变量的值,函数执行完毕返回对应的值.注意,这里f值的绝对值必须小于1E10.

原理:

这里设x>0
1、
Log[Γ(x)]≈Log(2*π)/2+(x-0.5)*Log(x)-x+1/(12*x)-1/(360*x^3)+1/(1260*x^5)-1/(1680*x^7)+1/(1188*x^9)-691/(360360*x^11)+1/(156*x^13)-3617/(122400x^15)+43867/(244188*x^17)-174611/(125400*x^19)+77683/(5796*x^21)-236364091/(1506960*x^23)+……

2、
Γ(1-x)*Γ(x)=π/sin(π*x)

参考:{Gamma函数<http://mathworld.wolfram.com/GammaFunction.html>}

例子:

BgGamma(3.456789)//回车后得到如下结果
ans =
"3.1696303952201826147892449721669169053011984092325722586758995041452031346627779310891273334896150153975784842086850367751890754061424887152354518450535131598050568728763924025325365163719535575807066094482024935072387624655361167272374670474858454927708996737014250951619122869850659421275439148129408"