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

问题描述

最近，我在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。对不起，我的英语不好。希望你能理解。.

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）