朱莉娅不完全伽玛函数的拉普拉斯逆变换

时间:2020-11-02 07:06:21

标签: python julia inverse mpmath

最近,我在julia有一个与反拉普拉斯有关的问题。

这是我在茱莉亚写的代码。

using SpecialFunctions
using InverseLaplace

value = Talbot(u -> gamma_inc(1,u,2)[2], 80)
value(2)

此代码会产生一些错误消息,例如

MethodError: no method matching gamma_inc(::Float64, ::BigFloat, ::Int64, ::String)

实际上,我已经可以基于mpmath在python中计算出这种功能了。

import mpmath as mp

laplace = lambda u: mp.gammainc(1,u)
result = mp.invertlaplace(laplace, t = 2, method= 'talbot',degree = 80)

print(result)

结果是

5.6730679553874e-43

我想在茱莉亚上解决这个问题,有人可以给我任何建议吗?

p.s。对不起,我的英语不好。希望你能理解。.

0 个答案:

没有答案