File:Lua Gamma Function in Chinese Wiki.svg
页面内容 不 支持 其他语言。
![File:Lua Gamma Function in Chinese Wiki.svg](https://upload.wikimedia.org/wikipedia/commons/thumb/c/cc/Lua_Gamma_Function_in_Chinese_Wiki.svg/539px-Lua_Gamma_Function_in_Chinese_Wiki.svg.png)
此SVG文 件 的 PNG预览的 大小 :539 × 482像 素 。 其他分 辨 率 :268 × 240像 素 | 537 × 480像 素 | 859 × 768像 素 | 1,145 × 1,024像 素 | 2,290 × 2,048像 素 。
摘要
描述Lua Gamma Function in Chinese Wiki.svg |
|
|
|
A2569875 |
许可协议
![w:zh:知识共享](https://upload.wikimedia.org/wikipedia/commons/thumb/7/79/CC_some_rights_reserved.svg/90px-CC_some_rights_reserved.svg.png)
![署名](https://upload.wikimedia.org/wikipedia/commons/thumb/1/11/Cc-by_new_white.svg/24px-Cc-by_new_white.svg.png)
![相同方式共享](https://upload.wikimedia.org/wikipedia/commons/thumb/d/df/Cc-sa_white.svg/24px-Cc-sa_white.svg.png)
- 您可以
自由 地 :共 享 – 复制、发行并传播本 作品 修 改 –改 编作品
惟 须遵守 下 列 条件 :署名 – 您必须对作品 进行署名 ,提供 授权条 款的链接,并说明 是 否 对原始 内容 进行了 更改 。您可以用任 何 合理 的 方式 来 署名 ,但 不 得 以任何 方式 表明 许可人 认可您或您的使用 。相 同 方式 共 享 – 如果您再混合 、转换或 者 基 于本作品 进行创作,您必须以与原 先 许可协议相 同 或 相 兼 容 的 许可协议分 发您贡献的 作品 。
code
local Reciprocal_gamma_coeff = {1,0.577215664901532860607,-0.655878071520253881077,-0.0420026350340952355290,0.166538611382291489502,-0.0421977345555443367482,-0.00962197152787697356211,0.00721894324666309954240,-0.00116516759185906511211,-0.000215241674114950972816,0.000128050282388116186153,-0.0000201348547807882386557,-1.25049348214267065735e-6,1.13302723198169588237e-6,-2.05633841697760710345e-7,6.11609510448141581786e-9,5.00200764446922293006e-9,-1.18127457048702014459e-9,1.04342671169110051049e-10,7.78226343990507125405e-12,-3.69680561864220570819e-12,5.10037028745447597902e-13,-2.05832605356650678322e-14,-5.34812253942301798237e-15,1.22677862823826079016e-15,-1.18125930169745876951e-16,1.18669225475160033258e-18,1.41238065531803178156e-18,-2.29874568443537020659e-19,1.71440632192733743338e-20}
--https://oeis.org/A001163 、 https://oeis.org/A001164
local stirling_series_coeff = {1,0.0833333333333333333333333,0.00347222222222222222222222,-0.00268132716049382716049383,-0.000229472093621399176954733,0.000784039221720066627474035,0.0000697281375836585777429399,-0.000592166437353693882864836,-0.0000517179090826059219337058,0.000839498720672087279993358,0.0000720489541602001055908572,-0.00191443849856547752650090,-0.000162516262783915816898635,0.00640336283380806979482364,0.000540164767892604515180468,-0.0295278809456991205054407,-0.00248174360026499773091566,0.179540117061234856107699,0.0150561130400264244123842,-1.39180109326533748139915,-0.116546276599463200850734}
function p._gamma_high_imag(cal_z)
local z = to_number(cal_z)
if z ~= nil and math_lib.abs(math_lib.nonRealPart(z)) > 1 then
local inv_z = math_lib.inverse(z)
return math_lib.sqrt((math_lib.pi * 2) * inv_z) * math_lib.pow(z * math_lib.exp(-1) *
math_lib.sqrt( (z * math_lib.sinh(inv_z) ) + math_lib.inverse(to_number(810) * z * z * z * z * z * z) ),z)
end
return nil
end
function p._gamma_morethen_lua_int(cal_z)
local z = to_number(cal_z) - to_number(1)
local lua_int_term = 18.1169 --FindRoot[ Factorial[ x ] == 2 ^ 53, {x, 20} ]
if math_lib.abs(z) > (lua_int_term - 1) or (math_lib.re(z) < 0 and math_lib.abs(math_lib.nonRealPart(z)) > 1 ) then
local sum = 1
for i = 1, #stirling_series_coeff - 1 do
local a, n = to_number(z), tonumber(i) local y, k, f = to_number(1), n, to_number(a)
while k ~= 0 do
if k % 2 == 1 then y = y * f end
k = math.floor(k / 2); f = f * f
end
sum = sum + stirling_series_coeff[i + 1] * math_lib.inverse(y)
end
return math_lib.sqrt( (2 * math.pi) * z ) * math_lib.pow( z * math.exp(-1), z ) * sum
end
return nil
end
function p._gamma_abs_less1(cal_z)
local z = to_number(cal_z)
if math_lib.abs(z) <=1.001 then
if math_lib.abs(math_lib.nonRealPart(z)) < 1e-14 and ( (math.abs(math_lib.re(z) - 1) < 1e-14) or (math.abs(math_lib.re(z) - 2) < 1e-14) ) then return to_number(1)end
return math_lib.inverse(p._recigamma_abs_less1(z))
end
return nil
end
function p._recigamma_abs_less1(z)
local result = to_number(0)
for i=1,#Reciprocal_gamma_coeff do
result = result + Reciprocal_gamma_coeff[i] * math_lib.pow(z,i)
end
return result
end
function p._gamma(cal_z)
local z = to_number(cal_z)
if math_lib.abs(math_lib.nonRealPart(z)) < 1e-14 and ((math_lib.re(z) < 0 or math.abs(math_lib.re(z)) < 1e-14)
and math.abs(math.floor(math_lib.re(z)) - math_lib.re(z)) < 1e-14 ) then return tonumber("nan") end
local pre_result = p._gamma_morethen_lua_int(z) or p._gamma_high_imag(z) or p._gamma_abs_less1(z)
if pre_result then return pre_result end
local real_check = math_lib.re(z)
local loop_count = math.floor(real_check)
local start_number, zero_flag = z - loop_count, false
if math_lib.abs(start_number) <= 1e-14 then start_number = to_number(1);zero_flag = true end
local result = math_lib.inverse(p._recigamma_abs_less1(start_number))
if math_lib.abs(math_lib.nonRealPart(z)) < 1e-14 and ((math_lib.re(z) > 1e-14 )and math.abs(math.floor(math_lib.re(z)) - math_lib.re(z)) < 1e-14 ) then result = to_number(1) end
local j = to_number(start_number)
for i=1,math.abs(loop_count) do
if loop_count > 0 then result = result * j else result = result * math_lib.inverse(j-1) end
if zero_flag==true and loop_count > 0 then zero_flag=false else if loop_count > 0 then j = j + 1 else j = j - 1 end end
end
if math_lib.abs(math_lib.nonRealPart(z)) < 1e-14 and ((math_lib.re(z) > 1e-14 )and math.abs(math.floor(math_lib.re(z)) - math_lib.re(z)) < 1e-14 ) then return math_lib.floor(result) end
return result
end
Reference
- ↑ Wrench, J.W. (1968). Concerning two series for the gamma function. Mathematics of Computation, 22, 617–626. and
Wrench, J.W. (1973). Erratum: Concerning two series for the gamma function. Mathematics of Computation, 27, 681–682. - ↑ Viktor T. Toth (2006). "Programmable Calculators: Calculators and the Gamma Function". Archived from the original on 2007-02-23.
- ↑ NIST Digital Library of Mathematical Functions.
说明
此文件 中 描述的 项目
描繪內容
18 11 2018
image/svg+xml
文 件 历史
缩 | 备注 | ||||
---|---|---|---|---|---|
2018 | ![]() | 539 × 482(10 KB) | A2569875 | User created page with UploadWizard |
文 件 用途
元 数 据
此文
如果此文
宽度 | 538.58px |
---|---|
481.89px |