文件名称:广义指数积分:计算正实参数 a 和参数 x 的广义指数积分 E_a(x)-matlab开发
文件大小:8KB
文件格式:ZIP
更新时间:2024-06-21 02:11:49
matlab
此函数计算正实参数 a 和参数 x 的广义指数积分 E_a(x)。 将其称为 y=genexpint(a,x) 或 y=genexpint(a,x,expscale)。 如果可选的第三个输入参数 expscale 设置为 true,则输出为 exp(x)*E_a(x),这对于大 x 是有限的,其中 exp(x) 溢出和 E_a(x) 下溢。 该代码使用来自 SLATEC 库的 FORTRAN 函数 DGAMIC 的 MATLAB 转换。 DGAMIC 使用 Gautschi 算法(ACM Trans. Math. Soft. 5(4) pp 466-481, 1979)计算负实参数的上不完全伽马函数。 对于 x>1,Legendre 连分数用于计算函数G(1-a)=exp(x)*E_a(x)。 对于 x<=1,广义指数积分由关系式获得E_a(x)=x^(a-1)*Gamma(1-a,x)。
【文件预览】:
genexpint.zip